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PREFACE. 



The present volume contains all that Professor Boole 
wrote for the purpose of enlarging his Treatise on Differential 
Equations, Had he lived to publish the second edition he 
would doubtless have Incorporated his more recent investi- 
gations with the original work, and it is therefore necessary 
to explain why another plan has been adopted. 

In some cases Professor Boole had indicated that certain 
portions of the original work were to be omitted and their 
places supplied from the manuscripts; but on examination 
it appeared that in subsequent passages of the work there 
were references and allusions to the portions thus marked to 
be omitted which would not apply to the substituted matter. 
Thus in attempting to carry out the directions it would 
have been necessary to accept the responsibility of making 
many alterations, and consequently to incur the risk of fail- 
ing in the attempt to improve the original form. 

Moreover the Treatise had been for some time out of 
print, and the long delay which must have been caused by 
the labour of reconstruction would have produced serious 
inconvenience to students at Cambridge and elsewhere. Pro- 
fessor Boole himself was always especially anxious to consult 
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the advantage of students, and those who had the charge 
of his manuscripts were naturally inclined to adopt a course 
of which they believed he would himself have approved. 

The design of reconstructing the Treatise was therefore 
abandoned; and it was resolved that the original volume 
should be reprinted, and that the manuscripts should be 
collected and published separately. This plan has the ob- 
vious recommendation of enabling those who are already 
familiar with the original work to turn their attention 
readily to the new investigations. It will be seen that 
many of the Chapters of the present volume may be re- 
garded as independent essays or memoirs which lose nothing 
by being separated from the other volume; and indeed no 
indications had been left by Professor Boole of the place 
which such Chapters were to occupy in the enlarged edition. 

I have printed all the unpublished matter relating to 
Differential Equations which I found among Professor Boole's 
papers. In a few cases it will be seen that an investigation 
is incomplete; such investigations have however been in- 
cluded in the volume, because I was unwilling that anything 
should be lost which so great a mathematician had written 
on a subject he had long and carefully studied. 

I trust that no serious error will be found in the volume, 
and that any faults which may be detected will be excused 
on account of the nature and difficulty of the task that had 
to be performed. Many of the manuscripts had not been 
finally revised; some of them were very obscure and had 
to be carefully and laboriously copied for the press. In 
general the equations were not numbered, and thus only 
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blanks occurred in place of references; this circumstance 
often caused great trouble and perplexity : I hope however 
that a satisfactory result has been finally attained. 

I may state for the benefit of those who are conversant 
with the first edition of the original work that the theo- 
rem which in the present volume is cited as contained in 
Chap. II. Art. 1 will be found in Chap. IV. Art. 2 of 
the first edition: the change was made by the direction of 
Professor Boole's interleaved copy. It was judged conve- 
nient to number the Chapters in the present volume in con- 
tinuation of those in the original work. 

All additions of my own are enclosed within square 
brackets. The sheets have been read by the Kev. J. Sephton, 
Fellow of St John's College, as well as by myself, and the 
volume is much indebted to his care and accuracy. Obvious 
mistakes in the manuscripts were of course corrected; thus, 
for example, the table at the end of the volume was calcu- 
lated by Mr Sephton, because the table in the manuscript 
was rendered erroneous by the use of a wrong sign in a 
formula. 
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CHAPTER XIX. 



ADDITIONS TO CHAPTER II. 



1. [In Chapter ii. Art. 9, two methods are given for 
solving the differential equation 

{ax + Jy + c) c&» + {ax + i'y + c) dy = 0.] 

But there exists another transformation by which the equa- 
tion may be reduced to, (because it may be constructed from), 
an equation in which the variables are separated. 

Assume as this equation 

{Ay' + C) dx'-^ {Ax'+ C) dy' = (1) 

and let x' = x + w^y, y' =:x + mj/. 

It will be seen that in these equations united we have as 
many constants as in the original equation. Now on substi- 
tuting in the assumed equation the values of x' and y', and 
comparing with the equation given, we deduce a system of 
relations equivalent to the following, viz.: 

The quantities wi^, m^ are roots of the quadratic 

am^ " {b + a)m + h' = 0. 

The quantities A, A\ (7, C are determined by the system 
of equations 

A-^A' = a, C+C'=^c, 

Am^ -f A'm^ = a, Cm^ + C'm^ = c, 
B.D.E. II. ^^ 
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from which we find 



^ aWj — a 






m^ — ^Wg 



Now (1) gives on dividing by {A'x\+ C) {Ay + C) and 
integrating 

^ log (^ V + C') + 2 ^^g (^2^' + ^) = ^''''^*-' 

or (-4 V + (?') ^' (^y + GY = const, 

which on substitution and reduction gives 

1 

{ (am. — a) (x + m.y) + cm. — c'|«"»i-«' ^ . . 

-^-^^ — ^^ — ^^-^ ■ — J— = const (2) 



{ (aTWg — a) {x + m^y) + cm^ — c'}«"*a-«' 

2. Under certain circumstances the general solutions of 
differential equations of the first order fail. This happens in 
the above example if w^ = w^, the solution then reducing to 

1 = const. 

The theory of the deduction of the true limiting form of 
the solution in such cases requires a distinct statement. 

Let the supposed general solution be represented by 

G being the arbitrary constant and u a function of x, y, and 
constants which are not arbitrary. Suppose too that when 
one of these constants k assumes a particular value k, the 
function u reduces to a constant v. Then we have 

?/ — V G — V 

fC ~~' /C f€ ~~' tC 

Now the second member being a function of an arbitrary , 
constant is equivalent to an arbitrary constant and may be 
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replaced by C. The first member is a vanishing fraction, the 

limiting value of which is (;7t), the brackets being used to 

denote that after the differentiation k is to be made equal to k. 
Hence the solution becomes 



QH- 



In applying this theory to the reduction of the general 
solution (2) in the ease in which m^ = m^, it must be 
observed that the numerator of the first member is the same 
function of m , x, y, as the denominator is o{ m^, x, y ; or 
attending solely to their functional character with respect to 
m^^ Wj, we may affirm that the numerator is the same function 
of m^ as the denominator is of m , Representing these func- 
tions by ^[m^i ^(^J respectively, we have 

But wi^, Wj being roots of a quadratic equation may be 
represented in the form 

m^ = «i + ^, Wj = m — A;, 
the roots becoming equal when h = 0. Hence 

<f>(m — k)* 

Therefore since 

d<^ {m + Ic) d^ (m + k) d<f) (m — k) __ d^{m-' k) 

dk "" dm ' dk dm ' 

■ 

we have 

du '^^'^-^^dic '^("^•^^^ dk . 

dk~ 16 (w- A)}' 
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,, ^ /au\ ^^ ^ dm dm 

therefore (^j j-j^ ^^ 

c 

Thus the solution becomes on putting (7 for — , 

£log0H = (7, 

/7 1 
^^ 5m ^^r=^' ^^^ {{am-a) {x + my) +07^1- c] = (7. 

3. [The next Article seems to have been intended to ap- 
pear in the enlarged form of Chap, ii.; but 1 cannot discover 
what precise position it would have occupied. I conjecture 
that " the above demonstration" refers to Chap. il. Arts. 2, 3; 
and I have accordingly supplied a reference to equation (3) of 
Chap. II. 

I had myself drawn Professor Boole's attention to Chap. ii. 
Arts. 2, 3. The geometrical process of Chap. li. Art. 3, ap- 
pears to have been first given by D'Alembert in his Opus- 
cules, Vol. IV. p. 255. D'Alembert calls it a demonstration; it 
seems to me only an illustration^ at least in the brief form of 
the text : and that such was Cauchy's opinion may perhaps 
be inferred from the elaborate investigation given by Moigno, 
to which Professor Boole refers in Art. 5 of the present 
Chapter. 

I had also drawn Professor Boole's attention to the state- 
ment at the end of Chap. il. Art. 12, that only one arbitrary 
constant was involved. Accordingly Article 5 of the present 
Chapter developes this statement, and Article 4 seems intended 
to bear on the same subject.] 

4. In the above demonstration the relation between y and 
X is regarded as one of pure magnitude, and the interpreta- 
tion of the differential equation becomes a limiting case of 
that of the equation oi finite differences (Eq. (3), Chap. ii.). 
But if we represent x and y by the rectangular co-ordinates 
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of a moving point on a plane the differential jsquation may be 
interpreted directly. For supposing it reduced to the form 

we see that the direction of motion is constantly assigned as 
a function of the co-ordinates of position. The entire motion 
is therefore determinate as soon as the initial point is fixed. 
The result of the motion is a line or curve wholly continuous 
or subject to irregularities according to the nature of the func- 
tion /(a?, y). That the arbitrariness of origin is geometri- 
cally equivalent to the appearance of a single arbitrary con- 
stant in the relation connecting x and y may be shewn thus. 

Let y = <A(^o»yo>^) 

be the relation between x and y indicated by the supposed 
motion, x^, y^ being the initial point of departure. Then this 
point being on the line of motion, x^^ y^ are particular values 
of X and y, so that we have from the above equation 

which establishes a relation between x^ and y^^ and shews 
that there exists virtually but one arbitrary constant. 

5. It is proved in»Art. 3, Chap, ii., that the constants 
a?^, y^, initial values of the variables a?, y in the solution of 
the differential equation of the first order, are necessarily 
equivalent to one arbitrary constant. I shall shew from the 
form of the above solution that this a priori condition is 
actually satisfied. 

Developing the expression for y [see Eq. (30) of Chap, ii.] 
in ascending powers of x, we have 



Ajx? Ajx 



.8 



y=J. + J.aj + -^ + ^-^ + &c., (32) 

in which J, = ^ /. (^j y.) (- a^.)""^ 

1.2...(n — rj 

the summation extending from n = r to w = oo. Forming 
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hence the differential coefficients of A^ with respect to x^ and 
y^, and reducing by (28), we shall find 

whence in particular 

^•H-/.(«.,ri^^-0. 



Eliminate between these equations^ (x^, 2^^, and we have 

= 0. 



dA^ (IAq dA^ dA^ 



^0 ^yo ^yo ^^t 

Therefore, by Prop/ 1., A^ is a function of -4^, so that tlie 
solution reduced to the form (32) contains but the single 
arbitrary constant A^. 

It remains to notice that the solution must be applied 
only under the conditions of convergency, i. e. under the con- 
dition that the ratio of the ri^ to the (n — 1)*^ term tends to a 
limit less than unity as n tends to infinity. For a discus- 
sion of the failing cases of this test see ' Finite Differences,' 
Chap. V. Generally it is desirable, in order to secure rapid 
convergency, to divide the interval x-^x^ into separate equal 
portions, to each of which the general theorem of solution 
may be applied. \i x — x^ be very small the theorem may 
be approximately represented by 

y-yi>=/(i»a»yJ(«^-^o)- 

On these principles Cauchy has founded remarkable methods 
of solution, which deserve attention from the commentary on 
the limits of error on their application by which they are 
accompanied (Moigno, Vol. 11. pp. 385 — 434). 
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CHAPTER XX. 



ADDITIONS TO CHAPTER VII. 



1. [This Article relates to Art. 2 of Chap, vii.] 

The sense in which (9) may be said to constitute the 
general solution of the differential equation is this. We 
obtain from it 

giving any particular value to C this will geometrically 
represent a curve consisting of two branches, and giving to 
G every possible value we obtain an infinite system of such 
curves, each consisting of two branches. The aggregate of 
branches thus obtained is evidently the same as the aggre- 
gate of curves given by the two primitives (5) and (6), un- 
restricted by any connexion between c^ and c,. In this sense 
then the solution (9) is general, that it includes all the parti- 
cular relations between y and x which are deducible from 
the original primitives (5) and (6). And it is only in this 
sense not general that it groups these relations together in a 
particular manner. 

To the expression of the complete primitive a certain 
variety of form may be given without affecting its generality 
in the sense above affirmed. Thus, if to the solutions of the 
component differential equations we give the forms 

ye"**- Cj = 0, logy + oa? - Cj= 0, 

we should have, by the same procedure, as the expression of 
the complete primitive, 

(ye"*" - c) (logy + aa; - c) = 0, 
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an equation which may equally with (9) be regarded as tlie 
complete primitive of the differential equation given, anrl 
which in geometry represents the same totality of branches of 
curves as (9), with this difference only, that they are differ- 
ently paired together. 

2. [This Article relates to Art. 3 of Chap, vii.] 

The question will here naturally arise, Since if F= c be 

a solution of one of the component differential equations, 

f[V) = c,m which /(F) denotes any function of F, is also a 

solution, by Chap. iv. Art. 3, why not give to the complete 

primitive the form 

• lAiVd - c] IMV,) - c] {/.(F„)-cl=0, 

or the stricter form 

/.(F;)/.(F,) ./.(F„) = {F), 

in which /i(F,), ^(FJ, ... jC(l^n) denote arbitrary functions 
of Fj, Fg,..., V^ respectively — stricter because the presence of 
arbitrary constants and functions in the previous form is a 
superfluous generality ? It is replied that though the form 
just given is analytically more general than (15), it is not 
more general than (15) with such freedom as is permitted 
in the interpretation of the arbitrary constants. In a physi- 
cal or geometrical application we should not only be per- 
mitted to assign a particular value to the arbitrary constant 
in (15), so deducing what in reference to its source would 
then be termed a particular primitive, but to combine the re- 
sults of different determinations of c together, so as to obtain 
every form of solution which is implied either in the func- 
tional equation {F), or in its component primitives 

The same considerations justify us in speaking of (15) as 
the complete primitive, and not as a complete primitive. 
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CHAPTER XXI. 



ADDITIONS TO CHAPTER VIII. 



1. [The Singular Solutions of Differential Equations of 
the First Order received great attention from Professor Boole, 
and the Chapter devoted to that subject is one of the most 
valuable and important in his work. He continued his re- 
searches after the publication of his first edition, and intended 
to reconstruct the Chapter with great improvements in the 
second edition. After carefully examining the manuscripts I 
came to the conclusion that it would be very difficult to re- 
write this portion of the work so as to connect the old matter 
with the new; and thus it seemed best to reprint the original 
Chapter viil. with corrections of obvious misprints, and to 
print the matter intended for the revised form in the present 
volume. The plan gives rise to some repetition ; but this 
seems unimportant, compared with the advantage of preserv- 
ing in the author's own language all that he left on an in- 
teresting and important point which he had carefully studied. 

2. It may be of service to the student to reproduce the 
substance of some remarks on his Chapter viil. which were 
sent to Professor Boole soon after the publication of his first 
edition ; for there is evidence in his manuscripts that he paid 
great attention to such remarks while engaged in the revision 
of his work, and thus the reason and the meaning of some of 
his additions and changes may be made more obvious. These 
remarks will occupy the next Article. 

3. The two pages beginning with " And these conditions 
are sufficient," and ending with " do not lead to confiictiuq 
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results,*^ forming part of Arts. 3 and 4 of Chapter viii., seem 
obscure and difficult. The following may perhaps be substi- 
tuted with advantage. 

The only ways in which 

dy ^ df{x, c) ^^^ dy ^ df{x, c) ^ df[x, c) dc 
dx dx dx dx dc dx 

can be equivalent when c is variable, are 

(1) wben^ = 0, 

(2) when ^^^=00; 

in the latter case -^ = go, and therefore -,- = 0, and this 

ax ay 

implies that the singular solution is of the form x = constant. 

Thus there can be no singular solutions except such as 

df (x c) 
are found from ^^—K- — - = 0, and such as are found from 

dc 

X = constant. 

Similarly, if the complete primitive be expressed in the 
form x = F{y, c), there can be no singular solutions except 

such as are found from ^ — = 0, and such as are found 

dc 

from y = constant. 

In Art. 8 of Chapter viii. we read, " We may pass over 
the case in which the above equation is satisfied independ- 
ently of c, because the relation obtained would involve x 

dt) 
only, while it is a condition accompanying the use of -t- = cc 

that it leads to solutions involving y at least." It is ob- 
jected, Why may we pass over this case? Such a case might 
occur and furnish a solution, and then we should want to 
know the character of that solution. Take for example 

p = x^'y ; here if n is negative, -^ is infinite when a? = 0, and 



j,«+i 



this is a singular solution. For the general solution isy =C€'»+i, 
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and 80 a; = is not a case of it. The words — while it is 

a condition.., at least — seem very diflScult, for by supposition 

dt) 
we are now investigating what is furnished by ^ = oo , 

Professor Boole met the objection in substance thus : 

" It will be found that the rules in the book are correct in 
this case. What is implied in the Chapter, though not stated 

with sufficient clearness, is that if -^ = oo leads to a solution 

which does not involve y in its expression, nothing is to be 
inferred whether it is singular or not. Then the proper test is 

dx \pj 

In this example we have 

-^ = CO gives a?" = GO ; no inference ; 
(Iff 



dx \pj 



00 gives x'^'^^^^y'^ = oo . 



Hence a? =. 0, provided n is between and — 1, or y = 0. 

Consider these separately : 

First. Let w be between and— 1, and 05 = 0. This is 
bv the test a singular solution. Substituting it in the com- 
plete primitive we get y = c, which confirms this. 

Second. Let y = 0. This satisfies the dificrential equa- 
tion; but from the fact that it comes from -7-^ f-j = co we 
have no inference ; from the fact that it does not come from 
-^ = Qo we have the inference that it is a particular integral : it 
corresponds to c = 0. 
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There remains the case of a; = when n is between — 1 

J / 1 \ 
and — 00 . As this does not satisfy -7- ( - ) = oo , we infer that 

it is a particular integral. To prove this we have 



When a? = this gives, since 1 + w is negative, 

c = oo orc = — 00, 

according as y is positive or negative. This is like Ex. 2 of 
Chap. VIII. Art. 8." 

The remark made by Prpfessor Boole in the above reply, 

that «y -p = 00 Uads to a solution which does not involve y 

nothing is to he in f erred... \s> important. It corrects the stater 
ment put too strongly in Chap. viii. Art. 7, " All we can affirm 

is that if -^ = CO gives a solution at all it will be a singular 

solution." 

From Art. 8 onwards it seems assumed that a solution for 

which -^ = is always to count as a singular solution, even if 

it should coincide with a particular integral. This does not 
seem to have been quite the view of the former part of Chap- 
ter VIII. : see Arts. 5 and 6 of the Chapter. 

In Ex. 3 of Art. 9 we read, " the second is obviously a 
singular solution." This means that since we have a solu- 
tion which makes -^ infinite, we conclude that it is a singular 
solution. 

So in Ex. 5 of Art. 11 we read, " is evidently a singular 
solution," when it seems better to say, " and is therefore a 
singular solution." 

4. The additional matter relating to Chapter viil. begins 
with another example which was to be placed at the close of 
Art. 3 of that Chapter.] 



ART. 4.] ADDITIONS TO CHAPTER VIII. 13 

Ex. The differential equation 

(Vu;' + y' - w' - y) ^ - a? = 

has for its complete primitive 

Va?'* + y^ — Tff — y — c = 0. 
Here ^^ = ^ -\ ^^ ^ ^ 

^ = -1. 
do 

dy ^ Va?* + y*—m* dx _ Vaj* -j-y^ — m^ 

dc y — Vic* •\- y^— 7f?^ dc X 

Both ^ and -j- vanish then if 
dc dc 

This therefore is the singular solution and it satisfies botli 
the tests, as both x and y are contained in its expression. 

Of the partial tests 

d<^ _ rw d<i> __ d4> __ 

dc ^ dx ^ dy * 

the first is not satisfied, the last two are satisfied. 

The determination of c as a function of x by the solution 

df(x c) . . 

of the equation -^ ^ ^ = is equivalent to determining 

what particular primitive has contact with the envelope at 
that point of the latter which corresponds to a given value 
of X. 

One important remark yet remains. 'The elimination of c 
between a primitive y=^f{x,c) and the derived equation 

dy 

— = 0, does not necessarily lead to a singular solution in the 
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sense above explained. For it is possible that the derived 
equation 

df[x, c) 



dc 



= 



may neither on the one hand enable us to determine c as a 
function of a;, so leading to a singular solution ; nor, on the 
other hand, as an absolute constant, so leading to a particular 
primitive. Thus the particular primitive 

dv 
being given, the condition -^ gives 

whence c is + oo if a? be negative, and — go if a? be positive. 
It is a dependent constant. The resulting solution y = 
does not then represent an envelope of the curves of particu- 
lar primitives, nor strictly one of those curves. It represents 
a curve formed of branches from two of them. It is most 
fitly characterized as a particular primitive marked by a sin- 
gularity in the mode of its derivation from the complete pri- 
mitive. 

All the foregoing observations and conclusions may be 
extended to the case of solutions derived from the condition 

dc 

5. We have seen that the equation ;^ = may be satisfied 

by an absolutely constant value of c, so leading to a particu- 
lar primitive and not a singular solution. In this case 
'^(x'\-}i^ c) as well as '^^ {x, c) would vanish, and the nume- 
rator of (9), instead of being the difference of a finite and an 
infinite quantity, would be the difference of two infinite and 
equal quantities. [See Chap. viii. Art. 8.] It would not there- 
fore be infinite. Hence we conclude that -y- would not become 

ay 

infinite for a particular primitive in the strict sense of that 



ART. 5.] ADDITIONS TO CHAPTER VIII. 15 

term, i. e. for a solution derived from the complete primitive 
by giving to c an absolutely constant value. 

This is one point of contrast between the conditions 

^2^ = ^ = 00. 
dc ^ dy 

There is another not less important. As the numerator 
of (9) may become infinite not only when -^ {x, c) = 0, but 
also when -^ (a?, c) = infinite, we see that a relation between 

y and x which makes -t-- infinite will not necessarily satisfy 

the differential equation. On the other hand, it is not a par- 
ticular primitive in the strict sense of that term. 

. . dx 
Exactly in the same way the condition -y- = 0, as relating 

to the complete primitive, leads to the condition 



iG)=^' 



as relating to the differential equation, with the same points of 
difference in the respective applications. 

Ex. Let -^^my'^ , and suppose m a positive constant 

greater than 1. 

dv - 

Here fy^^"^" ^^ '^^ 

which becomes infinite when y = 0. As this involves y and 
satisfies the differential equation it is a singular solution. 

To confirm this conclusion we may refer to the complete 
primitive 

y=(a;-cr, 

which does not give y = for any particular value of c. 

Now let m be a positive constant less than 1. We have 

dn 

still -j- — ^ when y = ; but this value of y no longer satisr 
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fies the differential equation. It is not a solution at all, nor 

would it result from the application of the condition ^- = 

to the complete primitive. The distinction of character of 
the two tests is here made manifest. 

6. We may express the most important results of the 
foregoing investigations in the following theorem. 

Theorem. Every solution of a differential equation of the 
first order which is derived from the complete primitive by 
giving to c a variable value will, if it involve y m its expres- 
sion, satisfy the condition 

dp 
and if it involve a;, the relation 



-("-1=00. 

dx \v) 



But relations satisfying these conditions will not neces- 
sarily be solutions of the differential equation. 

In applying this theorem the following points must be 
carefully attended to. 

1st. No conclusion can be drawn from the satisfying of 
the condition -^ = oo when the relation in question does not 
contain y in its expression, nor from the satisfying of 

dx \p) 

when the relation in question does not involve x in its ex- 
pression. For these conditions being respectively derived 

d'U dx 

from ;7-=0 and ;7-=0 are subject to the same limitations 

in their application. 

2ndly. It may be that -f- or j- i-A assumes for a particu- 
lar relation between x and y the indefinite form - . In this 

^ 
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case we must seek by the development of its terms or by 
other known modes its true limiting value or values. Finite 
values will indicate particular primitives, infinite values sin- 
gular solutions, and when such values emerge together out of 
the same relation between the variables, the solution will be 
a particular primitive possessing the geometrical properties of 
a singular solution. Its locus will be a particular curve en- 
velopmg other curves of the same family. 

See Examples 2 and 3 of Chap. vill. Art. 11* 
We have seen that the conditions 

dy ' dx \p) 

indicate in general the existence of a relation between c and 
a; or c and y. And when that relation is such as to Enable us 
to detdHnine c as a continuous function of one of the vari- 
ables, the corresponding solution of the differential equation 
is singular, and is geometrically represented by an envelope 
of the curves of primitives. But it may be, as we have seen 
in a particular example, that the relation does not determine c 
as a function of a? or y ; but according to the language already 
used, c is a dependent constant, or in some other way different 
from the constant of an ordinary particular primitive. Let 
us examine in particular instances the kind of singularity 
which may hence arise. 

Ex.l. Given ^ = ?^^^. 

"^ X 

This becomes infinite if a; = 0; but this not involving y 
must be rejected. Again, it becomes infinite if y = 0, and 
this proves to be a solution of the differential equation, the 
limiting value of the indeterminate function in the second 
member being (Todhunter's Differential Calculus^ Chap. x.). 
Now the complete primitive is y = e", discussed in Art. 4. 
The constant c is there shewn to be dependent, the sola- 

B.D.E. II. ' "L 
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tion y = emerging from the complete primitive by making 
c =5 — Qo if a; be positive, and c = oo if a; be negative. 

Ex.2. Given (|y-..2^| + yMogy=0. 

Here o^y ±y(:i^-4log.y)4 



? = 



h 



therefore ±^ ^ ±(^-^^o^y)' :f __L 

dy 2 (a;'-4logy)i 

and this is made infinite by y = and by ai* — 4logy=0, 
i.e. by 

y = 0, y = €*. 

Both satisfy the differential equation. 

Now the complete primitive is 

We see at once therefore that the second of the above solu- 
tions is singular. The first however is deducible from the 
complete primitive by making c = oo or c = — oo, irrespec- 
tively of the sign or value of a;, provided only that x he 
finite ; not so however if x be infinite. The value of c is not 
therefore in the most absolute sense independent of that of x. 
If from the complete primitive we seek the singular solution 

by the condition -^ = 0, we get the two equations 

€«^=0, aj-2c = 0. 

The second of these determines c as a function of x, and 
leads to the second of the solutions obtained above. The first, 
though it does not determine c as a function of x, still ex- 
presses a relation between c and x, which is the ground of 
th« fulfilment of the condition 

dp 
dy 
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We may further notice a peculiarity arising from this rela- 
tion. Supposing X finite and the solution y = a particular 
integral, it presents the singularity that it is the only case in 
which two particular integrals agree. We might in any com- 
plete primitive, by changing c into c', get two values of c for 
the same particular integral, but then it would be for every 
particular integral. 

One negative character seems indeed to mark all the cases 
in which a solution involving y in its expression satisfies the 

condition -^ = oo . It is that such solutions do not emerge 
from the complete primitive by the attributing of a single and 
absolutely constant value to c. The relation which makes -^ 
infinite satisfies the differential equation only because it satis- 
fies the condition -^"=^9 ^.nd this implies a connexion be- 
tween c and a?, which is the ground of a real though it may 
be unimportant singularity in the solution itself. 

At this point, then, the question arises, whether the term 
singular solution shall be confined to that class of solutions, 
the loci of which represent the envelopes of curves of primi- 
tives, or shall be extended to all solutions which, satisfying the 

condition -f- = oo , indicate the existence of a relation be- 

tween c and x, and possess an actual singularity arising from 
this source. While the all but universal consent of mathe- 
maticians is in favour of the former course, it is to be remem- 
bered that the question is solely one of definition. Not such 
is the question how singular solutions of the envelope species, 
or as would more generally be said true singular solutions, 
are to be distinguished from all other solutions. This we 
now propose to consider. The question is not an isolated one. 
It stands in close relation to a series of properties of singular 
solutions which admit of an orderly development. 



^— 'J. 
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Discrimination of singular solutions of the envehpe species. 



7. A negative test, which in the great majority of cases 
suffices for the present object, is suggested by the following 
consideration. 

The differential equation determining -^ as a function of 

X and y determines also t4 > tb > <^ if^f^'^ and the know- 
ledge of these enables us to construct in a developed form 
the complete primitive. See Chap, ii. Art. 12. 

The values of ^ , -j4 > &c. ad inf.y as derived from the 

differential equation, are the same as those derived from the 
complete primitive. 

But a solution deduced from the condition -^ = oo is only 

constructed so as to yield the same value of ~ as the given 

differential equation does. If it be of the envelope species, 
tlie curve it represents has in general no continuous contact 
with the curve of any particular primitive. It will not there- 
fore generally yield the same values for ^, -7^, &c. as 

the differential equation does. It will not therefore generally 
satisfy the differential equations of an order higher than the 
first, which would be derived from the given equation by dif- 
ferentiation. Hence we have the following Proposition. 

Proposition. If a relation which makes •—- infinite satisfy 

the given differential equation of the first order , hut do not 
satisfy all the higher differential equations obtained from ity 
such solution will be singular and of the envelope species. 
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Ex. 1, By comparison with its complete primitive we saw 
in Art. 5 that ^ = wy ~ has for a singular solution y = 
when m is a constant greater than 1, 

We will first suppose m a fractional quantity greater than 
1, and endeavour to deduce the character of the solution with- 
out making use of the complete primitive. 

From the solution we have 

g=0, g=0,.&c. <«?;«/. 

But from the differential equation 

d^v -- dv -2^ 

J = (m-l)y«^ = m(m-l)y"', 

I 

and generally 



^ = TO(»»-l)...(m-r + %". 
Hence, when r is less than m, the substitution of y = gives 

das' " 
as before. But if r is greater than m, it giyes 

We conclude that the solution is of the envelope species. 

Secondly, suppose m a positive integer greater than 1. 

In this case we find, when r is less than ?», the same series 
of values as before ; but fox r =m we have 

^jl = w(m-l)...l, 
and this also shews the solution to be of the envelope species. 
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Ex. 2. The differential equation 

is satisfied hy 

Is this a singular solution or a particular integral ? 

From the solution we find 

<?y __ X cPy __ 4 
dx^ y* do? y** 

From the differential equation we shall have 

(^aj*"* 2(y-a?p) ' 

substituting in which the value of -^ , obtained from the 
proposed solution, we find 

A^ C^ + y')'^ 8 

cfo« 2^* y'* 

Now this differing from the value before obtained, we con- 
clude that the solution is singular and of the envelope species. 

And this result is verified by comparing the solution with 
the complete primitive 

(aj-c)»+(2/-Vr:7y=;. 

As the test above exemplified is merely negative, it is in- 
sufficient. For it is conceivable that an enveloping curve 
should have an infinite order of contact with each of the curves 
which it envelopes, and this is also possible. Any test found- 
ed upon a comparison of the values of differential coefficients, 
any test therefore furnished by the Differential Calculus, would 
be insufficient for the discrimination of such cases. 



Ex. 3. Given :# = y (log yf. 



dx 



ART. 8.} ADDITIONS TO CHAPTER VIII. 23 

Here -^ = 00 gives y = 0, and this . satisfies tHe differen-r 
tial equation. 

From this solution we find 

^ = 0, ^ = 0, &c. <w? inf. 
From the differential equation we have 

S=y{aogy)* + 2(kgj,)'}, 

which consists of y multiplied by a rational and entire func- 
tion of logy. It is easy to see that all the higher differential 
coefficients of y hence derived will possess jthe same character. 
And all such vanish with y. 

We can therefore neither affirm nor deny that the proposed 
solution is of the envelope species. 

8. Before demonstrating a general Rule for the discrimi- 
nation of solutions of this character, we shall notice certain of 
their properties which serve to indicate in what direction the 
Rule IS to be sought. [See Chap. viii. Art. 14.] 

As the exact differential equation differs from the sup- 
posed given differential equation by having acquired a factor 
which the singular solution makes infinite, so the given dif- 
ferential equation may be said to differ from the correspond- 
ing exact one by containing a factor which the singular solu- 
tion makes to vanish. If we knew that factor, we could by 
rejecting it reduce the given differential equation to a form in 
which it would no longer be satisfied by the singular solution. 
Now Poisson has shewn on a particular assumption, which 
does not however affect the principle of the demonstration, 
that this factor can be found when the singular solution is 
known. His demonstration is in substance as follows. 

Let us represent the given singular solution of the dif- 
ferential equation by 

w = 0, 

u being a given function of x and y. Then introducing u and 
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X instead of y and x as variables, the differential equation 
after transformation will assume the form 

Now this equation being satisfied by t« = and the first 
member vanishing, the second must also. Poisson now 
assumes, and the assumption must be carefully noted, this 
second member to be capable of being developed in ascending 
positive powers of w. Supposing it so developed, the diffe- 
rential equation becomes 

-T- = -4w* -f Bifi + &c. , 

in which A, B,... are functions of x, and a, A--- ascending 
positive indices. 

Hence if w = be a singular solution we have, putting p 

-^ == Jatt«-i + 5)8w^i + &c. = GO . 

Bat this demands that there should be at least one nega- 
tive power of u in the development in the second member. 
Therefore a— 1, the lowest index, must be negative. There- 
fore a being already positive must lie between and 1. 

We may give therefore to the transformed differential 
equation the form 

dx ^' 

a being a positive fraction, and Q not vanishing with u. 
Hence,^ dividing by w*. 
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a differential equation which is not satisfied by t« = 0, since 
w = gives w^"" = 0, and the first member vanishes while the 
second member does not vanish. In its present form then 
the equation is not satisfied by w = 0. We see also that the 
property of being satisfied by m = has been lost not in 
reality through a transformation, but through the rejection of 
an algebraic factor w* from the transformed equation. It has 
been shewn in the treatment of Clairaut's equation, how in 
the ascent by differentiation to an equation of a higher order 
a somewhat analogous effect is produced, the singular solu- 
tion emerging out of a factor of that higher equation. 

If we inquire what is essential in Poisson's demonstration, 
we shall find it to consist in that the transformed equation is 
of the form 

£-«^. 

in which while Q neither vanishes nor becomes infinite when 
u=0, the functions 

both vanish with u. The question whether U is of the form 
w* as Poisson supposes, or is not, is wholly immaterial. 
This will ftdly appear from the demonstration of the follow- 
ing theorem, which is in effect Poisson's fireed firom arbitrary 
assumptions. 

9. Proposition, ij^ w = he a solntion of a differential 
equation of the first order between y and a?, and 

represent the form which that equation assumes when u and x 

are assumed as variables instead of y and x, then iff{x, u) be 

resolved into two factors Q, U, of which Q neither vanishes 

nor becomes infinite when w = 0, while the functions U and 

f^du 

I -jj both vanish when m = 0, then the differential equation can 

be reduced to a form in which it shall cease to be satisfied by 
w=0. 
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In the statement of this proposition x is supposed to be 
constant in the integration relative to u. 

The differential equation after the transformation whicli 
introduces u and x as variables becomes 



Let 






so that V is in general a function of x and w, the form of 
which is known by integration when that of U is given. 
And again, transform the differential equation by making v 
and X the variables instead of u and x. We have 

/dv\ __ dv dv du 
\dx) dx du dx' 

dt) 
in which -j- is the differential coefficient of v with respect to 

Xy on the above hypothesis as to the constitution of t; as a 

function of x and w, while (-r-jis the differential coefficient 

on the hypothesis that v is reduced to a function of x alone 
by the conversion of u into a function of x. 

oi. dv 1 du ^y^ 

the above equation becomes 

'dv\ dv 



\dx) dx 

Now if w = give v = for all values of a?, it will there- 
fore give 

and further, 

^_ d p du^_r. 
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since we are permitted to make m = before effecting the 
differentiation with respect to x. Hence the equation re- 
duces to 

And this is not satisfied, since by hypothesis Q does not 
vanish with w, 

-^=0, the transformed differen- 
tial equation will no longer have t« = for a solution. 

Cor. Assuming Q = l, which does not violate the hypo- 
thesis respecting Q, and gives 

we see that if 

be satisfied by w = 0, and if at the same time w = gives 

du 



f 

J A 



/('". «) 



= 0, 



the differential equation can be transformed so as to cease to 
admit of the solution w = 0. 

It is obvious however that it is best to assume Q so as to 
make the subsequent integration for determining v the sim- 
plest possible. 

It is manifest that a solution which can thus be made to 
cease to satisfy the differential equation cannot be a particular 
primitive. For the complete primitive of the transformed 
differential equation which it does not satisfy is convertible 
into the complete primitive of the original differential equa- 
tion which it does satisfy, merely by writing therein for v its 
expression as a function of x and v. It cannot therefore be a 
case of the complete primitive m an^ sense. It must be 
a singular solution of the envelope species. 

The converse proposition still remains to be proved. 
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10. Proposition. If u=^Ohe a singular solution of Hie 
envelope species of a differential equation of the first order ^ and 
if hy assuming u and x as the variables^ the differential equa- 
tion is reduced to the form 



then will 



du 



£ 



become when u = 0. 

Let the complete primitive be represented by 

F{x, u) = (7, 
then, since 

dF(xy u) dF{xy u) ^ _ /x 
dx du dx ^ 

we have if for brevity we represent F{x, u) by F^ 



dF 
du _^ dx ^ 
^""dF' 

du 



* 



« 



dF 



•^0 dx 

Now w = being a singular solution, F[x^ 0) is not a con» 
stant ; for if it were, the complete primitive would, on giving 
to G the constant value in question, yield m = as a particu- 
lar primitive. And this would equally be the case whether 
that constant were finite or infinite in value. We see then 

that F{x, 0) must be a function of a?, and therefore — \ ' ^ 

must either be a function of a;, or a finite constant differing 
from ; the latter if F(x, 0) be of the form eifa? + J, the former 
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if it be not of that form. Therefore the value of ^ — - 

ax 

when w = 0, since in this we are permitted to make w = 

before differentiating with respect to aj, will be a function of 

Xj or a finite constant differing from 0. 

Now it is manifest that in general 

r- 1 dF ^ ^rdF J 



J dx 



_ ^ du 

dx 



where H is some value intermediate between the greatest and 
least values which -r^ assumes within the limits of integra- 
te 
tion. When these limits are, as in the above case, infinitesi- 
mal, we have 

dx 

TT , du 1 rdF , 

Hence I ir. r = iw, I -r- du 



J J[x,u)'' dF{x, 0) 1 du 

dx J^ 



i-^{F{x,u)-Fix,0)], 



dx 

But we have seen that , ' does not vanish. Hence 

dx 

its reciprocal, the first factor of the right-hand member of the 

above equation, does not become infinite. Again, 

jP(a?, tt)-J?'(a?, 0) 
vanishing when w = 0, we have 

'• du 



L 



of{x, U) ' 

when u is made infinitesimal as was to be shewn. 



30 ADDITIONS TO CHAPTER VIII. [CH. XXI. 

It Will be observed that the previous general expression 
for I -J-. : becomes infinite if w = is a particular integral. 

For then, F{Xf 0) being a constant, — ^ * ' vanishes, while 

F{x^ u) —Fix, 0) does not vanish so long as u differs by 
however small a quantity from 0. 

These propositions form the ground of the following Bule 
for the discrimination of singular solutions of the envelope 
species from all others. 

11. Rule. The j)roposed solution being represented by 
tt = 0, let the differential equation, transformed by making u 
and X the variables, be 

-^ +/(», «) = 0. 



dx 
Determine as a function of x and u the integral 



f 

J n 



in which Z7 is either equal to f{x, w), or to f{x, u) deprived 
of any factor which neither vanishes nor becomes infinite 
when w = 0. K that integral tend to with u the solution is 
singular. 

Ex. 1. Determine whether y = b is a singular solution or 
particular integral of the differential equation 

Here, since w = y, no preliminary transformation is needed. 
"V^e have I — ^ ^^^ 

which tends to with y. Hence the solution is singular. 
To verify this we observe that the complete primitive is 



f 
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and this cannot be reduced to y = bj giving any constant 
value to c. 

We have seen in Art. 7 that the test which is founded upon 
the comparison of differential coefficients does not suffice to 
characterize the above solution. 

Ex. 2. The equation ^ = ^-^^ is satisfied by y = 0. 
Is this solution singular or particular? 

Here also no transformation is required. We have, reject- 
ing the factor - which neither vanishes nor becomes infinite 
when y = 0, 

^^^ = log logy -log logo 

^^^logO* 

and this being infinite, however small y may be, may 

be said to tend to infinity as y tends to 0. The solution is 

therefore particular. 

It will perhaps appear at first sight as if in the above ex- 
ample we ought to write 

when y is made equal to 0. But the course of the demonstra- 
tion shews that the value of the definite integral must be first 
obtained on the hypothesis that u (in this case replaced by y) 
is finite, and then the limiting value which its expression 
approaches to, as u approaches to 0, be sought. And in this 
case, since for all finite values of u however small the integral 
is infinite, its limiting value is infinite. 

The complete primitive in the above case is 

and the nature of the solution y = has already been dis- 
cussed in Art. 4. 
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History of the Theory of Singular Solv/ttons. 

12. It is remarkable that while the theory of enveloping 
curves and surfaces was at once founded and developed by 
Leibnitz in 1692 — 4*, the corresponding theory of the singular 
solutions of differential equations has been of very slow growth. 
The existence of these solutions was first recognised in 1715 
by Brook Taylor ; it was scarcely more than recognised by 
Clairaut in 1734. Euler, in a special memoir, entitled Eocpo- 
sition de quehues Parcidoxes dans le Calcul Integral, published 
in the Memoirs of the Academy of Berlin for 1756, first made 
them a direct object of investigation ; but the foundations of 
their true theory were only laid in 1768 in his Institutiones 
Calculi Integr alts. Laplace, Lagrange, Legendre, Poisson, 
Cauchy, and De Morgan have in various ways developed and 
extended that theory; but there has been so remarkable a 
want of unity and connexion in this long series of researches, 
that important portions of the theory appearing in a too 
isolated form have been neglected, forgotten, and rediscovered. 
I purpose here to give a brief account of what seems most cha- 
racteristic, rather than of what is most original in their several 
researches ; for the germs of nearly all subsequent discoveries 
on the subject are to be found in the great work of Euler. 

Taylor and Clairaut appear to have been led by accident to 
the noticing of singular solutions ; the former while directly 
occupied on the solution of differential equations, the latter 
while discussing a remarkable class of problems relating to 
the connecting properties of different branches of the same 
curve. Taylor gave them the name singular, while Clairaut, 
and Euler too in his memoir, regarded them as a species of 
paradox, not merely from their non-inclusion in the general 
integral, but from the mode of their discovery through a 
process of differentiation. The memoir of Euler, though it 
sheds no light on the real nature of these solutions, contains 

• Ada Eruditorum, 1692, p. 168 ; 1694, p. 311. Opera, Tom. ra. pp. 264, 
296. 

Meihodua IncremerUorum, p. 26. 

M&moiru de I^Acad4mie des Sciences, 1734, p. 209. 
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an interesting theorem concerning their connexion with the 
form of the differential equation, viz. If this equation can 
be brought to the form 

Vdz^ZiPdx-^Qdy), 

in which 2 is a function of x and y, and Zoi z, then will 

be a singular solution. In his Instituttones Calculi Inteffralts, 
Tom. I. p. 393, however, Euler gives a rule which is the 
counterpart of that of Cauchy. [See Chap. viil. Art. 12.] 
He shews that if w = be a particular integral, and if the 
differential equation be reduced to the form 

then I -7-7 r = 00 . 

The limits of integration are here supplied. The reasoning, 
which is not fully developed, is the following. From the 
transformed equation we have 

^^- ^^ 
ax=i 



Hence x=G+j 



<f> {^, u) ' 
du 



x_ l^C 

o~ ■*" c] 



<\> (a?, u) ' 
du 



(j) {x, u) ' 



If this be satisfied by a solution involving x and y, and if 
that solution be a particular integral, then on putting for x 
its value in terms of u and integrating, the above equation 
will be satisfied by giving some particular constant value to 
C. But if the supposed particular integral be w = 0, then x 
and u being independent, we may perform the integration 
with respect to w as if a? were constant. The resulting equa- 
tion cannot be free from x unless C be infinite, and tlv^^ ^^. 



B. D.E. IT. 
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will evidently not be satisfied unless I--, r be infinite. 

J <f>(x, u) 

We infer then that this is a necessary condition in order that 
w = may be a particular integral. 

This is Euler's fundamental theorem, and from this, by 
means of an hypothesis agreeing with that of Poisson con- 
cerning the form of the transformed diflferential equation, he 
arrives at the condition 

dp _^ 

dy" 

[In the passage to which Professor Boole refers, Euler 
does not undertake to discuss the nature of any solution, 
but only of a solution of the form x = constant. On his 
page 408 Euler proceeds to discuss the nature of any solu- 
tion. Professor Boole seems to me to attribute too much 
to Euler. For the convenience of those who wish to ex- 
amine the original, I will give the reference to the passages 
in the later editions of Euler's Institutioiies Calculi Integralis : 
Vol. I. pages 343 and 355 of the edition of 1792 ; Vol. I. 
pages 342 and 354 of the edition of 1824.] 

Laplace in the Memoirs of the French Academy for 1772, 
p. 343, established the tests 

dp_ —(}l\- 

dy^ ' dx\p)~~ * 

and shewed their respective uses. He established also the 
test which consists in the comparison of diflferential coefficients, 
and he supposes it universal. He adopts the hypotliesis of 
his predecessors as to the forms of expansion, but with some 
recognition of its insuflSciency. 

Lagrange in the Memoirs of the Academy of Berlin for 
1774, p. 197, and 1779, p. 121, appears first to have developed 
the theory of singular solutions in its two forms of derivation 
from the complete primitive and derivation from the diflferen- 
tial equation, and to have established the essential connexion 
of tliese. But supposing the diflferential equation to be ex- 
pressible in the rational form 

F{x,y,p)^0, 
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and employing the differential coeflScients of F{x^y^p) in- 
stead of those of p he was led to sacrifice rigour to symme- 
try. One of his results has often since been adopted as a 
test of singular solutions. It may be thus stated. 

Prop. A singular solution makes the general value of 
-~^ , deduced from the differential equation in its rational and 

integral expression, to assume the form - . 

[The demonstration is given in Chap. VIII. Art. 14.] 

This ambiguity of value of -y^ is evidently but an expres- 
sion of the fact that the contact of a curve of the complete 
primitive and that of the singular solution is not in general 
of the second order. 

The result given in equation (5) of Chap. vill. Art. 14 has 
also been adopted as the test of singular ^lutions. 

The researches of Poisson and Cauchy have already been 
noticed. It is certainly remarkable that the final test to 
which Cauchy's analysis led should be essentially the same as 
that which had been discovered by Euler so long before. 

Professor De Morgan has thrown an important light upon 
the nature of the conditions 

dp _ ^P __ 

dy'^ ' dx"^ ' 

which are fiilfiUed by all singular solutions in the expression 

of which X and y are both involved. He has shewn that any 

relation between x and y which satisfies these conditions will 

d\ 
satisfy the differential equation unless it make -t4, as derived 

from the differential equation, infinite ; that it may satisfy the 

d %J 
differential equation even if it make -t-^ infinite ; lastly, that 
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if It do not satisfy the differential equation, the curve it 
represents is a locus of points of infinite curvature, usually 
cusps, in the curves of complete primitives. 

The proof is as follows : 
Let p = ^ {x, y) 

be the differential equation. Then the proposed conditions 
are 

dy ^ dx ' 

therefore by differentiation, 

dxdy dy^ dx ' dx'* dxdy dx ' 

whence we have 

d'4> d^ 

dy _ d^cdy _ dx"" 

dx~^ d^<^ "" d^<f> ' 

dy^ dx dy 

These are two equivalent expressions for the same value of 
-^•- . The question now is, under what circumstances this 

value of -p will satisfy the differential equation. 

Now from that equation we have by differentiation 

d^y _ d^ d<f) dy 
dx^ dx dy dx^ 

whence 

d^y d(f> 

dy _ dx^ dx 

dx d<f) ' ' 

dy 
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If then -r? be finite we have, since -^ and -^ are both 
(Jbr ax ay 

infinite, 

d4 

dy _^ dx 

dx dt^ 

dy 

and this hj the rule for the evaluation of fractions of the 
form ^ is equivalent to the value in either of its forms before 

obtained for -r^. Hence, any relation which satisfies the 

given conditions and makes -vj finite, will satisfy the diffe- 
rential equation. 

And the same result holds even if -^4 be infinite, provided 

that -t4 -r -7^ vanish. 
dor ay 

Lastly, as when this result does not hold, the failure is due 

d^V 
to the infinite value of -tt* "^® see that the line in which the 

locus of the proposed relation intersects the curves of primi- 
tives will be a locus of their points of infinite curvature. 

[Transactions of the Cambridge Philosophical Society, 
Vol. IX. Part II.] 

Legendre's Memoir of 1790 throws but little light upon 
the subject of this Chapter. But it exhibits the theory of 
the singular solutions of differential equations of the higher 
orders, both ordinaiy and partial, in a form of great beauty, 
and will be noticed in the proper places. 
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CHAPTER XXII. 



ADDITIONS TO CHAPTER IX. 



1. By successive application of the second theorem of 
Chap. IX. Art. 13, a ^linear equation of the nf^ order may be 
reduced to one of the (n — r)^ order, if r distinct integrals of 
what the given equation deprived of its second term would 
be are known. 

The reduction may however be effected immediately by 
the method of the variation of parameters. In this and in 
most general investigations connected with differential equa- 
tions great advantages in point of brevity and of the power of 
expression are gained by the employment of the symbol of 
summation S, and of the language of determinants. I shall 
exemplify this here. 

Suppose the given equation to be 

^■*"^^^'*"^»^^ + ^. = ^ (1), 

and let y^, y^^...yr he r particular values of y, satisfying the 
equation 

Thus y = c,y, + c^, ... + c^r 

is a solution of the latter equation including these particular • 
solutions. We shall represent this by 

y=^tc(yi (3), 

and regarding the quantities c^, c„... c^, represented here by 
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Ci as variable parameters, shall seek to determine them so that 
the above value of y may satisfy the equation given. 

These r parameters, enabling us to satisfy r — 1 arbitrary 
conditions, besides satisfying the differential equation, we may 
choose these so that 

dt/ d^y dr^y 
3i' d^'''"d^' 

may be the same inform as if c^, Cj, ... c^ were constant. Now. 
from (3) 

whence 

dx ^ dx^ 
provided that the condition 

be satisfied. Differentiating the first of these equations, we 
find in the same way that 

provided that the condition 

dx dx 

be satisfied. And thus continuing we see that the system of 
r equations 

y^tc^. ^=2c,^, ^ = Sc.^, (4), 

will hold true provided that the r — 1 conditions 
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be satisfied. In each of these equations the symbol S is to 
be interpreted by giving to « the successive values 1, 2,... r, 
and taking the sum of the results. 

Differentiating the last of the equations (4), we have 

As we cannot impose the condition that the last term of 
this equation shall vanish, let z represent its unknown value, 
then 

%-^^' («)• . 

Now the system of equations (5), together with 

y^^ dci_ 

constitute a system of r simple algebraic equations deter- 
mining by solution the r quantities 

dc^ dc^ dc^ 
dx' dx '**' dx 

in terms of their coefficients and of z, and therefore in terms 
of X and z, since the coefficients are known as functions of 
X. It is evident also that as the second members of all the 
equations but one vanish, and the second member of that is 
Zf the values so determined will be of the form 

XjjXjj,... X^ being known functions of x. Thus th^e r un- 
known quantities 3-^>"«;t^ ^ire made to depend upon only 

one unknown quantity, viz, z. It remains then to deter- 
mine z. 

For this purpose we must complete the expression of the 
differential coefficients of y, and substitute in the given dif- 
ferential equation, and then seek to satisfy that equation. 
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Now difierentiating (6) we have 

<^'y _<, ^Vi ,^dy,dctdz 



-s.^...(^^)..| 



on substituting for -j^ the value XiZ as above determined. 

We observe that the coefficient of z is here a known function 
of X. If we differentiate this equation and in the result sub- 

stitute as above for -j^ , we shall have a result of the form 
dC^y ^ dr^Vi , ^ , ^,dz . d^z 

L and M being known functions of x. Ultimately then 
we have 



g = 2c.S + P.+ gJ + ? 



Z 



.n-^ 



etc* dx"" dx'* ' dod 

Thus, while y and the differential coefficients of y up to the 
(r — 1)*** are of the same form as if c,, c,,... c^ were constant, 
the succeeding ones differ in containing an additional portion 
consisting of z, and differential coefficients of z multiplied by 
known functions of x. The result of substitution of these 
values in the given differential equation will therefore consist 
also of two classes of terms, viz. terms under the sign of 
summation, which will be the same in form as if c^, c^j.-.c^ 
were constant, and terms involving the differential coefficients 
of z up to the {n — r)***, with multipliers which are known 
functions of x. We shall in fact have 

V, ^4.^ ^^ + A ^^ +A ) 
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Now yi being by hypothesis an integral of (2), the first 
line of the above equation vanishes, and there remains the 
linear equation of the (n — r)^ order 

^^- + ^^^^ + ^=^- 

Supposing z hence determined, we have In general 

c<= jXizdx, 
and hence 

j/ = yjx^z dx + yjx^z dx + yrj^r^ dx, 

and as z will have w — r distinct values, each involving an ar- 
bitrary constant, the above equation will furnish n — r distinct 
values of y, each involving an arbitrary constant. It is to be 
observed that no arbitrary constant need be added in the inte- 
gration of the terms XiZ dx, for the effect of such additioip 
would only be to reproduce the known integrals Ciy<. In 
this way, however, the equation would represent the general 
integral of the differential equation given. 

2. Let us examine the form of the result in the particular 
case in which r = w — 1 . 



Here we have 



dx"^"^' dx^ 



from 7W = 0to7W = w— 2, then 
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Accordingly the differential equation for z will be 

|^S(ip),.A.=X. (,). 



Now the equations for determining 

dx ' dx '"* dx 

become on putting XiZ for -r-^ , and writing for brevity i/\ for 
^y* v''for^' &c 



dx 



dx' 









Whence, by the theory of determinants, 



Y^l_ dM 






Jf standing for the determinant 



Vxy y^r 



% 



«-i 
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Now the determinant is ultimately a function of x ; and 
such indeed that 



^/ d-yA_ IdM 



For 



dM_ ^ dM dyi ^ dM dyl ^ dM dyj 
dx dyi dx dyl dx '" dy^"^^ dx 



(«-») 



dM , 



dM 



dM 



=^d^.y^^^Wi'""'''^W^' 



.(»-») 



(8). 



Now M being homogeneous and of the first degree with 
respect to the quantities y^, y„ y,^, we have 






dM 



Hence 2 ^^ yl is what M becomes when in its expression 

Vi^ yn •••y*-i ^® changed into y/, y^, ... y'„^^, therefore it is 
what M becomes when two of its rows of elements become 
identical ; therefore it vanishes. In like manner all the other 
sums in (8) vanish excepting the last, {ory^^'^%.,..y^^*^^^ is 
not a row of elements of the determinant if. Thus we have 



dM 



dM 



dyi 



C*-^)^* dx 



Hence 



tx,yn=tii 



1 dM 



.("->) = 



1 dM 



M dyt^^ y^ M dx' 



Thus the equation (7) becomes 



dz 



therefore 
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Hence, since 

[1 dM , 
whence ^^^jUdy^^^ 

^ [1 dM , 
we have y = ^Vij jf dy^^^ ' 

z being given above. 

In the case of X= 0, we have 

whence 
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CHAPTER XXIII. 



ADDITIONS TO CHAPTER X. 



1. The theory of singular solutions of differential equations 
of the higher orders has been presented in the most complete 
form which it has yet received by Legendre. {MSmoires de 
r AcadSmie Boyale des Sciences^ 1790, p. 218.) He determines 
first the possible forms of these solutions considered as emerg- 
ing from the complete primitive by the variations of its arbi- 
trary constants, and secondly the theory of their derivation 
from the differential equation itself. I shall follow the same 
order, and shall in the end endeavour to point out in what 
respect Legendre's theory may be regarded as complete, and 
in what respect it is imperfect. 

Suppose the differential equation to be of the n*^ order, 
and let it when solved with respect to the highest differential 
coefficient of y be represented by 

3^n=<^(a?,3^,2^i,y,, ...y«_J (1), 

in which, for brevity, 

_dt/ __ d^y __ d^ 

Let also its complete primitive, solved with respect to y, 
be represented by 

y=f{x, a„ a„ ... aj (2), 

a^j a^, ...a^ being the arbitrary constants of the solution. If 
we differentiate (2) with respect to x, regarding a^, a,, ... a» 
no longer as constants but as functions of a?, so to be deter- 
mined as to leave the expressions for y^, ^g, ... y^ as functions 
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of ttj, ttg, ... a„ the same as before, we shall have, on, repre- 
senting the second member of (2) by/, 



-^4.^^14.^ ^a . ^f ^^ 



+ 



+ 



whence 



^^ dx^ da^ dx ^ da^ dx ^ da^ dx ' 



provided that 



df^ da^ ^ 



y^" dx' 



df da. 



da, dx da^^ dx 



^d£ da^ ^ ^^ 
da^ dx 



Differentiating on the same hypothesis the first of these two 
equations, we find in the same way 



y^ da?' 



provided that 



ePf da^ <fy da^ d^f da„^^ 



dxdoj^ dx dxda^ dx dxda^ dx 

And continuing thus, it results that the system 



_dl _^ 



V = — — 

*'" dx* 



(3), 



will be satisfied, i.e., that ^1, ^2* • • • Vn will have the same ex- 
pressions when a^, Og, . . . a^ are variable as they have when 
these are constant, provided that the law of their variation be 
determined by the conditions 



^da, df^da^^ df da^^^^ 



da^ dx da^ dx 



da^ dx 



"n 



^da,^_^l_da, , cf/ da,^^ 



da^dx dx da^dx dx 



da^dx dx 



....(4). 



rfy da^ dy da^ d^ da^ ^ 

da^do^'^ dx da^dx"^'^ dx da^dx""'^ dx 
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In this system the coefficients of 

dx ^ da? ' '* ' dx 

are known functions of a;, a^, a,, . . . a„ when the form of y is 
known. 

Eliminating 

da^ da^ da. 



dx^ dx ^ '" dx * 

we have a relation between x, a^, a^y , . . a^; and this relation, 
with the given complete primitive and the first w — 1 of the 
derived and reduced equations, viz., with 

_/ -^ -^ -Q" 

will enable us to eliminate a^, Og, . . . a„, and to obtain a rela- 
tion of the form 

^(^'^''l' S'--P)=^ (^)' 

This is a diflFerential equation of the (n — 1)*** order. It dif- 
fers in its origin from the given differential equation, in that 
a new relation between a:, aj, a^, . . . a„ has been employed in 
place of the n^ equation, derived by differentiation from the 
complete primitive, for the elimination of the constants. 

The differential equation of the (n — 1)*** order thus obtained 
has an integral expressing y in terms of a*, and w — 1 arbi- 
trary constants. This is the most general form of a singular 
solution of the differential equation. 

It is possible that the elimination of a^ja^, ...a„ may 
lead to a resulting differential equation which, instead of 
being of the order w — 1, is of the order w — 2, w — 3, &c. 
The complete Integral of such equation would be a singular 
solution of the differential equation. These possible types of 
solutions are distinguished by Legendre according to the 
number of arbitrary constants which th^y contain. A solu- 
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tion containing n — 1 arbitrary constants is called by him a 
singular solution of the first order ; one containing n — 2 ar- 
bitrary constants a singular solution of the second order, and 
so on. 

Adopting this language we might term the complete primi- 
tive a singular solution of the order 0. 

Lastly, any relation between x and y, which satisfies the 
given differential equation, will constitute a particular case, 
either of the complete primitive or of one of the general 
forms of singular solutions above defined. In the case of 
differential equations of the first order it is seen that no arbi- 
trary constant can appear in the expression of the singular 
solvition* 

Ex* The equation 

lias for its complete primitive 

y==-— +5aj+a* + ff' (6), 

'leqtdred its singnlar solution! 

Proceeding su3 above, we find on the hypothesis of a and h 
being variable parameters, the same formal expressions for 

T^ , ;t^ as if those parameters were constant, viz. 

ax 






J 



(7), 



provided that the variation of a and & be such as to satisfy 
the conditions 



( 



dx dx ^ 



(8). 



B. D.E. II. 
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EUminatlng hence -p /tnd -7- , we have 

a? 
2a — — — 2ja? = 0. 

And from this, the complete primitive, and the first of the 
derived equations (7) elimmating a and h^ we find 

This is the difierential equation of the first order, by the 
solution of which the most general form of the singular solu- 
tions of the given differential equation will be determined. 

Reducing it to the form 

(I6y + 4a?+a;*)* ■ 

and integrating, we find 
, (16y + 4a:' + aj*)*=a?(l + aj^* + log{ajf V(l + aO}+ 01 

This then is the general expression for the singular solu- 
tions of the given differential equation. We see that it in- 
volved in its expression one arbitrary constant. 

The differential equation (9) may properly be termed a sin- 
gular first integral of the given differential equation. The 
singular first integral (9) has itself also a singular solution, 
viz. 

^ 4 16 ' 

but this is not a solution of the original differential equation. 
Nor have we any right to expect that it should be so. A 
singular solution of a differential equation of the first order 
does not necessarily satisfy the differential equations of higher 
orders derived from that equation. Chapter xxii. Art. 7. 

2. It remains to establish the theory of the derivation of 
the singular solution from the differential equation without 
the mediation of the complete primitive. 
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Kesuming the differential equation in its reduced form (1), 
and representing its second member by ^, suppose an infini- 
tesimal variation given to the arbitrary constants of its com- 
plete primitive, and let the symbol 8 be used to denote the 
corresponding derived variations of y, y^, ...y^. Then we 
have 

and BO on. Hence, substituting and transposing, 

^_^^_^ ^V_&,.^0 (10). 

Let us consider the real nature of this equation. 

K a value of y, suppose y = '^(a?), satisfy the given differ- 
ential equation, that value substituted in the coefficients 

# _#_ &c 



'n— 2 



of the above equation will convert them into functions of a;, 
and the equation itself will become a linear differential equa- 
tion, the solution of which will determine Sy as a function of 
X. K the differential equation (10) be really, as it is appa- 
rently, of the ri"^ degree, Sy will have n arbitrary constants, 
Oj, • . . a», and will be of the form 

P^jP^^.m.P^ being functions of x. Hence 

y + Sy = ^(aj) + a,P,... + a„P,. 

We see thus that the given solution y = '^(a?) will be a 
particular case of this general integral involving n constants. 
It will therefore be a particular integral of the ^to^c^'Si^i^ 

4—^ 
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If, owing to the constitution of its coefficients, the difiEerential 
equation (10) be of the degree w — 1, we shall have 



y + gy = -^(aj) +ajPj .. . + o^jP, 



«-i> 



and y = '^{x) will then be a particular case of a solution 
involving w — 1 arbitrary constants. It w^ill therefore be a 
singular solution of the first order. Even so, if the diflferen- 
tial equation (10) be of the degree w — 2, y = "^(a?) will be a 
singular solution of the second order. And generally, if 
the differential equation be of the r^ degree, y = '^{x) will be 
a singular solution of the order n — r. 

Eesuming the equation (10) it is evident that it cannot 
be of the degree w — 1, unless -r-^ be infinite. For, dividing 

by -7-^, we have 

d<l>'dx'' dor-' ^^--"^ 



rfy 



n-1 



in which the first term does not vanish unless t^ be infi- 

nite. This then is the necessary condition for a singular 
solution of the first order. For one of the second order we 
must have in like manner 



dd> d(b 

" ^-— = QO , ^^— 1=00* 



and so on. 



It follows hence that to find the singular solutions of a 
differential equation of the n^ order, we ought to differentiate 

the equation, regarding y, -# , ^ , &c. as varying through 

the variation of the arbitrary constants, to form in this way 
a linear differential equation for 8y, to examine the conditions 
under which this equation reduces to the (n — 1)*^, or to a lower 
degree, and to examine whether the most general relation be- 
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tween x and y which satisfies such condition, satisfies also the 
given difierential equation. If so it may be regarded as a 
singular solution* 

Besuming the last Example, viz, 

and operating with B we have 

(2 dor \dx da?j) dbr 

which reduces to a linear dififerential equation of the first 
order for determining 8y, provided that we have 

Eliminating -y^ from the given equation hy means of this 
there results 



■.(^»)|-(.+«l.-fe=«. 



and we find on difierentiating this that it does constitute a 
solution of the given equation. It is therefore a singular 
first integral of that equation. We see that it agrees with 
the result obtained under the same name in the previous 
Article, and the rest of the solution need not be repeated. 

3. Upon Legendre's theory, and upon its results, the fol- 
lowing ODservations may be made. 

1st. We learn from it that there may exist ^ different 
general forms of the solution of a differential equation of the 
n^ order, viz. the complete primitive involving u a^iVitecr^ 
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constants, and general forms of singular solutions containing 
fewer than n arbitrary constants^ A solution y = '^{x) of 
unknown origin being given, we construct a differential equa* 
tion for determining cty, and, solving it, form the expression 
for y+%, and from the number of infinitesimal arbitrary 
constants it contains, determine the nature of that general 
value of y of which the given value is a particuliur case. 
Now we are not to infer from this that the form of y + Sy will 
be the same as the general value of y in question. But we 
may infer that it will be a form to which that general value 
is reducible. And the actual reduction will be effected by 
expressing the general solution (as is always possible) in a 
form permitting its expansion in ascending powers of the 
arbitrary constants, and in the expansion making these con- 
stants infinitesimal, and rejecting all powers of them above 
the first. In fact, if 

be any general form of solution which, when we assign to 
a^, a^y.,.a^ particular values (e.g. make them vanish) re- 
duces to 

y = '^(a?), 
then we shall have • 

the brackets denoting that after differentiation we make 
a^, a,, ... a^ vanish. 

This is that limiting form of the solution which Legendre's 
method enables us to construct by the solution of a linear dif- 
ferential equation ; and the ground of the sufficiencv of his 
method consists in this, that the infinitesimal quantities 

which are in fact the arbitrary constants of that solution, are 
equal in number to the arbitrary constants of the general 
unlimited solution, the nature of which is thus made known. 
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2ndly. Legendre's tests for differential equations of the 
higher orders are in kind and effect analogous to the tests 

dp ^ cZ 1 _ 

dy ' dx p 

for differential equations of the first order. They enable us 
to decide whether a solution possesses singularity, not whether 
it possesses the envelope species of singularity. The comple- 
tion of Legendre's theory would consist in the discovery of 
those further tests dependent upon integration which corre- 
spond to the test of Euler and Cauchy for differential equa- 
tions of the first order. 
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[CH. XXIT, 



CHAPTER XXIV. 



ADDITIONS TO OHAPTEE XIV. 



• [Art. 1 was intended to follow Chap. xiV. Art 2.] 

1. As the condition of dependence of functions of two 
variables is of fdndamental importance in connexion with the 
theory of ordinary differential equations, so the generalized 
condition of dependence of functions of any number of vari- 
ables forms a fundamental part of the theory of partial differ- 
ential equations. This is contained in the following proposi- 
tion. 

Prop, I. If w,, u^,...Un are functions of w^j a7j,...a?«, 
but are as such so related that some one of them is expressi- 
ble as a function of the others, or more generally that there 
exists among them some identical equation of the form 

F(u^, w„ ... u^) = 0, (1), 

so that as functions of aj^, a?,, ...ar^ they are not mutually 
independent, then, adopting the notation of determinants, the 
condition 



du. du^ 

du^ 
dx' 



du^ 
dx^' 



du^ 
dx. 



du^ 
dx. 



dUn 



du^ 
dx. 



= 



(2), 



Is identically satisfied. Conversely, if the above condition be 
identically satisfied, the functions Wj, w,, ... w^ are not mutu- 
ally independent in the sense above explained. 
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First let it he noticed that the Proposition is hut a general- 
ization of that of Chap. Ii. Supposing U and m to be two 
fiinctions of a; and y, the condition of their dependence is 
affirmed to be 

\ dU dU] 



-0, 



(h>^ dtf 

du du 
dx* dy 

i. e. it is the result of eliminating dx^ dy^ from the equations 



du J du ^ 



and therefore it is 



dU du _ dU du _ 
dx dy dy dx" ' 

as expressed in Chap. ii. 

We proceed to the general demonstration. 

Let the first member of (1), considered as a function of 
Wj, Wj, ... M^ be represented for brevity by F; then differen- 
tiating, we have 

dF . dF. ^ dF ^ . 

from which it follows that if rfu^, c?w,, ... du^^ are equal to 0, 
then is du^ equal to ; or, since w^, w,, ... u^ are functions of 
a?i, a?„ ...a?„, that if 



du. J , du. J . du. , 

du^ y . du^ , du^ y 



^'-^d^^ + ^dx,...+^dx. = 



dx^ 



dx. 



dx. 



J 



(3), 
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then is 

Thus the last n equations, linear with respect to 

are not independent, and therefore by the theory of linear 
equations the determinant of the system vanishes identically. 
Now this is expressed by the condition (2), 

It remains to prove the converse, viz. that if the condition 
(2) be identically satisfied, the functions w^, w,, ... w^ will not 
be mutually independent. 

First, the n — 1 fdnctions w., w„ ... w,^ are either mutually 
independent or not mutually independent. 

If not, then the n fimctions w^, w,, ... u^ are not mutually 
independent, and the Proposition to oe proved is granted. 

If Wj, Mj, ... w^j are mutually independent as functions of 
a?j, x^y ... a?^, they may be made to take the place of w — 1 of 
these quantities, e. g. ar^, a?,, ... x^^ in the expressing of w^, 
i.e. we ma^, by means of the expressions for Wj, w,, ... u^^^ 
eliminate from that of u^ the quantities a?i, a?,, ... a?^j, and so 
express t«^ as a function of u^y w,, ...w„_j and a?„. Suppose 
this done, then the system (3), (4) will be converted into 

Jwj = 0, {?t«j| = 0, {?M^^j = 0, 

Now, the determinant (2) vanishing, the equations of the 
linear system (3), (4) are not independent; therefore those of 
the transformed system, as written above, are not independ- 
ent; therefore the last equation of that system must be a 
consequence of the others which manifestly are independent. 
But from the form of that last equation we see that such can- 
not be the case unless we have 
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which implies that w„ is a function of «;^, w,, ... w^-i merely. 
Hence the functions u^^ u^^ ...u^ are not independent, as was 
to be shewn. 

The first member of the equation of condition (2) is com- 
monly called the functional determinant of w^, w,, ... w, with 
respect to x^^x^j...x^^ The proposition may therefore be 
expressed as follows. 

The condition of dependence or independence of any sys- 
tem of functions of as many yariables is the vanishing or 
non-vanishing of the fonctional determinant of the system. 

On account of the great importance of this proposition it is 
desirable to illustrate it by an example. 

Ex. Are the functions 

« + 2y + «, a? — 2y + 3;3?, 2ary - a?i5 + 4y« - 2a" 
mutually independent or not? 

The equation of condition is 

1, 2, 1 

1, -2, 3 =0, 

2y — Zf 2aj + 4i5, — a? + 4y — 4« 

that is, 

- 4 (- a? + 4y - 4^) + 8 (2y - «) - 2 (2x + 4«) = 0, 

which is identically satisfied. Hence the functions are de- 
pendent. In fact, representing them by Uy v^ w, we have 

4:WJ = ti* — t;*. 

[Art. 2 was intended to follow Chap. xiv. Art. 4.] 

2. As it has been shewn that a primitive 

u = <l>(v) (1) 

leads to a linear partial differential equation of the form 

Fp +. Qq^B .^V 
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provided that m = a, v = J, are integrals of the system of ordi- 
nary differential equations 

dx^dy dz ' f^. 

T'~Q^M ^^^' 

it is evident that we shall obtain a solution of the partial dif- 
ferential equation (2) by constructing the system of ordinary 
differential equations (3), deducing their general integrals 

and then constructing from these the primitive (1). 

But the question arises, Will this be the most general solu- 
tion of the partial differential equation given? 

That it will be so, may be shewn by means of Ihe general 
proposition. See Art. 1. 

For let w = represent any solution whatever of the given 
partial differential equation. Differentiating this with respect 
to X and y, we have 

dw dw ^ dw dw _^ 

dx dz-^^ '* dy dz^" ' 

substituting the values of ^ and q formed from this in the 
given equation, we have 

dx ay dz 
which must be identically satisfied* 

In like manner, u = a^ i? = & being solutions of the same 
equation, we find 

which must be identically satisfied. 
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Eliminating P, Q, jB from these three equations, it results 
that the Junctional determinant qfwj w, v, with respect to x, y, z, 
will identically vanish. Hence w ia a function of u and v, 
and the equation w = is A' particular case of ' 

which is thus shewn to be the general integral of the given 
equation. 

We are thus led to the following general Rule. 

EuLE. Ih integrate iihe emmtion Pp-^- Qq^R we must 
form the system of ordmary differential equations 

dx ^dy ^dz 

deduce their general integrals in the form 

t« = a, i? = &, 

and construct the equation 

F{u, v) = 0. 

This will be the general -solution sought. 

[Art. 3 was intended to follow Chap. XIV. Art. 5.] 

3. The above theory may be extended to linear partial dif- 
ferential equations of the first order, without regard to the 
number of the vaa^iables. • • 

First, the theory of the genesis of such equations is ex- 
pressed in the following proposition. 

Prop. A primitive equation of the form 

F{u,, w„ ... w«) = (1), 

in which w^, t«2,...w„ are any given functions of the vari- 
ables z, dependent, and x., a?,, . . . x^ indejjendent, will satisfy 
the linear partial differential equation obtained by eliminating 
dz, dx^y cfa?a, ....da?^ from 

(ZWj=0, rfwa = 0, (Zu,^«0, 
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expressed as total differential equations with respect to the 
primitive variables, and the equation 

Of this important proposition I propose to give two distinct 
proofs. 

1st proof. Forming the total differential of the given 
equation we have, on representing its first member bj F^ 

dF 7 dF m dF m 

Xow this cannot be true fof all forms of the function unless 
we have the separate conditions 

du^ =s 0, du^ = 0, du^ = 0. 

Strictly to prove this, suppose i^j, -f^, . . . -F^ to be any n 
distinct and independent fimctions of ^9^a9***t«»9 ana as 
such, distinct and independent forms of F. Then the above 
equation gives 

dF. J , dF^j .dF^j 



dF^jj . dF^jj dF^j 

Now jPj, i^,...jP^ being independent, their ftmctional 
determinant with respect to w^, w„...w^, does not vanish. 
This again is the condition necessary and sufficient that the 
above system^ of linear equations may be independent ; and 
this lastly being the case, their only possible solution mil be 

as was to be shewn. 
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These equations in their developed expression 

da:, * ^ dz ' 

. du^ ^ . du^ , - 

^'^^ +^^^=0' 

enable ns to determine the ratios of (2a?,, dx^^ • . • dx^^ dz in 
the form 

X, "" X, ••• " x^r R ••••• w» 

where X,, -Z,, ...X^, 5 are functions of the original vari- 
ables. And now, forming the equation 

p^dx^ +p^dx^ . . . + p^dx^ — (fe = 0, 

and eliminating the differentials, we find 

X^^+ Xjp^ . . . + X^p^ — B 
for the partial differential equation sought. 

2nd proof. Differentiating the given primitive with respect 
to a?„ as contained explicitly in the functions w, , t«j, . . . w„, and 
also implicitly in the same through z^ we have, on represent- 
ing the first member of the equation by F^ 



dF/du^ du\ , dF (du^ du\ 

duXdxjP^d^J'^di^XdxjP^dz}''' 

dF (du^ du^ ^ 



du^ 



or 



dFdu^dFdu^ dFdu^^dF ^^ 

du^ dx^ du^ dx^ "* du^ dx^ dz ^^ ' 

since 

dFdu^ dFdu^ dF du„ ^ dF 

du^ dz du^ dz "* du^ dz da * * 
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Differentiating thus with respect to the remainijig inde- 
pendent variables, we obtain finally the system 

dFdiA^dF^du^ .dFdu^dF ^ 
du^ dx^ du^ dx^'" du^ dx^ dz^^ ' 

dF du. dF du^ . dF du^ . dF ^ 



duj^ da?, dwj, dx^ 



du^ dx^ da 



dFdu, dl^du^ dF dti^ dF 

du^ dx^ du^ dx^"\ du^dx^ dz^^ ' 

from which, in combination with the equation 



dF du^ dF du^ 
du^ dz du^ dz 

we can eliminate 

dF dF 
du^^ du^^" 

The result will be 



dF du^ dF_ 



dF dF 
du^ ' dz * 



du^ 


du^ 


du^ 


dx^' 


(&,'••• 


'" dx^ 


du^ 
dx,' 


du^ 

dxr' 


du^ 


du^ 
dz' 


du^ 


du^ 
dz 



iift 



p. 



il^n 



,-1 



= 0, 



or, converting rows into columns, 

du^ du^ 



du^ du^ 
dx^ ' dx^ 



dx^ ' dz 



dUn du^ 
dx^' da?a'" 



du^ du^ 

' dx„ ' dz 

p.* -1 



ni 



= a 
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which is the determinant form of the result affirmed in the 
proposition. 

The second of the above forms of demonstration seems to 
be preferable to the first, in that it rests pnlj upon the consi- 
deration of the one general form of the function F. I have, 
however, given the two proofs, chiefly in order to illustrate 
an important remark, vtz. that, in nearly aU general re- 
searches connected with partial differential equations of the 
first order, two modes of procedure, the one involving the 
use of differentials, the other that of differential coefficients, 
may be employed, and that between the forms to which these 
respective modes give rise, a certain law of reciprocity will be 
found to exist. 

The theory of the solution of the partial differential equa- 
tion 

follows immediately firom that of its genesis. If we repre- 
sent yy 

the integrals of the system of ordinary differential equations 
(2) a solution of the given partial differential equation will 
be represented by (1). That this will be also the most gene- 
ral solution may be shewn by the argument of Art. 1. For 
if t£7 = represent any solution, then since 

dw dw - dw . dw 

^/^«^=<^' ^„+^-^=«' 

we find 

^dw ^ ^dw ^dw j.dw 

firom which, in combination with the corresponding equations, 



^ du^ yrdu^ y. du^ j^du^ _ 



B.nE. II. 
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eliminating X, Xj,, ... -Z», JJ we obtain a result which ex- 
presses that the functional determinant of w, w., ... w,» with 
respect to the original variables is virtually 0. Whence w is 
a function of t^^, Wj, ... t«„, and the proposed solution is in- 
cluded in the one to which the above method of solution 
leads. 

That method may therefore be stated in the following Kule. 

B.ULE. To integrate the linear partial differential equation 

^ dz y. dz . yr dz ^ 
^^^ + ^«^/-+^»^„ = ^ 

form the system of ordinary differential equations 

dx^ __ ^2 _ ^n _ ^z 

and deduce their general integrals 

then % 

-F(w„Wj,...w„)=0 

will be the general integral sought. 



[The general observations were intended to follow Chap. 
XIV. Art. 6.] 

General observations. 

4. The relation which exists between a proposed linear 
partial differential equation and its auxiliary system of ordi- 
nary differential equations should be carefully studied. While 
it is prober to say as above that the general integral of the 
one requires the knowledge of all the integrals of the other, 
it is also proper to describe that general integral simply as 
the most general form under which an integral of the auxi- 
liary system can appear. If 

are integrals of that system, then 

F{u^,%,...u^^A 
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is the one general form of an integral of that system, and 
due regard being had to the arbitrariness of F^ this is equi- 
valent to 

-F(2«i, t^j,, ... w») = 0. 

5. The form which the auxiliary system assumes when 
the given partial differential equation is deficient in any of its 
terms should be noticed. 

If X, = 0, the auxiliary equation 

becomes, on clearing of fractions, 

<irj = 0. 
And thus, if Xj, Xj,, ... X^ vanish, the given equation being 

y. dz y. dz . T ^^ — V 

^'dx:''^^'^''d^,'''^^''dx,''^' 

the auxiliary system will be 

dx^ = 0, dx^ = 0,... (&?y = 0, 

dx^^ __ ^'^r+2 _ ^n _ ^^ 

X^^j Xy^j X^ X' 

and the integrals of this system being of the form 



\ 



^r+1 •-- ^r+1 > ^n "~ ^»»> 



the general solution of the given equation will be 

Fix^ Xr, W^+i, Wn) = 0. 

This conclusion would follow also from the principle laid 
down in Chap. xiv. Art. 2. 

Linear partial differential equations in which the absolute 
term is wanting, and which are therefore of the form 

■y^ dz -y dz -y dz 

5— Z 
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may he termed homogeneous. As in this case one of the 
auxiliary equations is 

<fo = 0, 

the general integral will be 

F{u^y%, w^n«) = 0, 

Wj, w,, w^i being found by the integration of the remain- 
ing auxiliary equations 

When X,, X,, -X"„ do not contain z, the solution is best 

exhibited in the form 

« = ^K, w,, O' 

6. Every linear partial differential equation can be converted 
into a homdgenepus one containing one additional variable. 
For it is shewn in Art. 3, that if w = be any integral of 

^ dz ^ dz ^ dz ^ 

then is 

^ du ^ dw du ^du 

^^^,^^^d^r ''<^n'^ dz^ * 

a homogeneous equation with a new variable. 

From the general integral of this equation, that of the 
former one may be deduced by making w = 0. 

7. The solution of partial differential equations is some- 
times facilitated by introducing a new system of independent 
variables. The actual transformation is greatly facilitated 
by the following symbolical theorem. 

Theorjsm. If the partial differential equation 
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be expressed symbolically in the form 
in which 

then, i{ ffiij/^j y* l>® a ^^^ system of independent vari- 
ables given in expression as functions of the old ones, the 
transformed equation will be 

(Ay,)|+(A,o|...+(A,j|.-x 

For, regarding 2; as a function of y^, y,, y«i we have 

dz _ dz^ dy^ dz^ dy^ dz dy^ 

d^.'^dy.dx, ^^"' dy^l^,' 



dz ^ dz dy^ dz dy^ dz dy^^ 

dx^''dy,dx^ dy^ dx^^'^ 'dy^d^J 

whence, substituting in the given equation we find, as the 
total coefficient of -r- , the expression 

or symbolically, Ay^; and so on for the other coefficients. 
The result then is 



It remains only after calculation of Av^, Ay,, Ay„, as 

functions of x^j x^^ x^^ to express these functions and X 

in terms of y„ y„ y^. 

[It appears from the. manuscript that an example was to 
have been supplied here.] 
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[The next Article may be considered supplementary to 
Chap. XIV. Art. 10.] 

Singular Solutions of partial Differential Equations. 

8.. Legendre's theorjr developed in Chap, xxiil. for ordi- 
nary, may be applied also without essential change to partial, 
differential equations.^ Eegarding the independent variable 
z as receiving an infinitesimal change hz through infinitesimal 
change, not in the values of the independent variables 

but in the values of the arbitrary constants of the complete 
or in the forms of the arbitrary functions of the general inte- 
gral, and performing unon the given equation the operation 
denoted by S, we shall obtain a linear partial differential 
equation for determining the general value oihz corresponding 
to any particular given value of z. If that linear equation be 
of a lower order than the differential equation given, then the 
equation expressing the value oi z + hz will be a limiting 
form of a solution less complete or less general than the com- 
plete or general solution of the differential equation given^ 
and the given solution, formed by making the infinitesimal 
constants in the limiting form actually 0, will be singular. 

Conversely, to deduce singular solutions without the know- 
ledge of the complete or the general integral, we ought to 
construct the equations of condition for the reduction of the 
equation determining hz to a lower order than the equation 
given, and the most general solution of the differential equa- 
tions of condition so formed, will be the most general expres- 
sion for the singular solutions of the differential equation 
given. 

Ex. (^w? — gyY s 4- ^wa?' {z — ap) = o, 

in which 

dz 3a 
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Eepresenting the first member of the equation by jP, we 
have, on operating by S, 

dFdSz dFd^z dF^ _ 
dp dx dq dy dz "" ^ 



and the conditions 



^=0 ^=0 
dp ' dq ' 



necessary to reduce the equation for Sz to a lower order give 

{jKc - jy) J — 2mx^ = 0, 
ipx-qy) {px-Bqy) = 0. 

From these we find 

definite and simultaneous values of ^ and j, which being sub-, 
stituted in the given equation lead to 

z = 2m*a?fy*, 

and this, as it gives the same values of p and q as those 
obtained before, will necessarily satisfy the given equation. 
It is therefore a solution, and from the nature of the analysis, 
a singular one. 

Legendre shews that this singular solution is also dedu- 
clble from the general integral of the given partial differen- 
tial equation. That integral is the result of the elimination 
of a from the two equations 

{<l> {a)Y — 2ax<l) (a) + o« — mocy = 0, 

{^ (a) - ax] j> (a) — 2xj> («) + « = 0. 

To deduce the singular solution he supposes j> (a) to be not 
simply a function of a, but a function of a and of one or 
both of the independent variables. He expresses the varia- 
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« 

tion of <^(a) derived from this new source by S, and operating 
on the first equation with S, finds 

{2^ (a) - 2aa?} S^ (a) = ; 

therefore ^ {a) = <ix. 

Substituting this in the equations of the general integral^ and 
eliminating a, we find 

as before. 

Legendre states his theory of the derivation of the singular 
solutions of partial differential equations from the equations 
themselves with great brevity, but still as a general theory. 
And there is nothing in the statement that carries with it any 
apparent restriction upon either the order or the degree of the 
equations given. Until however we are in possession of a 
perfect theory of the genesis of partial differential equations 
we shall not be entitled to say that Legendre's theory of 
their singular solutions is a perfect one; for until then we 
cannot even define, in a perfectly general way, the nature of 
the operation denoted by S. 
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[The next three Chapters all relate to the subject of partial 
differential equations of the first order. The manuscripts do 
not appear to have received their final revision from Professor 
Boole. It is certain that he intended the contents of Chapter, 
XXV. to form a part of the new edition; and it is highly 
probable^ althougn not certain, that the contents of Chapter 
x!xvi. and Chapter xxvii. were also to be included. 

The three Chapters are mainly derived from two memoirs 
by Professor Boole, published in the Phihaophical Trans-- 
actions^ 

The first memoir is entitled On Simultaneous Differential 
Equations of the First Order in which the Number of the 
Variables exceeds hy more than one the Number of the jEqua- 
iions: it occupies pages 437 ...454 of the Philosophical Trans- 
actions for 1862. 

The second memoir is entitled On ths Differential Equa- 
tions of Dynamics. A sequel to a Paper on Simultaneous 
Differential Equations: it occupies pages 485... 501 of the 
Philosophical Transactions for 1863. 

The first memoir was finished before Professor Boole had 
seen Jacobins researches, which are cited at the beginning 
of Chapter XXVI ; these researches indeed could only just 
have been published. In his second memoir Professor Boole 
describes Jacobi's methods, refers to his own already pub- 
lished, and points out the nature of the connexion between 
them.] 
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CHAPTER XXV. 

ON SYSTEMS OP SIMULTANEOUS LINEAR PARTIAL DIPPEREN- 
TIAL EQUATIONS OP THE FIRST ORDER, AND ON ASSO- 
CIATED SYSTEMS OP ORDINARY DIFPERENTIAL EQUATIONS. 

1. The term simultaneous is here applied to a system of 
partial differential equations, to signify that in that system 
there is but one dependent variable, the general expression 
of which, as a function of the independent variables satisfy- 
ing all the equations at once, is the object of search. All 
linear partial differential equations of the first order being re- 
ducible to the homogeneous form, we shall presuppose this 
reduction here. Under this form indeed the problem actually 
presents itself in Geometry, in the theory of partial differential 
equations of the second order, and in Theoretical Dynamics. 

We are sometimes led, in connexion with the same class 
of inquiries, to systems of ordinary differential equations 
marked by the peculiarity that the number of the variables 
exceeds by more than one the number of the equations. Such 
systems are intimately connected with the former'— stand 
to them indeed in a similar relation to that which the 
Lagrangean auxiliary system bears to the single partial dif- 
ferential ecjuation from which it arises. The theory which 
explains this connexion, and grounds upon it the method of 
solution of both systems will form the subject of the present* 
Chapter. 

Connexion of the Systems. 

2. Prop. I. The solution of a system of simultaneous 
linear partial differential equations of the first order may be 
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made to depend upon that of a system of ordinary differential 
equations of the lirst order in which the number of the vari- 
ables exceeds by more than one the number of tlie equations. 

The system of partial differential equations being reduced 
to the homogeneous form, Chap. xxiv. Art. 6, let n be the 

number of the equations, x^, a?,, x^^ the independent 

variables, and P the dependent variable. 

Then from the n given equations determining 

dP 



dP dP 
dx, ' dx^ ' 



dx' 



we obtain an equivalent system of equations which, by trans- 
position of its terms to one side, assumes the reduced form 



dP^ . dP ^j dP 



dx, • ""rfa?«+i * "^^«+9 

dP . dP . dP 

— + ^ +-4 ■ 

dXg ^^n+i ^^n+2 



+ ^K-5 = 



dx, 



n+r 



4 dP ^ 

+ A^- — =0 



"■ar 



dx 



n+r 






(1), 



dx. '^dx.., "dx. 



+ A^4^^0 



'«+l 



M+a 



dx. 



n+r 



Multiplying these equations by the arbitrary constants 

\> \> \ 

respectively, and adding the results, we have 
dP ^ dP .. dP 

'l 2 

dP 



* dx. ' dx„ 



+ 



dx 



•H-l 



dP 



+ (x,Ar+M» +\A^) J-- = 



dx. 



(2), 



n^ 



I 
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a single partial differential equation which, on account of the 

arbitrariness of X,, \,, \, is equivalent to the system 

from which it was formed. 

Of this equation the Lagrangean auxiliary system will be 



\ \ 



^fi 



\,-4„ -f \^ai + \^ 



ni 



^n^ 



\^ir + \^ar + \^m' 



(3), 



W, 



whence, eliminating \, ,X8, X^, we have the system of 

ordinary differential equations 

dx^^ — A^^dx^ — A^^dx^ .. . .--A^^dx^ = 0' 
dx^,-- A.^dx,- A„dx, ....--A^dx^^o 

dx^-AJbn^ - AJoi^ ....-A^dx^ = 

These equations being included in the previous system (3), 
any integrals 

u^a^ t? = i, ti? = c, &c. 

of them will be integrals of it. Therefore w, v, t(?,... will be 
values of P satisfying the partial differential equation (2). 
For they will be the only values which can satisfy it inde- 
pendently of Xj, X,, X^. Hence they will satisfy the 

equivalent system (1), and the general integral of that system 
will be 

F{u, % w, ...) =0 ..........;• (5), 

the form of F being arbitrary. 

Thus the relation of the system (4) to the system (1) is the 
same as the relation of the auxili^iy system of a single linear 
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partial differential equation to that equation. And the ground 
of this relation is seen to be the same in both cases. The 
one form necessitates the other, 

3. Instead of employing the above mode of deducing 
the auxiliary system, we might employ the following which 
is practically more convenient. 

Since any value P which satisfies the partial differential 
equations determines P=c as an integral oi the ordinary sys- 
tem, the latter must be consistent with dF=0 in its de- 
veloped form 

g(fo. + g&,. +^dx^ = o. 



'8 -""^fH-r 



■p,. . ^. dP dP dP 
Ehmmatmg ^, ^, ^^ 

by means of the n given equations (1), we have 
dP 



dF 



dx. 



n+r 



Whence, equating to the respective coefficients of 

dP dP dP 

we have the system (4). 

In the same way we can pass from the system of ordi- 
nary to that of partial differential equations. From the equa-* 
tton dP = 0, in its developed formf we must eliminate a number^ 
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of differentials dlr,, db,, ••• equal to that cf the given equa- 
iionsy and then equate to me coefficients of the remaining 
differentials, 

4 Lastly, the formal connexion of the two systems should 
be noticed. The partial differential equations being given in the 
reduced form (1), the ordinary system may be constructed as 

follows : For any differential coefficient, as -^ — , in any 

column after the first, write the corresponding differential 
dx^^, subtract from this the sum of dx^, dx^, dx^, mul- 
tiplied respectively by the descending coefficients of that 
column, and equate the result to 0. The system of equations 
thus successively formed will be the auxiliary system sought. 

The transition from the ordinary to .the partial system may 
be effected by the same rule, substituting only differentials 
for differential coefficients. 

[It appears from the manuscript that an example was to 
have been supplied here.] 

Up to this point the theory of systems of partial differen- 
tiial equations is in analogy with that of single equations. 
But here a difference arises. We do not know beforehand 
what number of integrals a system of ordinary differential 
equations, in which the number of variables exceeds by more 
than one the number of the equations, admits. 

The theory which removes this diflSculty will be developed 
in the following sections. It will be shewn that a system of 
linear partial differential equations which admits of solution 
by the assigning to the dependent variable a value which 
satisfies all the equations in common, must either itself satisfy 
a certain condition, or be capable of being developed into a 
new but equivalent system which will satisfy that condition. 
It will be shewn that when that condition is satisfied, the 
auxiliary system of ordinary, is capable of expression as a 
system of eocact differential equations determining the inte- 
grals sought 
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It will be found convenient to express by a single symbol 
the aggregate of the operations to which tne dependent vari- 
able is subject in the expression of a partial differential equa- 
tion. Thus the equation 

dz dz dz ^ 
dt dx ^ d2/ 

may be expressed in the form 

Az = 
if we assume 

A — ^ _L ^ _1_ ^ 

dt dx dy* 

Under this convention the following proposition is to be 
imderstood. 

5. Prop. II. If AP = 0, A'P=0 represent any two 
homogeneous linear partial differential equations of the first 
order, then will 

(AA'-A'A)P=0 

also be a homogeneous linear partial differential equation of 
the first order, and it will be satisfied by all the common 
integrals of the equations from which it is derived. 

First, the equation will be linear. For, let a?, y represent 
any two variables whatever, or the same variable repeated, 
out of the set a?! .... a;„, and let A, B represent any functions 
of the variables a?^ .... a?«. Then A may be represented by a 

series of terms of the form A -7- , and A' bv a series of terms 

tic' -^ 

of the form B -j- . Hence (AA' — A' A) P can be expressed 
by a series of terms of the form 



dx\ dy) dy\ dx) ^ 
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which, on effecting the differentiations, becomes 

dx dy dy dx* 

the terms containing the second differential coefficients of P 
mutually destroying each other. Hence the equation 

(AA'-A'A)P=0 

will be a homogeneous linear partial differential equation of 
the first order. 

The constitution of the coefficients of this equation is easily 
determined. For suppose the given equations to be 

dP dP .ji^P_n 

80 that 

then the equation 

(AA'-A'A)P=0 

may be written in the form 

and, since terms involving second differential coefficients of 
P will disappear, this becomes 



+ (A5»-A'^J^ = 0.- 
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We see from this that the f(yrm of the result is the same as 
if the A or A' from either equation operated only on the coeffi" 
dents in the other equation. 

Secondly, the above equation will be satisfied by all the 
common integrals of the equations from which it is derived. 

For, let (^ = c be a common integral of 

AP=0 andAT=0, 
then 

A<^ = 0, A'<^ = 0. 

Performing on these the respective operations A' and A, 
operations which involve only diflferentiation together with 
algebraic processes, we have 

A'A<^ = 0, AA'<^ = 0, 

whence, by subtraction, 

AA'(^ - A'A(^ = 0, 

or ( AA' - A'A) (^ = 0, 

from which it appears that <^ is also an integral of the 
equation 

(AA'-A'A)P=0, 

as was to be shewn. 

6. Prop. III. If by the above processes of reduction and 
derivation we convert a system of partial diflferential equa- 
tions into a new system, such that if expressed in the form 

A,P=0, AjrP=0, A^P=0, 

the condition 

(A,A,-AA)i'=0 

shall for each pair of equations be identically satisfied, then 
the system of ordinary diflferential equations corresponding 
to this new system will admit of reduction to the form of 
exa^t differential equations, the integration of which will 
enable us to construct the general value of P satisfying the 
system given. 

B. D. E. IT. . 6 
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1st. Suppose the given system of n equations reduced to 
the form (1), marked by the peculiarity that n of the differen- 
tial coefficients appear only in successive equations and with 
the coefficient unity. Tnen taking any two of those equa- 
tions (we select the first two), we have 



dx,^-^dx^, -""•<&, 



R+r 



A ^ i A ^ . /I ^ 



ar J^ I 



<^. "«fe.*t "dx, 



n^ 



from the forms of which we see that the derived equation 

(A,A. - A.A.) P = 

cannot contain either 

dP dP 
dx^ dx2 

It can only, as appears from Art. 5, contain the differential 
coefficients 

dP dP 



aiid must be of the form 

It cannot therefore be an algebraic consequence of any of the 
equations of the system (1) from which it was derived. It is, 
unless by the vanishing of B^,....B^ it present itself as an 
identity, a new equation algebraically independent. Com- 
bining this with the former ones, we have a system of w + 1 
equations admitting of the same reduction as to form fol- 
lowed by the same subsequent process of derivation. And 
the result of each of these completed steps is to convert the 
system into one containing one equation more than before ; 
but containing in each of its equations one term fewer than 
before. The process must then end either in the genesis of a 
system of partial differential equations such that the further 
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application of the process of derivation of Prop. Ii. shall only 
lead to identities, or in the emerging of the system 

— = — = — = 

dx^ ' dx^ ' dx^^ 

The latter supposition would imply that P is a constant 
The consequences of the former we proceed to examine. 



The final system of linear partial differential equations 
will be of the same type (1) as the original system, but will 
differ from that system in that n will be increased, and t 
diminished by the same amount. We sliall therefore simply 
state the form (1), only under the condition 

(A,A,-AA)i'=0, 
and with the altered values of n and r. 

First, then, the common integrals of the new system will 
be the same as those of the original system. This is evident 
from Prop. II. 

Secondly. If we write 
the first equation of the system (1) will be 

and the auxiliary Lagrangean system of this will have m — 1 
independent integrals 

among which the n — 1 known integrals (Chap. XXIY. Art. 5) 



^«=^a> ^8 = ^8> »n = ^,| 



6—2 
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are Included. And the general value of P satisfying the 
above first equation will be 

The assumption P=x^ would not satisfy the said equation, 
for It would lead, on substitution, to 1 = 0. Hence we infer 

that while the functions u^, u^^ u^. are independent with 

respect to each other, they are also independent with respect 
to a?i, so that the m functions u^, w^, ...... w^^, ajj, are mu- 
tually independent in the sense explained in Chap. xxiv. 

Let us now transform the equations of the system (1) after 

the first by introducing u^, w^, w^^, x^^ as independent 

variables. Those equations being 

A«i'=0, A,P=0, 

the result of the transformation will be (Chap. xxiv. Art. 7) 

/A ^^P . /A \^P . /A N ^^ . /A ^^P ^ 



But P = a?j. being an integral of each of the equations of 
the system (1) except the first, as appears from their forms, 
we have 

^^1 = 0, A„a?j = 0, 

thus the last terms in the transformed system vanish. Further, 
the coeflBcients of the remaining terms reduce to functions of 
t*j , Wj, ...... . w^j merely. For, considering the coeflScient A^u^ , 

we have 

^ (AA-AA)wi = 0, 
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which, since A^u^ = reduces to 

\A^u^ = 0. 

. Hence A,Wj must be a solution of AiP= 0, and therefore a 

function of w,, u^, w«,_j- And so for the others. It 

results therefore that the transformed system is 

/A \^P fA ^dP /A \ ^P r. 



(^^^) ^ + ^^-^^^ ^ •*'^^""-^) a:::== '' 

ttj, w,, w^j being the oc^tmZ independent variables of the 

system. 

But the transformation having involved no loss of gene- 
rality, for a new system of m independent variables was 
simply substituted for an old one, the condition 

(AA-AA)-P=o, 

satisfied before, will continue to be satisfied in the new sys- 
tem represented symbolically in the form 

A.P=0, A3P=0, A,P=0. 

Any common inte^als of this system will also be common 
integrals of the previous system. For as functions of 

«*!»««,» w^i 

thev will satisfy the first equation of that system, and they 
will satisfy the other equations, because the present system is 
but a transformation of those. The converse is equalfy mani- 
fest. 

Thus a system of n partial differential equations contain- 
ing m independent variables and satisfying a certain condir 
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tion, has in virtue of that condition been converted into a 
system of w — 1 equations between w — 1 independent vari- 
ables, and satisfying the same condition. This then is con- 
vertible into a similarly constituted system of n — 2 equations 
containing w — 2 independent variables, and so on till we 
arrive at a final single partial differential equation containing 
m — n + 1 independent variables. This equation has m — n, 
that is, r integrals, and these are the common integrals of the 
system (1). 

But the system of ordinary differential equations corre- 
sponding to (1) is in number r, and is satisfied by all the 
common integrals of that system. Hence these differential 
equations must admit of reduction to the exact form. 

7. We may deduce from the above investigation the fol- 
lowing Rule. 

To integrate a system of simultaneous linear partial diffe- 
rential equations of the first order. 

EuLE. Eeduce the equations to the homogeneous form 
(I), express the result symbolically by 

A,P = 0, A^=0, A,P = 0, 

and examine whether the condition 

(AA-AA)P=o 

is identically satisfied for every pair of equations of the sys- 
tem. If it be so, the equations of the auxiliary system. 
Prop. L, will be reducible to the exact form, and their inte- 
grals being 

w = a, t; = 5, w = c, 

the complete value of P will be F{u^ v, w, ...), the form of JP 
being arbitrary. 

If the condition be not identically satisfied, its application 
will give rise to one or more new partial differential equa- 
tions. Combine any one of these with the previous reduced 
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system, and again reduce in the same way. Witli the new 
reduced system proceed as before, and continue this method 
of reduction and derivation until either a system of partial 
differential equations arises between every two of which the 
above condition is identically satisfied, or, which is the only 
possible alternative, the system 

^ = ^ = 0... 

appears. In the former case the system of ordinary equations 
corresponding to the final system of partial differential equa- 
tions will admit of reduction to the exact form, and the gene- 
ral value of P will emerge from their integrals as above. In 
the latter case the given system can only be satisfied by sup- 
posing P a constant. 

Ultimately then the determination of P depends on the 
solution of a system of ordinary differential equations reduci- 
ble to the exact form. This does not mean that each equation 
of the system is reducible to the exact form, but that the 
equations may be combined together so as to form an equal 
number of equivalent equations of the exact form. Generally 
when we know this combination to be possible it is easy 
to effect it, and best to endeavour to do so. We might how- 
ever employ the method of the variation of parameters as fol- 
lows. Supposing j> the number of differential equations make 
all but j> + 1 of the variables constant, integrate the reduced 
system, and then seek to satisfy the unreduced system by the 
same series of integrals with the arbitrary constants as new 
variables. The successive integrations and transformations 
of this method would amount to the same thing as those 
upon which the second part of the demonstration of Prop. ill. 
rests*. 

Lastly, given a system of ordinary differential equations 
containing a superfluous number of variables without know- 
ing how many integrals they admit, we must, supposing 
P= c to be any integral, construct the corresponding system 

* It was thus indeed that the author was first led to that theory. 
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of homogeneous partial differential equations satisfied by P, 
and apply to them the foregoing Eule. 

8. Ex. Keqaired tbe Integrals of tlie simoltaneous par- 
tial differential equations 

dP , -.dP , ^dP ^ 

dP^, ,^ .dP , , , .dP . 

Representing these in the form AiP=0, AaP = 0, it will 
be found that the equation 

(AA-AA)P=o' 

becomes, after rejecting an algebraic factor, 

dPdP . 

and the three equations prepared in the manner explained in 
the Rule will be found to be 

f+(3a^4-0f=O, 
dP^ dP . 

dt dz 

No other equations are derivable from these. We conclude 
that there is but one final integral. 

To obtain it, eliminate 

dP^ dP dP 
dx* dy ^ dt 
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from the above system combined with 

dx dy '^ dz dt ' 

dP 

and equate to the coefficient of -r- in the result. We 

find 

dz — {t+ 3a?) dx -r- ydtf — ocdt = 0, 

the integral of which is 



.8 y - 



« — aj^ — a? — ^ = c. 



An arbitrary function of the first member of this equa- 
tion is the general value of P. 

[It appears from the manuscript that another example was 
to have been added here.] 



( 90 ) [CH. XXVI. 



CHAPTER XXVI. 



HOMOGENEOUS SYSTEMS OF LINEAR PARTIAL DIFFERENTIAL 

EQUATIONS. 

1. The theory of homogeneous systems of linear partial 
differential equations in which when expressed in the sym- 
bolic form 

A.P= 0, A.P= 0, A„P = (1), 

the condition 

(AA- AA) P= (2) 

is for all combinations represented by i and j satisfied in 
virtue of the constitution of the symbols A<, Ay, forms the 
subject of important researches by Jacobi {Nova Methodus... 
Crelle's Journal, Vol. LX. p. 1). The following are the most 
important of his results. 

1st. An integral of any one equation of the system being 
found, other integrals of the same system may be obtained 
without integration, by a process of derivation founded upon 
the condition (2). 

Let ^ be an integral of the first equation of the system. 
Then is the equation 

A..^ = 
identically satisfied. 

Also the condition (2) being satisfied in virtue of the con- 
stitution of the symbols, we have 
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and in particulary making t «= 1, and separating the terms, 

which reduces hj a prior equation to 

AjA^ = 0. 

It appears from this that A^, if it do not reduce to a con- 
stant, is an integral of the iirst equation A^^ = 0, and, if it 
prove to be not a mere function of ^, a new integral. 

This process may be repeated upon the new integral with 
a similar alternation of results. It will be evident from this 
that if we confine our attention to the two equations 

and suppose, as before, ^ to be an integral of the first, then 
-wiU 

A..^, A,(A..^), A.{A,(A.^)},... 

or, as these maj be expressed, 

Aa^, Aa'^, Aj'^, ... 

be also integrals of the first equation ; and this process of 
derivation may be continued until we arrive at an integral 
A,**^ which is not independent, but is expressible as a func- 
tion of prior integrals 

and, sooner or later, such a result must present itself, since 
the number of independent integrals is finite. 

It is further seen that the most general symbolic form of an 
integral derivable from the root integral ^ is 

A,*A3^ A,'*<^, 

aj fif fly being positive integers. 

The above remarkable theorem was in some degree antici- 
pated by the researches of Poisson. 
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2ndly. Jacobi shews how by the aid of such derived in- 
tegrals of the first eqttation of the system a common integral 
of the first and second equation may be found, and how from 
this integral and its derived series a common integral of the 
first three equations of the system may be found, and so on, 
until a common integral of the entire system has been as it 
were built up out of previous integrals of less general appli- 
cation. 

Let <^, ^', ^", ^<'*-^> represent a series of independent 

integrals of the equation AjP=0, of which <f> is the root in- 
tegral, and the rest are derived from it by successive applica- 
tions of the operation denoted by A,, so that 

^' = A,^, <^^-i) = A/-i<^; 

also let Aj^^ be not a new integral but a function of 

^,f ^^-«. 

Now ^, <f>\ ^('*"i) being particular integrals of AjP=0, 

the function F{(f>^ ^', ^^'^)) will also be an integral of 

the same equation irrespectively of its form. Let us inquire 
whether the form of the function can be so determined as 
to render it also an integral of the second equation AsP=s 0. 

We have then to satisfy the equation 

\F{<f>,<l>\ <^(^-i))=0. 

By the principles of the Diflferential Calculus this equation 
assumes the form 

A f ^^ . A f r dF A ,1 i\ dF _ 

But A,^ = ^', A,<^' = 4>", Aj^CM-i) ^ ^0* -1) . 

lastly, A,^^-^> may bv hypothesis be expressed in the form 
/(^, <l>', (f)^"^^). Thus the equation to be satisfied is 
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+/(^,.^', ...... ^0.-1))^ = 0, 

a linear partial differential equation of which the auxiliary 
STStem is 

Now the integration of this system may be made to depend 
npon that of an ordinary differential equation of the {ji — 1)"* 
degree between the two variables ^('*"^) and j>. 

For we have 



d^ <f) d<f> <f> 

d<^ (f}' 

Differentiating the last equation with respect to 6, and attend- 

ing to the former ones, we shall be able to express J in , 

terms of the variables ^, ^', ^^-^). Proceeding with 

this in the same way and continuing the process we shall be 
able to express the series of differential coefficients 

d(t> ' e?<^« ' d(l>^-^^ 

in terms of <^, ^', ^^"^>, From these /a — 1 equations, 

eliminating ^', ^", ^^""^j we shall have a final equation 

between 

that is, a differential equation of the (ji — 1)"* order between ^ ' 
and <f)^\ 
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The complete integral of this equation will be of the form 

^^"?=/(*>.Ci><5«, <5m-i)- 

Diflferentiating this /a — 2 times in succession with respect to 
(f), and continually substituting for the differential coefficients 
of (f}^"'^'^ their values as before assigned in terms of 

^, f , ^<^-«, 

we shall have a system of /ifr — 1 equations connecting the 

above variables with the constants c^, Cj, c^-i. Finally, 

solving these equations with respect to the constants, we shall 
possess the integrals required in the form 



J^M-l(*.f r-'') = CM-l. 



and each of these will be a common integral of the first two 
equations of the given system (1). 

[On the back of a page of the manuscript the following 
paragraph occurs, which seems to have been intended as a 
simplification of the preceding argument which begins with 
"The complete integral."] 

Suppose that a first integral of the equation can be found. 
Its form will be 

F[6 A^'*-^^ -5 ^ y ]-c 

Substitute in this for the differential coefficients of ^^'*-^J 
their values before assigned in terms of <f>, ^', <f>',...^^f^\ and 
we have an integral of the system (3), and therefore a com- 
mon integral of the first two equations of the system (1). 

[We now return to the place at which we inserted a para- 
graph.] 

Just in the same way Jacobi deduces a common integral of 
the first three equations of the system (1). For representing 
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any one of the first members of the above system by >|r, and 
deriving thence the new independent integrals Ag>|r, A^yfr,... 
he substitutes an arbitrary function of these for P in the 
equation 

A,P= 0. 

It is evident that the solution of the partial differential 
equation so found will again be reducible to that of an 
ordinary differential equation between two variables. And 
so the process is carried on till all the equations are satis- 
fied. 

2. The above remarkable process was developed by Jacobi 
in connexion with the theory of non-linear partial differential 
equations of the first order. In that particular connexion it 
admits of certain reductions tending to diminish the order of 
the differential equations to be integrated. But these do not 
affect the general principle of the method. It was in this 
special form that the theory of the solution of simultaneous 
Imear partial differential equations originated. Jacobi does 
not consider the theory of equations in which the condition 
(2) is not satisfied ; but the language in which he refers to 
the condition shews that he had speculated upon the general 
problem-^and it is difficult to conceive that he should have 
meditated upon it and not arrived at its complete solution. 

[The manuscript here gives the first two words of 
the passage firom Jacobi's memoir which is quoted in the 
FhilosopKical Transactions for 1863, page 486.] 
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CHAPTER XXVII. 

OP NON-LINEAR PARTIAL DIFFERENTIAL EQUATIONS OP THE 

FIRST ORDER. 

1. In treating the present subject we shall first consider 
that class of non-linear partial differential equations of the 
first order which involves two independent variables, and 
then proceed to the general theory. The reason for this 
procedure is that the particular theory, though of course in- 
cluded in the general one, rests upon a somewhat simpler 
basis, and it was in fact developed by the labours of 
Lagrange and Charpit long before the general theory was 
known. The latter we owe to the independent researches 
of Cauchy and Jacobi. 

[Here the manuscript refers to the matter contained in 
Chap. XIV. Arts. 7 to 12 inclusive; and then passes on to 
the general theory.] 



General Theory. 

2. Given an equation of the form 

« = ^ (a?j, a?2, ... (Tn, «!, ^j,, ... «»)> 

the number of arbitrary constants 6l, a^, ... «« involved being 
equal to the number of the independent variables a?j, a?^, ... Xn, 
we obtain by differentiation and elimination of the constants 
a partial differential equation of the first order. Of this the 
proposed equation is said to constitute a complete primitive. 
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The form of the above process which it seems best, as 
throwing light upon the inverse problem of deducing the 
complete primitive from the partial differential equation, to 
employ, is the following. Let the given primitive, solved 
with respeict to one of the arbitrary constants a^, be presented 
in the form 

/(a?j, ... x^,Zf a^j ... aj —a^^ (1). 

Differentiating with respect to each of the independent vari- 
ables we have a system of n equations of the, forms 

'(2). 

f^ (iTj, ... iCn, Z,Pn, «a, ••• «n) = J 

These n equations enable us first to eliminate the w — 1 

constants a,, a^, and so deduce the partial differential 

equation sought in the form 

F^{x^,... x^,z,p^,...pn)=0 (3); 

secondly to determine the w — 1 constants as functions of 
x^y ... x^, ^9pi9 ... ^« in the forms 



F^ (a?j, ... x^, z,p^, ... pO =«s 



(4). 



As the system formed of these w — 1 equations, together 
with the previous one, is merely another form of the system 
(2) obtained bv directly differentiating the primitive, it follows 
that if from these equations we deduce the values otp^, >.. Pn 
as fiinctions of aj^, ... a?„, a,, ... a^, and substitute them in the 
equation 

dz=^pjdx^+p^dx^+ ...+Pndx^ (5), 

they will render that equation integrable, and its integral 
will be the complete primitive (1), the constant a^ being re- 
gained by integration. 

B.D.E. IT. 7 
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Examining the system (3), (4) we see that the first mem- 
bers of all the equations which it contains are functions of 
ajj, ... aJn> ^> i^i> •••i^*> while the second members are con- 
stants. The question then arises, What mutual connexion 
exists among these functions in virtue of which they yield 
values of p^, ...i?n> which render the equation (5) inte- 
grable? 

The answer to this question must involve the entire theory 
of the solution of partial differential equations of the first 
order, so far as relates to the determination of a complete 
primitive. Given a partial differential equation of the form 
(3) it is evident that if we can construct a system of associated 
equations (4) possessing the character above described, the 
final value of z obtained by integration of (5) will both 
satisfy the given equation and contain the requisite number 
of arbitrary constants. It does not follow tirom this that 
it will be the only complete primitive, but it will be a 
complete primitive. 

3. The relation sought is expressed in the following 
Proposition : 

Proposition. If 

<E>(a?i, ... x^yZ,p^, ...p^=b 

represent any two out of a system of n independent equations 
such that the values ofp^^ ••• Pn> thence determined would make 
the equation 

dz =]pjdx^ + p^dx^ + ... +j[>nidxn 

integrablef then the first members of these eqiuztions being 
represented for simplicity by F and <E>, the condition 

^ (fdF , dF\ d^ dF (d^ d<bS\ ^ 

the sumrMition extending to all values of i^frcm, Ifon inclusive^ 
ivill be satisfied identically. 
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Reciprocally^ if the above condition he satisfied identically 
for each binary combination of Junctions in the proposed system 
of eqimtions, and if these Junctions be independent, then the 
values ofp^ ,.,. p^, as functions ofx^^ ...x^, z, which they yield, 
will make the equation 



dz=pj^dXj^+p^dx^ + ... +Pndx^ 



integrabjs. 



It will be convenient to begin with the particular case in 
which the proposed equations do not explicitly contain z, the 
particular pair to be considered being represented hj 

Jf (x^ , »,.x^, p^, ... Pf^) = a, 
^{x^,...x^,p^,...p^) = b. 



Differentiating with respect to Xi, and regarding ^j, ...p» as 
functions of the independent variables, we have 



dF dF^ dF dp^ 

dxi dp^ dxi "* dpn, dxi 

d^ d^ dp^, 4,^ ^ 

dXi dp^ dXi *" dp^ dxi 



=0, 



= 0, 



(6), 



to which we may give the form 

dF ^ dF dpj 

^ - _ V ^^ ^Pi 
dXi "" ^ dpj dXi^ 



(7), 



the summation with respect toy extending from j=l to J=n 
inclusive. 

d^ 
From the first of equations (7) multiplied by -j- subtract- 

1— 'L 
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the second multiplied by -7- , and sum the result with respect 
to i from t = 1 to t = w inclusive. We have 



* \dXi dpi dpi dxj 



^^tt (—— ^^— — ^^ (^) 

* ^\d][)j dpi dxi dpi dpj dxJ 

The expression under the double sign of summation in the 
second member vanishes when t=y; we may therefore re- 
strict the summation to unequal values of i and j. Now 
as for any particular combination of values, e.g. 2, 3, there 
would exist in the completed member both the terms cor- 
tesponding to « = 2, J = 3, and those corresponding to y = 2, 
t = 3, it is evident that if we employ the symbol S^y to denote 
summation with respect to different combinations of i andj, 
the second member of the last equation may be expressed in 
the form 

^ fdF^ c?<t> dpj dF d^ dpj 
''^Kdpj dpi dxi dpi dpj dxi 

dFd^^^dFd^ dpA 
dpi dpi dXj dp^ dpi dxjj ' 



+ 



or 



S f/^^_^^W^i_^M 



\\dpi dp J dp J dp 
so that the equation (8) becomes 

* \d^i dpi dpi dxJ 



^[Kdpidpi dpjdpj\dxj dxJ] ^ ^' 

The number of terms of which the second member ex- 
presses the sum is thus — ^- — - , and it will be observed that 
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as to any particular term it makes no difference in what order 
the numerical values of i and j are assigned to these quan- 
tities; e.g. whether for the combination 2, 3 we make /=2, 
y==3, or « = 3, j = 2; but we must confine ourselves to one 
order. 

Now when the equation 

is integrable in the manner here supposed, we have for all 
combinations of t andy, 

dxj dxi* 

All the terms in the second member of (9) therefore vanish, 
and we have 



2 (^(^_dFd^\ 
* \dxi dpi dpi dxj 



This is the direct form of the Proposition under the parti- 
cular limitation supposed. 

As F, ^ represent, under the same limitation, any two of 
the first members of the n equations (3), (4), which determine 

Pit »..pny there will exist ^ equations like the above. 
It is usual to employ for brevity the notation 

/dF d^ dF d^' 



2 ((^d^_dFd^\^.„^ 

^\dxi dpi dpi dxJ *■ ^* 



and this being done the above system of equations expresses 
the — ^^ — -functions of the form [F^F^ as linear homogeneoits 

functions of the — ^-^ — - quantities of the form ^ — -r-^ . 

It is hence that the vanishing of the latter series of quantities 
secures the vanishing of the former. 

The converse truth will therefore be established! by shewing 
that the -^-^ — - quantities of the form -J-^ — -^ are, when 
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F^y F^y ... F^ axe independent with respect to p^^ j?,, ...^„, 

expressible as linear homogeneous functions of the — - 

functions \FiF^. 

To avoid complexity of expression I shall establish this for 
the particular case of n = 3, and shall shew that the reasoning 
is general. 

The functions F^^F^^F^^ being independent with respect 
*^Pi> P^^Psi ^^^ determinant 



dF^ 


dF, 


dF, 


dp,' 


dp,' 


dp. 


dF, 


dF, 


dF, 


dp,' 


dp,' 


dp. 


dF, 


dF, 


dF, 


dp,' 


dp,' 


dp. 



does not vanish. This determinant we shall denote by A. 

In (9) writing for i^and * first F and jPg, secondly F^ and 
i^j, thirdly F^ and jF,, we have on changing signs the system 

^ "^"^ \dp^ dp^ dp^ dpj\dx^ dxj 

* *^ \dpt dp, dp^ dpj\dx^ dxJ '" J 

Multiply the first equation by --r-^ , the second by -v-- , the 
third by -j-^ and add. Then 
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whence as A does not vanish we have, on dividing "by it, the 

n/n (IT) 

function ^ "*nr^ expressed as a linear homogeneous function 
of [JTFJ, tp;Fj' and [i^,i^J. 

In like manner multiplying the equations by j-^ , -r-^ , -r-^ 

Pit Pi ■ Pa 

respectively, and dividing by A, we obtain -^ — P^ as a 
similar linear homogeneous function, iand lastly, multiplying 
by ^, ^-^, ^-^, and proceeding as before, we obtain 

-P^ — -j-^ as a similar linear homogeneous function. 

From all which it follows that when [FJF^, K^Jj [^1^2] 
vanish, then 



^^^Pz ^^^Pi iPx_^ 
dx^ dx^^ dx^ cfccj' djc^ dx^ 

will vanish also. 



The reasoning is general in its nature. If i^,, F^, ... jP„ 
are independent with regard topjjjp^, ...p^, the determinant 



dt\ dF, 



dF. 



dF^ 



dp^ ' dp. 



= A 



(11), 



does not vanish. This deterniinaiit is from its constitution 
as a determinant linear and homogeneous, not only with 
respect to any row or column of elements, but also with 
respect to the possible binary combinations which can be 
formed of two rows or columns, ternary out of three rows or 
columns, &c. provided that these combinations are themselves 
of the form ot determinants. In the language of the theorjr 
such combinations are called minor determinants. Hence if 
we construct liie system of equations teptesented by (10), and 
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observe that the coefficients of any particular term of the form 

-^ — -^ in the several equations form a system of such 

minors to the general determinant (11), it will be plain that 
the equations can by multiplication and addition be brought 
to a form in which the coefficient of that particular term will 
be A. At the same time the coefficients of all the other 

terms of the form -^ — p will vanish. For a little atten- 

tion will shew that they will be what the determinant A 
would become on making two of its columns or rows of ele- 
ments equal, and therefore will be identically equal to 0. 

Thus the Proposition is eenerally established for the case 
in which z does not explicitly appear in the functions 

^19 ^iJ ^n- 

When z does appear in those functions the equations (6) 
will be replaced by 

iE+„^+€:±^+ +^l^=o 

dxi ^* dz dp^ dXi djpn dxi ' 

^ + p ^ + ^& + .....+ ^^« = 
dxi ^* dz dpi dXi dp^ dx^ ' 

from which it is seen that the theorem above established will 
only need to be changed into the form employed in the state- 
ment of the general Proposition, 

As. the above is one of the most important propositions in 
the entire theory of Differential Equations, it may be desire- 
able to illustrate it by examples. 

[There are no examples in the manuscript.] 

4. We resume the general theory. 

The integration of non-linear partial differential equations 
may be effected by two distinct methods, both resting upon 
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the ground of the above PropositioD. The first of these 
methods, originally established hy a different analysis from 
that which will here be employed, was discovered by Cauchy 
{Exet^cices d^ Analyse), and rediscovered by Jacobi {Crelle's 
Journal). The second method, discovered by Jacobi at a 
later period, forms the subject of his posthumous memoir^ 
Nova MeihodiLs 



Cauchy* 8 Method. 

We will, as before, begin with the case in which z does 
not appear explicitly in the proposed partial differential equa- 
tion, which we shall represent in the form 

•^i(a?i, a?„, ^„ i?0 = (1). 

We have seen that to find a complete primitive of the 
equation it is necessary and sufficient to construct a series of 
equations 

-^2(^1* a?»» Pi^ P^-<^^ 

• (2), 

i^,(a?j.....a:«, jp,, i>^) = aj 

such that not only shall the conditions 

[F,F2 = Ki^J = (3), 

connecting the new functions jP,, jP^ with F^, be identi- 
cally satisfied^ but also the series of conditions 

[^«i^J = (4), 

Fa and Fj, representing any two of the new functions re- 
ferred to. 

The first of the above series of conditions amounts to 

this, that F^, F^ must be integrals of the partial different 

tial equation 

Ki^=0 (5). 
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It IS the peculiar aim of Cauchy's method to determine the 
integrals so as to cause the second series of conditions to be 
satisfied also. And it is shewn that this will be attained if 
the integrals of (5), which form the first members of (2), are 
such that the particular values which o^, ....^n_i assume when 
x^ is made to receive any constant value, as 0, are difierential 

coefficients with respect to aj^, x^^ of any single function 

of those variables, the form of which may be arbitrarily 
assigned. 

The necessity of this condition is obvious. If the general 

values of ^j, p^ are differential coefficients of a function z 

with respect to a?^, a?„, then the particular forms which 

j9^, p,^ assume when x^ receives any constant value are 

simply differential coefficients with respect to a?., x^^ of 

what z becomes under the same circumstances. To prove its 
sufficiency we must shew that when it is satisfied the condi- 
tions represented by (4) will be satisfied also. 

Since F^ and Fj, are integrals of [F^P] = 0, 

[F,F,] = 0, [F,F,]^0 (6). 

Also, since if in (I) and (2) we give to x„ a particular con- 
stant value, as 0, and then in (2) regard ^„ as a function of 

^l> ^n-iJ JPi5 Pn^l 

determined by (1), the system (2) will virtually contain only 

of which »j, j?^j are differential coefficients of a single 

function with respect to a?j, aj^j, it follows firom the pro- 
position of Art. 3, that any two fanctions F^ and F^ will 
satisfy mutually the condition 

2i=n.i (dF, dj\ _ dF, dF,\ ^ ^ 
**^^ \dxi dpi dpi dxi) ' 

the differentiations having reference to 
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explidtlj as they appear in F^ and Fj,^ and implicitly as 
involved injp^. Thus the developed form of the above equa- 
tion is 



i=n-l 



St=n 
i=l 



\dxi dp^ dxj\dpi^ djp^ dp J 

\dp^ dp^ dp J \dxi dp^ dxj) ' 
the forms of -p^ and -^ being determined from (1). 

Performing the multiplications, the above equations will be 
reduced to the form 

^"^ \dxi dpi dpi dXiJ 

^ dK 2i=n-i (dpn dF[_^ dF]\ 
dp^ *"^ \dxi dpi dpi dxJ 

_ dl\ ^i=^.i fdp^ dFa _ dp^ dF^ \ ^q /^n 
dp„ *^^ \dx4 dpi dpi dxi) ""^ 

But from the form of the total differential of (1) we see that 

dF[ dF[ 

dpn__dxi dp^_ dpi 

dx^ TF^' dp^ dF\' 

dpn dp^ 

Hence , 

2<=«-i f^n dF^ _ dp^ dFA 
*~^ \dxi dpi dpi dxJ 

\dpj ^"^^ \dXi dpi dpi dxJ* 
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Now since by Art. 3 
we have 



*"^ \dxi, dpi dpi dxJ 



fdF^ dF, dFi dF, 
"" \dx^ dp^ dp^ d^J 

therefore 



^■. 



^*^^ \dXi dpi dpi dxJ 



- /'^V' '^^^' ^^« ^"^^ ^^"' 



( 



\dpj \dx^ dp^ dp^ dx. 



In the same way 



^i=n^i(dPndF\yp^dI\\ 
^^-^ \d>Xi dpi dpi dxJ 



_ (dFX-'/dF, dFl_dF\ dFl\ 
[dpj \dx^ dp^ dp^ dxJ ' 



The substitution of these values in (7) gives 



^i=n^i(dK dFl^dF^ dFf\ 
*'^ \dxi dpi dpi dxJ 



fdFX' m. fdF^ dF, _ dF\ dF;^ 
\dpj \dp^ \dx^ dp^ dp^ dxJ 

dpn W« dp^ dp^ dxJ) ' 
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or 



or 



*'^ \dxi dpi dpx dxj 

dF^ dt\^dJE\ ^*^Q 

^,^,fdKdF[_dK dF\\^ 
^=^\dxi dpi dpi dxj ' 



which is precisely the equation 

We see therefore that to solve the partial differential 
equation 

■^i(^, ^n, JPp i>«)='0, 

it is only necessary to construct the linear partial differential 
equation 



*~^\dxi dpi dpi dxJ ' 



and to obtain w — 1 independent integrals of this 

such that if we determine from these conjoined with the given 

equation the values of p,,....^„, then those of^^, jo»^- 

shall, when x^ is made constant, be the partial differential 

coefficients of one and the same function ot a;^, x^^ with 

respect to these variables in succession. 

Now provided that we can find all the integrals of the 
above partial differential equation the particular determination 
required may be effected in the following manner. 

The Lagrangean auxiliary system consists of 2w — 1 ordi- 
nary differential equations 

^1 dx^ dp, _dp^ .gv 

^[t\ '"I^""^t" ^M 

dp^ dp^ dx^ dx^ 
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These admit of 2w — 1 Integrals, one of which will be i^j = c^ ; 
and this will agree with the given equation if we make c^ = 0. 
We have therefore, besides the particular integral jPi = 0, 
271 — 2 integrals of the form 

(9). 

Now suppose it required to. find a value of « as a function 
of ajj, ...aj„, which shall satisfy the given partial differential 
equation, and shall reduce when a;„ = (or any numerical con- 
stant) to a particular given ftmction of a?^, ...aj^i, which we 
will represent by -^(ajj, ...»«_ J. Then on the assumption 
that x^ = 0, we have first the given equation 

« = i/r(a:„...a!«.J (10), 

secondly the derived equations 






_ d^lr(x^,...x^_^) 
^"-^ - dx. 



(11). 



n-i 



Make in the 2n — 1 integrals a?„ = 0, and suppose at the 
same time a?j, ... a?^^, p^, •••i?n-i to receive therem the same 
values as in the above derived equations. Then firom the 
3/1 — 2 particular equations which we thus possess in the two 
systems united (particular because under the assumption that 
^n = ^)j "^6 can eliminate the 2n — 1 particular values of 
x^, ...a;^i, Pjr*»Pn^iy and so obtain n — 1 equations among the 
constants. These express the conditions which are necessary 
and suflScient in order that the values of^, , ... p„.i thus derived 
from the integral equations may, when aj„= 0, aeree with the 
values assigned in (11). Accordingly if we substitute in these 
equations of condition for Cj, . . x^^^ the general values^^^, . . .^t«-i 
we shall obtain n— 1 equations between aj^, ... a?^, p^, >"Pnj 
which will at once be particular integrals of the system (8), 



ART. 4.] BQCTATIONS OP THE FIRST ORDER. Ill 

and will possess the property that the values ofp^,.,. p^^ which 
they in conjunction with -F, = give will when aj„ = reduce 
io the values given in (11). Hence these values with that of 
p^ derived from the same equation will make 

dz^p^dx^^ —p^dx^^O 

an exact differential equation. In the integral of this it will 
only remain to determine the constant so as to make the value 
of z agree with that given in (lO), All the conditions will 
then be satisfied. 

We may collect the results of the above investigation into 
the following Kule : 

To obtain an expression for 2; as a ftmction of the inde- 
pendent variables a?^, ...aj„, which shall satisfy the partial 
differential equation 

F{x^, ... x^j jpi, ...p^) = 0, 

and shall when a;^ is made equal to (or to any numerical 
value) reduce to a given function of aj^, ... x^^, which we will 
represent by ^jr {x^, ... aj^J. 

Rule. Construct the linear partial differential equation 



^,_^/dF dP dF dP\ ^ 

v*-*^f 1 = 

^^^\dxi dpi dpi dxj ' 



and forming its auxiliary Lagrangean system deduce its in- 
tegrals 

in addition to the known particular integral jF= 0, 

Between the above integrals and the equations 

_ d^(a?t,...a;^t) ^ _ d'^{x^,...x^ ^ 

^^"" d^, "•• -P-^"" dx^, "' 

eliminate, after making a?, = 0, the quantities 

« 

ajj, ...a*„_j, p^y ••*Pn^i* 
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In the resulting n — 1 equations replace 

and we shall have a system of equations which with F=0 
will determine values of j?i, ...jp«, which will render 

dz ^p^dx^ — ... ^Pn^n 

an exact differential. The integration of this will give the 
integral sought. 

In the case in which the given partial differential equation 
is of the form 

F{x^,...x^y z, jp,, ...^J = 0, 
z being contained explicitly, the linear equation to be solved is 
^•^^(fdF dF\dP dF/dP dP\[ ^ 

^-^\[d^^p'-d-z)d^r^^^ 

and the argument by which it is shewn that the integrals of. 
this to be employed in conjunction with F= for the deter- 
mination of ;t7j, ... p„ need only be so conditioned as to make 
i^i5 •••2^n-i differential coefficients of one and the same function 
of ^j,...a^^j when ic„=0 is in character the same as that 
already developed in the present Article. It is only necessary 

to substitute in its exposition -v- +jpi -j- for -,- , and so 

for the other functions. 

But as the auxiliary system 

dx^ dx^ dz 



dp, dp, P^dp, - -P-^ 

dx,^'dz dxj^^^dz 

r 

virtually includes the equation 

dz—p^dx^— ... —Pn<^Xn = ^f 

the ultimate expression of the Bule will be as follows i 
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To obtain an expression for ^5 as a function of x^, ... x„ 
which shall satisfy the equation 

F{x^, ... x^, z, Pj, ... p^) = 0, 

and shall when x„ is made equal to (or to any particular 
constant value) reduce to a given function '^ (a?j, ... x^^ of 
the independent variables a?^, ... x^^. 

Rule. Let 

be the 2w— 1 integrals of the auxiliary system (12) which are 
additional to the particular integral F= 0. Make in these 2n 
equations a?n = and forming the further equations 

z = ylr{x^, ... x^^), 

_ dyl^(x^, ... ar^J 
^^"" dx. 



^ dy^jx^, ... a?«_i) 



eliminate the 2n quantities 0?^, ... a?^_i, z, p^^ ... j?^. We thus 
obtain n equations among the constants c^^»'*c^* 

■ 

Substitute in these equations ^^ for c^, ... ^j^ for Cjn, and 
we have n equations connecting x^, ... x^^, z, p^f.,.p^, from 
which with the aid of the given eq^uation p^, .*. pn may be 
eliminated, and there will result a single equation connecting 
ajj, ... x^ with z. This is the integral sought. 

[It appears from the manuscript that an example was to 
have been supplied here.] 

5. Cauchy's method is evidently a general one. But its 
generality is not of the same kind as that which belongs 
to Lagrange's solution of linear partial diflFerential equations* 
It conducts us, not to a form embracing every possible 

B.D.E. II. & 
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solution, but to a system of results from which every possible 
solution maybe derived, by arbitrarily varying the form of 
the function which expresses the initial state of the dependent 
variable, that is the value of z when x^ = 0, and then per- 
forming certain eliminations. To obtain a complete primi- 
tive we should only have to assume as the form of z when 
a?n = a function of the variables x, ... x^^ involving n 
independent constants. The form of this function is arbitrary. 
Each distinct determination of it under the conditions leads 
to a distinct complete primitive. The number of such com- 
plete primitives is infinite. 

There are some most important problems in which the 
knowledge of a single complete primitive is all that is re- 
quired. For this purpose the method of Jacobi which we 
shall now give may be employed. 



Jacobt*8 Last Method. 

6. Supposing z to be not explicitly involved in the given 
partial differential equation 

i^,(a?„ ... Xn,jp,,pn)^0, 

9 

which we shall as before represent by F^ == 0, the problem of 
the discovery of a complete primitive consists in the finding 
of n — 1 equations 

such that between any two functions FiFj the relation 

lFiF;\ = (1) 

shall be identically satisfied. The values of jp^ ,.../?„ deduced 
from the equations, by rendering 

dz —p^x^ — ... — jpn^'^«= 

integrable lead us to the complete primitive expressed by its 
integral. 
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Now the idea upon which Jacobi's later methods rest is 
that of directly solving the different systems of linear partial 
differential equations flowing from the general condition (1), 
not of solving, as in Cauchy's method, one of those equations 
and then limiting that solution by conditions which virtually 
involve the satisfaction of the others. 

It is evident that the entire series of — ^^ conditions 

(1) will be satisfied if we determine F^ to satisfy the single 
equation 

then F^ to satisfy the system of two simultaneous partial 
differential equations 

then F^ to satisfy the system of three simultaneous partial 
differential equations 

and so on, until finally F^ is determined by the solution of 
the system of n — 1 partial differential equations 

iF,F;[ = 0, \F^;\ = 0, [f^,f;\ = o. 

Now all these are particular cases of the general problem 
of determining a function P which shall satisfy simultaneously 
the equations 

[F,F] = 0, [F,P] = 0, [F,P] = 0. (2) 

F^, jF„ ... F^ being given functions between each pair of 
which the equation 

= 



is identically satisfied. Here P will represent in succession 



the series F^^ F^y ... F^. 



The given system is one of homogeneous linear partial 
differential equations. It belongs to the class of systems the 
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general theory of which is discussed in Chap. xxvi. But it 
is not necessary to apply the theory in its general form. We 
need only a single integral ; for a single value of each of the 
functions i^j, JP3, ... i^„ suffices in combination with the given 
value of F^ for the determination of a complete primitive. 
Now it may be shewn that the system is 01 the class dis- 
cussed in Chapter xxvi. If expressed symbolically in the 
form 

A,P=0, A,P=0,...A^P = 0, 

the condition 

(AA-AA)P=0, 

will be identically satisfied. Hence Jacobi's method for the 
treatment of systems of this kind may be applied. 

That the system is of the kind asserted is a consequence of 
the following proposition. 

Proposition. If the equations 

[wP] = 0, [vP] = 

are expressed in the symbolic form 

AP==0, A'P=0, . 

then the derived equation 

(AA'-A'A)P=0 (3), 

will be equivalent to 



rMPi=o. 

For A = 2{:rf|?^;^-|?i/), 

A' = S'"**/'— — -— — ^ 
*"^ \dxi dpi dpi dxJ * 
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Hence since -r- is the coefficient of -j- in A'P, and -7— 

dXj ^ dpj djL'j 

its coefficient in AP, its co efficient in the derivedequation (3) 

will be (Chap. xxv. Art. 5), 

. dv . , da ^ 
dxj dx/ 



or 



^i^nfdu d^v du d\ dv d^u dv d^u \ 
\dxi dpidxj dpi dxidxj dxi dp^dxj dpi dxidxjj' 



or 



dx 



^i^nfdu dv^ du dv\ 
:j *"^ \dxi dpi dpi dxJ ' 



In like manner the coefficient of -^— is 

dXj 

- 1; M- 

Hence (AA' - A'A) P = Sft? f ^ dP_dJuv\ ^\ 

' ^ \ dXj dpj dpj dxjj 

= [m p\ , 

whence the Proposition is established. 

Applying this to the system (2) we see that any derived 
equation will be of the form 

But \FiF^ = by the conditions given ; hence the condi- 
tion (A^Ay— A^Ai) P = 0, is identically satisfied. 

The results of Chapter xxvi. being thus directly applicable 
to the system under consideration, we see that a common 
integral of the system (2) may be found by a series of alter- 
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nate processes of integration and derivation. We begin by 
seeking an integral of the first partial differential equation. 
B7 a process of derivatian, always possible, followed by the 
integration of a differential equation between two variables, 
we arrive at a common integral of the first two partial diffe- 
rential equations. Again, by a process of derivation followed 
by the solution of a differential equation we obtain a common 
integral of the first three partial differential equations. And 
so on, until a common integral of all is obtained. 

7. Another solution of the above problem has recently 
been given. Beginning as in Jacobi's method by finding an 
integral of the first partial differential equation, a process of 
derivation agreeing in principle with Jacobi's, only more 
extended, may lead us without further integration to a point 
at which the discovery of a common integral of the entire 
system will depend only upon the solution of a single diffe- 
rential equation of the first order susceptible of being made 
integrable by a factor. Failing this, it will enable us to 
convert the given system of partial differential equations into 
a new system possessing the same general character, but con- 
taining one equation less. Upon tliis the same process may 
be tried with a similar final alternative — and so on till the 
required integral is discovered. {On the Differential Equa- 
tions of Dynamics. PhihsopMcal Transactions^ 1863). 



( "9 ) 



CHAPTEE XXVIII. 



PARTIAL DIFFERENTIAL EQUATIONS OF THE SECOND ORDER. 



[This Chapter is a reconstruction on a larger scale of part 
of Chapter XV. At the end of the Chapter reference will 
be given to other writings of Professor Boole on the subject 
here discussed.] 

1. The general form of a partial differential equation of 
the second order is 

F{x, y, z,p, q, r, «, ^ = (1), 

where 

^"■^' ^"Ty' ''"^' ^^db^' d^^ 

It is only in particular cases that the equation admits of 
integration, and the most important is that in which the 
differential coefficients of the second order present them- 
selves only in the first degree ; the equation thus assuming 
the form 

Rr-\-S8+Tt^ V. (2), 

in which By 8^ T^ and Fare functions of ar, y, z^p and q. 

The most important part of the theory of the solution of 
this equation is due to Monge, and was extended by Ampfere 
to the more general equation 

J&- + &4-in5-l-Z7(«*-r0=*r (3). 
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nate processes of integration and derivation. We begin by 
seeking an integral of the first partial differential equation. 
B7 a process of derivatian, always possible, followed by the 
integration of a differential equation between two variables, 
we arrive at a common integral of the first two partial diffe- 
rential equations. Again, by a process of derivation followed 
by the solution of a differential equation we obtain a common 
integral of the first three partial differential equations. And 
so on, until a common integral of all is obtained. 

7. Another solution of the above problem has recently 
been given. Beginning as in Jacobi's method by finding an 
integral of the first partial differential equation, a process of 
derivation agreeing in principle with Jacobi's, only more 
extended, may lead us without further integration to a point 
at which the discovery of a common integral of the entire 
system will depend only upon the solution of a single diffe- 
rential equation of the first order susceptible of being made 
integrable by a factor. Failing this, it will enable us to 
convert the given system of partial differential equations into 
a new system possessing the same general character, but con- 
taining one equation less. Uppn this the same process may 
be tried with a similar final alternative — and so on till the 
required integral is discovered. {On, the Differential Equa- 
tions of Dynamics. Philosophical Transactions^ 1863). 
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CHAPTER XXVIII. 



PARTIAL DIFFERENTIAL EQUATIONS OF THE SECOND ORDER. 



[This Chapter is a reconstruction on a larger scale of part 
of Chapter xv. At the end of the Chapter reference will 
be given to other writings of Professor Boole on the subject 
here discussed.] 

1. The general form of a partial differential equation of 
the second order is 

F{x^ y, e,p, J, r, e, <) «0 ..*.. .(1), 

where 

da ^dz d^z _ d^z ^__^^ 

^"^' ^'^dy' ''^d^' *^^^' df* 

It is only in particular cases that the equation admits of 
integration, and the most important is that in which the 
differential coefficients of the second order present them- 
selves only in the first degree ; the equation thus assuming 
the form 

5r + &+ 2i= V. (2), 

in which -B, 8^ jP, and V are functions of a?, y, «, p and j. 

The most important part of the theory of the solution of 
this equation is due to Monge, and was extended by Ampfete 
to the more general equation 

Br-\-88^Ti^U{s^--rt)^V (3). 
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This equation, together with the particular equation of 
Monge, and the equation 

both which though falling under Ampere's general form 
possess peculiarities demanding special notice, 1 propose to 
consider in this Chapter. I shall in conclusion make some 
observations on the theory of partial differential equations of 
the second order with more than two independent variables. 

Monge's method, and Ampfere's in so far as it is an exten- 
sion of Monge's, consists in a certain procedure for discovering 
either one or two first integrals of the form 

«=/W (4), 

u and V being determinate functions of x, y, «, p, and q ; and 
f being an arbitrary functional symbol, i rom these first in- 
tegrals, singly or in combination, the second integral involving 
two arbitrary functions is obtained by a subsequent inte- 
gration. 

Now this procedure involves the assumption that the pro- 
posed equation admits of a first integral of the form (4). But 
such is not always the case. There exist primitive equations 
involving two arbitrary functions, from which by proceeding 
to a second differentiation both functions may be eliminated 
and an equation of the form (2) obtained, but from which it 
is impossible to eliminate one function only so as to lead to an 
intermediate equation of the form (4). Especially this hap- 
pens if the primitive involve an arbitrary function and its 
derived function together. Thus the primitive 

-3^ = ^ (y + ») + '^ (y - ») - » {0' (y + a?) ^ '^' (y - a?) } . . . (5) , 

leads to the partial differential equation of the second order 

r-t = '-S (6), 

X 

but not through an intermediate equation of the form (4). 

It is necessary therefore, not only to consider the case in 
which the assumed condition is satisfied, but also to notice 
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what has been done in those cases which do not at present 
fall under the dominion of anj known method. 



Oenesis of the Equation^ 

2. Prop. i. A partial differential eqvxition of the first 
order of the form u =f(v), or its symmetrical equivalent, 

in which u and v are any functions of a?, y, «, p^ y, always 
leads to a partial differential equation of the form 

Br-hSs+ Tt+Uis^-rt) = F. 

For, differentiating the proposed first integral with respect to 
X, and with respect to y, we have 



dF/du du du , ^ ^ 
du \dx dz^ dp dq I 



dFfdv dv dv ,^y^/^ 
dv \dx dz^ dp dq J" ^ 



dF fdu du du ,du\ 
du \dy dz^ dp dq J 



dF/dv , ^^ , ^ ^^ /A _ A 
dv \du dz^ dp dq J 



^dy dz^ dp dq 

For brevity, write 

fdu\ ^ du du , fdu\ j, du ^ du 

and then eliminating 

dF dF ^ 

. rfa ' dv^ 



+ 



s 
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we have 

ydu\ du du ] {fdv\ , dv , dv } 

{fdu\ du du) {fdv\ dv dv ] 

which, on eflfecting thie multiplication, gives 

(du fdv\ ^ (di£\ dv\ 
{dp \dy) \dyj dp] 

(fdu\ dv_^du rdv\ _ rdu\ ^ du /dv\\ 
\\dx)^ dp\dx) \dy)dq dq\dt/J) 

fdu dv du dv\ , , ^ 
\dq dp dp dqj ^ 

=(|)(l)-(l:)(|) •••('). 

a result which, since u and v are by hypothesis given func- 
tions of X, y, z, p, J, is seen to be a particular case of the 
general form (3). 

We may hence deduce also the conditions under which 
particular forms included in the general form (3) arise. Thus, 
m order that the equation u =f{v) may give rise to a par- 
tial differential equation of the second order of Mongers form 

Er + Sai-Tt^-V, 

it is necessary that the condition 

du dv ^w ^^_/^ 
dq dp dp dq 

should be identically satisfied. This requires, by Chap. ii. 
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Art 1, that u and v, considered as functions ofp and q^ should 
not be independent. 

3. The geometrical relations of the equation (3) are also 
remarkable. It may in particular be shewn that an equation 
of this form will be satisfied by the equation of any surface 
which constitutes the envelope of any system of surfaces 
formed by the variation of three parameters in subjection 
to two arbitrary conditions. For let the common equation of 
the enveloped surfaces be 

«=/(^>3^j«j^^) W> 

the parameters a, J, c varying in subjection to the conditions 

conditions which, determining b and c as functions of a, may 
be reduced to the form 

b^^{a), c = ylr{a) (9). 

Now the values of p and q being the same for any point 
in the envelope as for the same point in the generating surface, 
we have for all such points 

dfipc^y^a^hcl df{x,y,a,h,c) .^^. 

^ dx ' ^ dy ^ ^* 

These two equations in conjunction with (9) enable us to 
determine a, J, c as functions of x, y, z, p, q. Let these 
values be 

Then substituting in (9) we have 

equations which hold for aZ? such points. These are then the 
partial differential equations of the first order of the envelope. 

Now each of these equations is of the general form (4) ; 
whence by Prop. i. the partial differential equation of the 
second order is of the form (3), as was to be proved. 
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Let US actually construct this equation. 

Differentiating the first of the equations (10) with respect 
to X and to y, and regarding therein a as a function of those 
variables, and b and c as functions of a, we have 

d'f . MV . ^y ^* dY dc\ da 
etc* \dadx dbdx da dcdx da) dx ' 

^^ ^y ^(dj , dy db ^ dydc\da 
dxdy \dadx dbdx da dcdx da) dy ' 

from which we readily derive 

\ da? ) dy \ dxdy) dx 

Proceeding in the same way with the second equation of 
the system (10) we have 

(^_d^\da_f_dj\da^^^ 
\ dxdy) dy \ d^ ) dx ' 

Hence, eliminating -v- and -7- , we have 

('-^)'-('-S)('-f)"<- 

d]^ dxdy dx^ dj? dy^ \dxdy) ' 

the equation sought. 

Comparing this with the general form (3) we have the equa- 
tions 

dy* dxdy _ da^ _ 1 _ da? dy* \dxdy) 
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whence eliminatinff v4> -t«j ^^^ , ^ we arrive at the 

° dor ay axdy 

eqaatiou, 

£f^+4(trF-iZr)=0. 

This then is the condition which must be satisfied in order 
that the equation (3) may admit of an integral representing 
the envelope of a system of surfaces in which three parameters 
vary in subjection to two connecting conditions. It is only 
proved however to be a necessary^ not to be a sufficient^ con- 
dition. 



Solution of the equation Br + Ss ■{- Tt -{■ U{s^ — rt) = V, when 
a first integral of the form F{uy v) = 0, eansts. 

4. In the following sections we propose 

1st. To shew that when a first integral of the above form 
exists, its discovery depends upon the solution of two simul- 
taneous partial differential equations of the first order re- 
solvable into linear equations. 

2ndly. To shew' how from such first integral or integrals! 
the second integral is to be obtained. 

Prop. ii. If the equation 

Itr+8s+Tt+U{s'''rt)^V 

admit of a first integral of the form F{u, v) =0, in which u 
and V are functions of x^ y, «,^, y, then will F{u, v) considered 
as a function ofx, y, z, p, q, and represented as such for brevity 
by F satisfy the two partial differential equations of the first 
order^ 



\dx) dq \dy I dp \dx)\dy) 

dp dq ' 
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\dq/ dq dp Kdp ) 



}\dx) dp \dy J dq] ' 



in which 

dF . dF fdF\ dF . dF 

dz ' 



(dF\dF dF (^\^dF 
\dxJ^dx -^ dz^ \dy ) dy ^ 



Eegarding the function F in the proposed integral F= 
simply as a function of ic, y, «, jp, q^ we have 



'dF\ dF dF 
dx J dp dq 



+ -7— r + -J- « = 



^rfF\ di^ dF ^ ^ 



h% 



y (11). 



dy J dp dq 



On the other hand, regarding F as a function of x, y, z, p, q, 
mediately through u and v, we have the system 



du \\dx) dp dq J dv ]\dx) dp dq ) ^ 



dF {fdu 



H12), 



\ldu\ du ,du\ dF (fdv\ dv ^^ 1 ^ 
du\\dyj'^^''^ dq^rdi;\[d^)'^'^'^^^^ 

and these systems are equivalent 

dF 
Now if from the second of these systems we eliminate -j- 

dF 
and -J— , we obtain (Art. 2), a result which must be equiva- 
lent to the proposed partial differential equation, 

Br + 88+ Tt+ U{s' ^rt) ^r (13). 

This equation then considered as a relation between r, 5, t, 
must be an algebraical consequence of the relations (12), and 
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therefore of the equations (11). If then we determine alge- 
braically two of the quantities r, s, t, (we select r, t) from the 
system, and substitute their values in (13), that equation 
ought to be satisfied independently of the value of the re- 
maining quantity s. Now supposing^ and q to be both con- 
tained in jP, so that neither -7- nor -v- vanish, we have 
from (ll)y 






\dx) dq 



(dF\^dF^ 



dF 
dp 



\dyJ dp 
dq 



substituting which in (13) there results 



\dx J dq \dy J dp \dx J \dy J 

I \dq/ dq dp \dp J 



dF dF 
dp dq 



dq dp ' ^ \dp 

^-(S)f-^(f)f|---<»)- 



Now as this equation is to be satisfied in virtue of the con- 
stitution of jB, 8, T, U, F, and the function F, and indepen- 
dently of «, both the coefiicient of s and the absolute term not 
containing s must be separately equated to 0. Thus F con- 
sidered as a function of x, y, z, p, q, and containing p, q, at 
least must satisfy the partial differential equations 



^fdF\dF^^fdF\dF 
\dx) dq \dy ) dp 



-\ 



+ 



^©(f)^"^^"- 



dp dq 



R i^jy .8^' ^ + T (^S\' 



\diJ 



dF dF 
dp dq 



V 



\dpJ 
\\die J dp \dy J dq ) 



(15). 
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This result may also be established by forming the equa- 
tions of. condition which express the proportionality of 
i?, 8, ..,V, to the corresponding quantities m the constructed 
equation (7). From these equations of condition it is actually 
possible to eliminate in two distinct ways the quantities 

l^J' l^)' TF^ TT' *^^ result being the formation of two* 

partial differential equations for u agreeing in form with those 
above given for F, (See the memoir Ueher die partielle Diffe- 
rentialgleichung,,. Crelle's Journal, Vol. 61.) The actual 
transition from the former to the latter rests upon the con- 
sideration that the equation F {u, v) = 0, when F is arbitrary, 
is not really less general than the form ^ {F {u, v), v] = 0, m 
which the <I> is arbitrary. And here u has been replaced by 
F{u,v). 

The only condition respecting the application of the above 
equations is that we do not admit any relations which make 

either -t- or -r- to vanish. 
dp dq 

5. Prop. hi. The solution of the system of partial diffe- 
rential equations established in the last preposition may in all 
cases he made to depend upon that of simultaneous linear partial 
differential equations of the first order. 

In demonstrating this proposition we shall consider first 
the case in which U= 0, then the case in which F= 0, lastly 
the case in which neither of these quantities vanishes. The 
ground of this division will appear in the investigation. 

Case 1. Suppose ?7= 0. The equation then is of Monge's 
form, 

Rr+8s+Tt^V. 
The second equation of the system (15) becomies 

\dq J dp dq \cip J ' 
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and therefore breaks up into the equations 
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dF dF „ 



dF dF „ 
d^-'^'dp-^'^' 



m^ and m^ "being the roots of the quadratic equation 

Rm^-8m'\-T={i (16). 

As each of the above constituent equations is of the form 

dF^ dF 
dq dp ^ 

the system (15) may be reduced to the form 

Jt (^\iE. T(—]~ + Vm — — =0 
\dxj dp \dy J dp dp dp 



which breaks up into the equations 



dF 
dp 



= 0, Bm 



\dx ) \dy ) 



TT dF ^ 
dp 



The former of these we must reject (Art. 4). There re- 
mains for the determination of F the system of linear partial 
differential equations 



dF dF ^■^ 

-J m-Y- =0 

aq dp 



Bmf^"^ 



\dx) 



Kf) 



dp 



(17), 



and there will exist either one or two systems included under 
this form, according as the roots of the quadratic (16) airife 
equal or unequal. 



B.D.E. !!• 
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Case II. Let F=0. The system (15) then becomes 
\dx) dq^ \dy) dp \dxj \dy ) ' 

\dq) dq dp \dp J 

\\dxjdp \dy) dq]" 

Eliminate Uhj multiplying the first equation by 

(dF\dF (dF\dF 
\dx) dp \dy J dq ' 

the second by 

'dF\ fdF\ 
4x)\dy)' 



(; 



and subtracting; we obtain, after rejection of the common 

^ ^ dF dF 
factor -7- -J- 7 
dp aq 

\dxj \dx/\dy) \dy ) " 

We shall put this equation in the place of the second equa- 
tion of the system. This we are permitted tp do under the 
restriction that in seeking to satisfy the system so changed 
we do not make use of any relations which would cause either 
of the two factors employed in the process of elimination to 
vanish or become infinite. 

The new etjuation reduces to one equation, or breaks up 
into two equations of the form 

. (S)-»(D -•:■•■■• (-). 

m being determined by the quadratic equation 
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Making (^t— ) = wi i-j-j in the first equation of the system 
(15), we get 

which breaks np into 

But if we combine the first of these with (18), we obtain 

(f)=»- ©='• 

and this combination causing both the factors employed in 
the elimination of U to vanish must be rejected. There 
remains then the combination 

and this will represent either one or two systems of equations 
according as the quadratic determining m has equal or un- 
equal roots. 

Case III. Let neither Cr= nor F= 0. 

Multit)ly the second equation of the system (15) by 
indeterminate quantity Z, and add to the first ; then we ni 



(19). 



an 
ave 



-'(f)'^-'(f)-^(S)( 



dF\ fdF\ 
dy) 



\dxj dp \dy ) dp 

\dy) dq \dx) dq 

*(''-^ff-« ■•■•■••••» 
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We shall enquire whether it is possible so to determine I 
as to resolve this into linear factors. 

We might investigate this by resolving the. equation as 
a quadratic with respect to -j- or -p-. . But the form of 
the equation suggests what the forms of the linear factors 
must he if the resolution be possible. For as the squares of -v- 

J rp 

and -7- both appear, and these squares alone, in the fiinc- 

tion to be resolved, it is clear that -j- and 1— will be the 

dq dp 

only differential coefficients of F which will appear in both 

linear factors in common. The most general supposition 

possible is then that one factor shall contain -7- and -r- with 

dq dp 

-7- J , the other the same with ( ;?- J . 
Assuming then one factor to be of the form 



I 



dF ^ dF /dF\ 
dq dp \dx J ' 



it is seen from the. form of the coefficients of the first three 
terms of (20) that the other factor must be of the form 



j^dF I^dF U/dF\ 
dq m dp n \dy)^ 

and the resolved form of (20) must be 

\ dq m dp n \dy J) \ dq dp \dxj' 

Multiplying out and equating coefficients, we obtain the 
conditions 



0. 
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m 

n 

n 

rpn 

m 

The third and fourth of these conditions are equivalent, 

and give n = l. The first and second are also equivalent, 

T 
and give w = -^. These values reduce the last equation of 

condition to 

so that I is determined hj a quadratic. The resolved form 
of equation (20) now becomes 



{-f--f--(D}{'f4f- 




= 0. 



To these results we may give a somewhat simpler form by 
making Ul=m; not the m used above. We have then as 
the quadratic for determining «i, 

m'+Sm + RT-- UV= (21), 

and as the resolved form of (20), 
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Let TWj and m^ be the values of m. Then we have from the 
last the two distinct equations 

dF 



1 



dq ' dp \dy 



))f.f--f--©}-. 



^dF dF r^,^ 



dq ^dj> 



(f)} f ^ 



<^\T^^-^um\^o, 



dq dp 



\dx J 



and it is evident that these will be together equivalent to the 
equations (15) from which they were derived. 

Now to satisfy these equations simultaneously it is neces- 
sary that we should equate to one linear factor from each of 
their first members. If we equate to the first linear factors, 
we have 

n^F dF (dF\_ 



dF dF 



dq *dp ' ^ \dy 



^KI)' 



0; 



whence, by subtraction, 



/ sdF ^ 

This combination must therefore be rejected (Art. 4). For 
the same reason must the combination formed by equating to 
the second linear factors in the left-hand members of the 
above two equations be rejected. There remains then only 
the combinations formed by equating to the first factor of 
one of these members, and the second of the other. 

Thus we should have the combination 



^dF^ dF^jj(dF\ ^ 



dp \dy ) 



dF 



dF\ 



V 



(22), 



dq dp 



^(i)- 
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with the combination which would be obtained from tliis by 
interchanging m^ and m^. 

• 

6. It results from tlie foregoing investigations that the 
ftmction F is in all cases to be determined by the solution of 
two simultaneous linear partial differential equations with 
Jive independent variables. Now the theory developed in 
Chapter xxv. shews that the number of integrals of such a 
system cannot exceed three. That theory enables us both to 
determine what the number, of integrals is, and to construct 
the system of ordinary differential equations, reducible to the 
exact form, upon which their discovery depends. 

We iiave seen that the knowledge of two integrals u = a, 
V = J of the system enables us to construct a general first 
integral 

of the partial differential equation (3). And the solution of 
this first integral would lead us to the second integral which 
is the final object sought. But the direct solving of a partial 
differential equation of the first order which is not linear and 
which involves in its actual expression an arbitrary function 
is diflScult, and happilv it may be avoided here. The fol- 
lowing propositions will enable us to accomplish the virtual 
solution by a different solution, founded however upon the 
same general principles. 

7. Prop. iv. The integrals of the respective systems of 
simultaneous linear partial differential equ^ations upon which 
the determination of F depends are so related that y^ from two 
such respective integrals the values of p and q are determined^ 
they witl render the equation 

dz =pdx + qdy 

integrable. And in the particular case in which the two systems 
become identical^ any two integrals of the system stand in the 
same relation. 



136 PARTIAL DIFFEBENTIAL EQUATIONS [CH. XXVIII. 

For, let ^ be an integral of the system (22), and 'S? an 
integral of the associated system obtained by interchanging 
m. and m^ in the case in which these quantities are different. 
Then ^ satisfies the equations 






and "9 satisfies the equations 

But the necessary and sufficient condition in order that the 
values of j? and q derived from the equations 4> = 0, "^^ = 0, 
may render dz —pdx — qdy integrable, is 



© 



d^\ d9_d^ /dT\ 
dp dp \dxj 



Qf-?Q=» <->■ 



See Chap. Xiv. Art. 11, Equation (36). 

Now if from the previous equations we determine the values 
of 



(S). (f ). '^ '^ 



\dxj' \dyj' 

and substitute them in the above equation of condition it will 
be identically satisfied. 
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The determination of (-3-) > ••• fro°i the previous systems 

requires that U should not vanish. Hence the proposition is 
established except in the case of Z7=0, which is left doubtful. 

To examine this case let us revert to the system (17) which 
is proper to it. To that system since ^ 

5w"-/Sm+r=0, 

whence Bm^m^ = T^ 

we may give the form 



dF dF ^ (dF\ . fdF\ V dF ^ 
dg[ ^dp * \dxj '\dy J JR dp ^ 

or the form obtained from this by interchanging m^ and m,. 

Substituting in these respective forms ^ and '9^ in succes- 
sion for F, we find 

d^ d^ fd^\ fd^\ V d^ 



dg[ 

d9_ d9 f^_ fd9\ V d^ 

dq "^^ dp ' [dxj " "^^ \dy) R dp ' 



and these values substituted in (23) reduce it to an identity. 
Thus the proposition is established geperally. 

Lastly, as in the case in which the two roots of the quad- 
ratic for determining m are equal, the two systems of partial 
differential equations for determining ^ and "9 become one, it 
follows that if from two integrals of that one system we can 
deduce values of ^ and ; these values will render the equation 

dz^pdx-^qdy^^O 
integrable. 
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8. Prop. V. When the system of simultaneous linear par- 
tial differential equations determining F admits of two integrals 
u = aj v = J, it will admit or will not admit of a third inte- 
gral w = Cy according as the roots of the quadratic determining 
m are equal or unequal. 

The system in question, (22), becomes when we divide by 
Z7and write for (-i— ) and ( ^j- ) their full expressions 

dF dF m^dF ^dF_ 
dy'^^ dz'^ U dp'^ U dq" ' 

dF dF TdF rn^dF_ 
dx^^ dz^ U dp'^ U dqT ' 



or \F=0, A,J'=0, 



in which 



A _^ d d m^ d R d 

. _ rf d T d m^ d 

Hence the equation 

(AA-A,A,)i?'=0 • 



becomes 



m, — m, dF 

-^ * + 



U dz 






In this expression the coefficient of the first term only has 
been calculated. 

Now, by the theory developed in Chap. xxv. in order that 
the two simultaneous partial differential equations should 
have their full complement of integrals (three) it is necessary 
that the above equation should be satisfied identically. This 
involves three conditions, namely, 
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7^^-7^2 = 0, 

the first of which is the one aflSrmed in the Proposition to be 
necessary. 

Secondly, it is to be shewn that if this condition be satis- 
fied and if the system of given linear equations admit of two 
integrals w — a, t? = J, it will admit of a third. 

Replacing w^ and m^ by m the system becomes 

dF dF m dF E dF_ 
dy'^^ dz'^ U dp^ U dq"^' 

dF dF TdFmdF^ 
dx ^ dz U dp U dq 

Now if we construct from this the corresponding system of 
ordinary difierential equations, we shall find it to be 

dz — j?daj — qdy = 0, 
dp- yjdx—jjdy^^ 

dq — 7/^ ■" Yt ^y^ 

Now it is impossible that the first of these equations should 
be integrated without a previous determination of p and q as 
functions of a?, y, «, seeing that dr, dy, dz are the three differ- 
entials entering into that equation. Such determination can 
only come firom the integration of the second and third equa- 
tions of the system. But if these equations can be integrated 
in the forms w = a, v = J, then u and t? being particular values 
of F satisfying the partial differential equations, it follows 
from the last Proposition that the values of p and q which 
they will yield will make the first equation integrable. 
Hence if the system admits of two integials it will admit of 
three ; as was to be shewn. On the basis of these Proposi- 
tions the theory of the second integration rests. 
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TTieory of the Second Integration. 

* 

9. First suppose the values of m unequal. 

Then Wj = Oi, Vi = Ji being the two integrals (and we have 
seen that there cannot be more than two) of one of the systems 
of linear partial differential equations, and w^ = a,, v, = i, those 
of the other, the general first integrals of the given system 
will be 

4>(wj, t?i) = 0, ^(t^„ t?,)=0. 

The values of p and j determined from these will by 
Proposition iv. render 

dz '-pdx — qdy = 

integrable, and the integral of this will be the general integral 
of the proposed partial differential equation. For it will in- 
volve explicitly or implicitly two arbitrary functions derived 
from those in the first integrals. 

It suflSces however, following herein Charpit's method, to 
combine one general first integral derived from the one sys- 
tem with a particular first integral derived from the other 
system, e.g. the integrals 

4>(Wj, vj = 0, Wj = a. 

The values of p and q hence derived, and employed as 
before, will lead to a second integral involving one arbitrary 
function and containing two arbitrary constants. This con- 
stitutes a complete primitive from which the general solution 
will be obtained by converting one of the arbitrary constants 
into an arbitrary function of the other, and eliminating the 
latter between the equation and the one derived from it by 
differentiation with respect to that constant. 

Secondly, suppose the values of m equal. 

In this case we have but one system of partial differential 
equations so constituted however that if it admits of two inte- 
grals it will admit of three. 
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Let u=a, v = J, w=^c represent these integrals. Then if 
from these we eliminate j? and q we shall obtain a final inte- 
gral of the form 

.and this constitutes a complete primitive from which we shall 
deduce the general integral by making J = 0(a), c ='^(a), 
and eliminating a between the equations 

Q ^ df{x, y, a, <l> (g), ^fr (a)} ^ 

da 

To prove this let us combine the general and particular 
first integrals 

The values of p and q hence obtained make 

dz—pdx^qdy=^0 

integrable, and the result can be no other than the remaining 
integral «? = c, or rather what this would become on eliminat- 
ing^ and q from it. But since the equations by which this 
integration are to be effected are equivalent to 

w will become a function of a?, y, z, a and (f> (a). Also by 
Charpit's method c is to be treated as a function of a, so that 
ultimately we have the result above assigned. 

We have here supposed U not to vanish. If it do the 
theory assumes another but simpler form. Let 

be the two general first integrals. Then, since by the con- 
dition at the close of Art. 2, if ^ be eliminated from these 
equations q will also disappear, it suffices to eliminate them 
together in order to obtain tne general second integral. 
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10. Although the cases in which ?7= and F= have in 
the foregoing sections been treated for simplicity apart, their 
theory might have been deduced from that of the case in 
which neither Z7nor F vanishes. 

Thus to deduce the equations for the case of Z7= elimi- 
nate from the general system (22) -^ and -^ in succes- 
sion, and we find 

But from (21) RT- m^m^ = UV. 

Substituting, and then dividing by ?7we find 

the equation determining m^, m^ being 

This is equivalent to the results of Art. 5, Case i. 

11. We found it necessary (Art. 3) in order that the gene- 
ral partial dififerential equation of this Chapter should be satis- 
fied by tlie envelope of a system of surfaces the equations of 
which contain three parameters varying under two conditions 
that the relation 

fif« + 4(Cr7-i?T) = 

should be satisfied. 

It appears from Art. 8 that this is but one of three condi- 
tions necessary and together suflScient for this purpose. The 
formal conditions for eveiy form of ultimate solution con- 
sistent with the existence 01 a general first integral -F(tt, t?)— 
can be deduced in the same way. 
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[In the Bulletin de VAcocdimte Imp^nale des Sciences de 
8t P^iersbourg, Vol. IV. 1862, there is an article entitled Con- 
siderations sur la recherche des integrales premiires des iqua- 
tions diffSrentieUes partielles du second ordre^ par G. JBoldt 
(Lu le 7 Juin 1861). 

The article occupies pages 198 — 215 of the volume. Al- 
though the name does not quite correspond, I consider that to 
be a misprint, and I attribute the article to Professor Boole, 
partly from the nature of the contents, and partly because it 
IS known by his friends that he was engaged at a time corre- 
sponding to the date here given in the preparation of a mathe- 
matical article in French. 

The object of the article is to determine the conditions 
necessary for the existence of a first integral of the equation 

j^d^,g_^^ T^^W^O 
dix? dxdy dy^ ' 

where jB, 8, jT, and TFare any functions of x, y, z, -y- and --r- ; 

and also to determine the conditions which must hold in order 
that Ampfere's method of integration may be employed. 

In Crelle's Journal, Vol. LXI. there is an article by Pro- 
fessor Boole, entitled tleber die partielle Differentialgleichung 
zweiter Ordnung Itr+ 8s+ Tt+ U{s*'^ rt) = V, 

The article is dated 1862 ; it occupies pages 309 — ^333 of the 
volume. 

Among Professor Boole's manuscripts I found a memoir 
very closely resembling the article in Crelle's Journal; it 



144 PARTIAL DIFFERENTIAL EQUATIONS &C. [CH. XXVIII. 

would appear that the memoir was drawn up with a view to 
publication in the Transactions of some English Scientific 
Society, and that this design was afterwards abandoned in 
favour of the article in Crelle's Journal. 

After some hesitation I have resolved to print this memoir. 
Even if the memoir had been identical with the article in 
Crelle's Journal it would have been convenient to the English 
reader to be able to avail himself of the investigations ; and 
the memoir contains remarks which do not occur in the article, 
and which are interesting in connexion with the history of the 
subject. There is some repetition of matter which has already 
been given in Chapter xxvili. ; but I was unwilling to impair 
the completeness of the memoir by abridgment or omission. 
Accordingly the memoir forms the next Chapter of the present 
volume. 

In Article 2 of the next Chapter will be found the pro- 
cess to which there is an allusion towards the end of Article 4 
of Chapter xxvill. 

It is obvious that the subject of partial differential equa- 
tions of the second order was mucn studied by Professor 
Boole. The chronological order of his writings on the sub- 
ject appears to be as follows : 

1. Chapter xv. of the first edition of his work. 

2. The article in the Bulletin of St Petersburg. 

3. The memoir which forms Chapter xxix. .of the pre- 
sent volume. 

4. The article in Crelle's Journal. 

5. The Chapter xxviii. of the present volume.] 
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CHAPTER XXIX. 

ON THE SOLUTION OF THE PARTIAL DIFFERENTIAL EQUATION 

Iir+8s+Tt+ Uis^'-rt) = F, IN WHICH B, 8, T, U, V 

ARE GIVEN FUNCTIONS OF X, y, Z, p, J. 

1, The equation, the theory of the solution of which 
I propose to consider in this paper, is remarkable from its 
connexion with Geometry. If the equation of a surface 
contain three constants which vary as parameters in sub- 
jection to any two conditions connecting them, the gene- 
rated envelope will satisfy a partial differential equation of 
the above form. In other words any envelope of the surface 

F{x,y,z, a, &, c) = 

formed by the variation of a,h,cm subjection to two con- 
necting conditions 

<^j (a, J, c) = 0, <^,(a, J, c)=0 

is necessarily an integral of a partial differential equation of 
the form given above. 

Now this theorem is the more important, because it is 
only when three parameters in the equation of a surface 
vary in subjection to two relations that the envelope pos- 
sesses, irrespectively of the form of the connecfing relations, 
any definite character. If there be but one connecting rela- 
tion it is possible to determine that relation so as to make 
the envelope assume the form of any surface whatever, and 
therefore the possible system of envelopes is in such case 

B.D.E. II. 10 
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unlimited. If there be three connecting relations the para- 
meters become absolutely constant and no envelope exists. 

The partial differential equation 

is remarkable also as including all the cases in which a 
partial differential equation of the second order admits a 
first integral of the form 

w=/(i?), 

u and V being definite functions of a?, y, z^p^ j, and/(f) 
arbitrary in form. 

Neither of these statements is suflSciently general to con- 
stitute a theory of the genesis of the partial differential equa- 
tion under consideration, but the second one is more general 
than the first, and is indeed sufiiciently so to serve as the 
grotmd of an investigation which connects the solution of 
the equation in all cases with the satisfaction of a system 
of simultaneous ordinary differential equations of the first 
order and degree. And this is the ground upon which the 
method of the paper will rest. I propose to shew, Ist that 
the solution of the given equation on the assumption that 
a first integral of the form u=f{v) exists requires the satis- 
faction of a system of two partial differential equations of 
the first order and second degree ; 2ndly that this system* may 
be resolved into four systems, each consisting of two partial 
differential equations of the first order and first degree, two 
of which systems 'are irrelevant and the other two relevant ; 
3rdly that the solution of the two relevant systems ulti- 
mately depends on the solution of a system of ordinary 
differential equations of the first order, and that from these 
ordinary differential equations the given equation of the 
second order may be deduced independently of the assump- 
tion abeve mentioned. 1 shall also discuss the theory of the 
second integration. And I shall exemplify another method 
of solution connected by a remarkable law of reciprocity with 
the above method. 



\%. ^ 
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First Investigation. 

2. Prop. I. lfu=f{v) he a first integral of the equation 

Br + 8s+Tt+U{s''-rt)=V. (1), 

then will u and v, considered as functions of x, y, z, p, j, each 
satisfy two partial differential equations of the form 



^du\^ 



fdu\* 



^ "IKS +4' (!)}-» 

P /du\ du ^/du\ du 
\dxj dq \dy) dp 

jjfdu\ (du\ y.du du ^ 
\dx) \dy) dp dq 



m, 



in 



or 



du . du 



du . du 



dx'^Pr.'''"'^^+^-dz 



In which f-^j and (-7-] stand fc 

respectively. 

To demonstrate this proposition we shall form directly the 
partial differential equation of the second order of which 
u^f{v) is an integral and, comparing that equation with (1), 
deduce the conditions for the determination of u and v. 

Differentiating u=f{v), first with respect to x and secondly 
with respect to y, we have 

du du dz du dp du dq 
dx dz dx dp dx dq dx 



— -f't \ i— JL^^ dz dv^dp dv dq] 
""•^ ^ [dx dz dx dp dx dq dx) ' 



du du dz du dp du dq 
dy dz dy dp dy dq dy 



_ ^, . . (dv dv dz dv dp dv dg 
[dy dz dy dp dy dq dy) * 



10—2 
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.- ^ du , du y (dv\ du du , fdu\ 

or, rf we represent s+i. 5 V (^J . ^+5^ l? (^j . 

S ''5' ••• I "■* S ''J' '• ""* I '^ '• 

(|)+'|+'|-/»{(|)+'|h-'|}- 

Eliminating/' [v) we arrive at the partial differential equa- 
tion of the second order, 

(du (dd\ _^ dt? /^M 
Ic^ \dy) dp \dy)) 

(du (do\ _ dv /du\ dv /du\ _du fdv\) 
\dq \dy) " dq \dy) dp \dx) dp \dxJ) 

(dv (du\ du (dv\\ (du^ dv du dv\ . ^ . 

= (I) ©-©(!) (»)• 

It is seen that as respects the mode in which the quan- 
tities r, s, t are involved this equation is of the same form 
as the given equation (1). That it may be equivalent, its 
coefficients must stand to those of (1) in a common ratio fi. 
This gives 

du (dv\ dv fdu\ t, , v 

^l^J-^U^j='*^ («)' 

du /dv\ ^ dv fdjA dv^ fdu\ du (dv\ _ „ .^. 
dq\dy/ dq\dy) dp\dx)^dp\idx)~^ ^ '' 

dv fdu\ du (dv\ ^ , . 

^ te; "^ \^ ='*^ (*=)' 
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du dv ^u ^_ jj /7v 

dqdp dp dq ^ '' 

©©-(D©-"" («)■ 

As we have here five equations which are homogeneous with 
respect to the four differential coefficients of v and to fi, it is 
clear that we can, by the elimination of these quantities, 
obtain a relation connecting the differential coefficients of u 
with B, S, T, &c. But the peculiar cyclical form of the 
iunctions in the first members of the above system enables 
us to effect this elimination so as to lead to Jwo final equa- 
tions independent of v and fi. 

Thus .nltipl^^g (.) b, (I) |. (0) ,y (|) |, 

^^ ^y (^) (^) ' *^* ^*) ^y ^ ^' *°^ ^^^°S' ^^ ^''^' ''^ 
rejecting the common factor fi, 

n /du\ du rp(dv\ du —./du\ fdu\ 
\dxj dq \dy) dp \dx) \dy) 

+ r^^=o (4). 

dp dq ^ 

Again, multiplying (a) by (^) , (h) by (^) (^), (c) by 

y ' and (e) hj[^)^ + {^)^, addmg, aad agam re- 
jecting the common factor /a, we have 



+ 



-{©I*©!}-" (')• 
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Hence, u considered as a function of x, y, z, p, q satisfies 
the two partial differential equations (4), (5), botn which are 
of the first order and second degree. 



As u and v enter symmetrically into the system (a), (J), 
&c., V will also satisfy two partial differential equations of 
the same form, viz. the equations 

\ax/ dq \dy/ dp 



\dxj \dyj dp dq 

\\dxj dp \dyj dq) J 



> ... (6). 



Further, these two systems of equations constitute the 
complete system of equations resulting from the elimination of 
fi from the five equations (a), (J), (c), &c. ; for in their deter- 
mination, no factor involving either the differential coefficients 
of w and V, or the quantities -B, 8^ T, &c. has been rejected 
directly or indirectly. 

I am not aware that the above results of elimination have 
been noticed before. 

3. Prop. ii. The system of partial differential equations 
above obtained for the determination ofu, viz, 

J. (du\ du rp (dv\ du 
\dx) dq \dy) dp 

4. TT (^^\ (^^\ 4- F— — — 
\dx) \dy) dp dq~ 

^(l)'^«(S(|)-^(|)' 

-rrUduX du fdu\ du[ _ 
'^^^)d^'^\d^)Tq)^ 







" (7) 







ABT. 3.] OF THE SECOND 0£DE£c 151 

admits of resolution into four systems^ each consisting of two 
linear partial differential equations of the first order. Of these 
systems two only are relevant to the solution of the problem. 

For, multiplying the second by an indeterminate quantity 
X, and adding the result to the first, we have 



a. P (^^\ ^4. V— —M t(^^\ ^^ 
\dx) dq dp dq \dy) dp 

+'^{(s)l+($)l)-« «• 

Now let us see if it is possible to determine X so as to 
make the first member of the equation resolvable into linear 
factors. We cannot say (i priori that such resolution is pos- 
sible as we should be able to do if that member were homo- 
geneous and of the second degree with respect to three instead 
of with respect to the/owr subject variables 

ldu\ (du\ /du\ fdu\ 

[d^j' \dy)' \d^)' v^r 

A du « du J. 

Observing that the squares of -j- and ^ are wanting 

in the first member of (8) while those of f-r-j and i-jA 

appear, we are led to assume as the proposed equivalent of 
that member an expression of the form 

Multiplying the factors of this expression together and 
then equating the coeflSicients with those of the first member 
of (8) we have 
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Bm+\m^U+8\ (a), 

miti^TX (J), 

Xn = XF=7?in' (c), 

Hn =R {d)^ 

nm'^T (e), 

nn'^V (/), 

From (J), (c), (e?), we find 

T 
n= F, 7i' = l, w = XF, m' = -p., 

values which will be found to satisfy {e) and (/) also, and 
which reduce (a) to the form 

Supposing \ thus determined, the equation (8) becomes 

The. result is a little simplified if we retain m in place of X. 
We thus find as the resolved form of the given equation 

{^©-»(|)+''|}{Kl)^^(l)^''|}=''-'"' 

m being determined by the quadratic 

w«-/Sm + jBr-D'F=0. 

If Wj, wij be the values of m thus found, we have 

{*(l)^»'.(|)-''|lh©-<|)-''|}=»- 
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and these two equations are manifestly together equal to the 
system (7). 

Now these equations can only be simultaneously satisfied 

by equating to 0, one factor in the first member of each ; and 

the different combinations which are thus possible give rise 

to four binary systems of linear equations. Let us examine 

these systems separately. 

« 

If we simultaneously equate to the two first factors of 
the left-hand members of the last two equations, we have the 
systems 

a system which, when m^ and m^ are different, is reducible to 
the system 



\dx/ dp ^ ' \dy/ "" 



dp ' \dy/ 

It is clear that this cannot lead to a value of u satisfying 
the given differential equation (1), because it takes no 
account of the forms of 8, ?7, and T. Indeed if we actually 
eliminate 

du\ fdu\ du du 



(dAj\ fdv\ du du 



<dy/ ^ dp^ dq 
from the above equations by means of the system 



/dv\ du ^'^ ^fS 
\dx/ dp dq ~ 

/du\ du du 



(10), 
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(derived from the assumed first integral u =f(v) by making 
f^v) = c, and differentiating the result first with respect to x, 
then with respect to y), we find as the result 

Again, if we equate to the two last factors of the right- 
hand members of (10), we have 



©^^©^"S-". 



dy/ dg 
which, if Wj and m^ are difierent, reduce to 

^^)^''l-». ©-»■ 

And it is evident that neither are these equations consistent 
with the given equation (1), because they take no account of 
yS, Z7) and -B. The equation of the second degree to which 
they actually lead is 

Fr+r(s'-rO = 0. 

There remain then the two systems formed by combining 
the first factor of each one of the first members with the 
second factor of the other, viz. 



Ml) 



ld\j(\ ^du ^ 



•^©H-^d)^"!-" 



> 



(11), 



^(1) 



(; 



-"•'I) -"I-" 



* \dxj \dj/) dq "" 



.(12). 
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That these systems are relevant to the solution of the pro- 
blem under consideration may be shewn by eliminating from 
either of them by means of (10) the quantities 

/diA fdu\ du du 
[di)' \dy)' 'dp' dq' 

The actual result will be 

F{i?r + &+3T?+D^(«»-rO-F} = (13), 

which, except in the particular case of V— 0, reduces to the 
given equation. 

More generally, if in the equation 

u =/(v) 

u and V are any distinct solutions of the system (11), the 
same result of elimination may be deduced. For v by hypo- 
thesis satisfies the equations 






Subtract these equations multiplied by/'(t;) from the corre- 
sponding equations of (11), and representing m— /(v) by TT, 
we have 

which being of the same form as (11) it follows that 

Tr=0 or w-/(v) = 

, also leads to the partial differential equation of the second 
order (13). 



m 
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4. Prop. hi. To reduce the determination of the first inte- 
grals of (1) to the solution of a system of ordinary differential 
equations. 

Each of the systems (11), (12) presents u as satisfying 
simultaneously two linear partial diflferential equations of the 
first order. 

To deduce the value of u thus conditioned it will obviously 
suffice to multiply in each system one of the partial differen- 
tial equations by an indeterminate multiplier \, to add the 
result to the other equation so as to form a new equation 
which will, like those firom which it is formed, be linear and 
of the first order, and which on account of the indeterminate 
character of \ will be equivalent to the two. From the 
auxiliary equations which we obtain in the process of solu- 
tion, \ must be eliminated. 

If in this way we combine the equations of the system (11), 
we have, on arranging the resulting equation according to the 
differential coefficients of w, 

(5+).„,)|+K+xT)| 

dp dq 

Hence we have the auxiliary equations 

dx ^ dy dp dq dz 

Ii + \m^~m, + XT^'V''\V''Iip + m,q + 'K{Tq + m^p) ' 

du = Of 

and it is to be remembered that 7?i„ m^ are the roots of the 
equation 

m^-Sm + RT'-TJV^^O. 
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Eliminating X from the first four of the above equations we 
have 

Udq + m^dx — Rdy = 0] 

Udp + m^dt/ - Tdx =^ oi (I). 

dz — j[hIx — qdy = o) 

This then is the system of ordinary differential equations 
deduced from (11) upon the integration of which the determi- 
nation of u will depend. 

A similar system, differing from the above only in the 
mutual transposition of mj and wj,, is given by (12), viz. 

Udq + m^dx — Bdy = 0] 

Udp + m,dy-Tdx = 0> (II). 

dz — pdx — qdy = O' 

1£ from either of these systems we can deduce two inte- 
grals of the forms 

t^ = a, V = J, 

it is obvious, from what precedes, that 

u =/(v) 

will constitute a first integral of the proposed (1), and there 
being two systems in question, two such first integrals, each 
involving an arbitrary constant may coexist. 

6. Pbop. IV. To deduce the second integral of (1), 

It will be necessary to consider separately the cases in 
which 7^1 and m^ are equal and unequal. 

First let m^ and m^ be equal. 

Both the systems (I), (II) reduce to a single system which 
may be expressed in the form 



dp = yx^aj — Yj^y 



m 



dz = pdx + qdy ^ 



y (14). 
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Now, since the condition ;/ =^ is l^cre satisfied, it is 

manifest that if from any two integrals of the above system 
of the fonns w = a, v = J, simultaneous values of p and a 'be 
detennined, these values will render the third equation ol the 
system integrable, and the effect of its integration will be 
virtually to determine « as a function of a?, y, and three arbi- 
trary constants, viz. a, J, and a constant c introduced in the 
last integration. Let us represent the result in the form 

« = ^(a?, y, a, J, c) (15). 

Now what relation will this result bear to the general solu- 
tion of the partial differential equation given, to the solution 
which we should obtain by integrating, not the particular 
equations w = a, v = J, but the general first integral u =/(v), 
which includes them both. 




dz —pdx — qdy = 



integrable, and to integrate the last equation regarding the 
arbitrary constant of integration as an arbitrary function of 
b. The result is a complete primitive in wnich, by the 
variation of 6 as a parameter the general integral is implicitly 
involved. 

Now either of the equations u = a, v = J will, in conjunc- 
tion with u =f{v) determine^ and q so as to make 

dz —pdx — qdy = 

integrable. Take the equation v = J, then u =f{v) reduces to 

u =/(5). 

Thus, in place of the equations w = a, v = J, of the previous 
section, we have 
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for the detennination of p and j. The constant c Introduced 
in the final integration becomes also, according to the above 
theory, a function of J, and the complete primitive is of the 
form 

^ = ^{a;,y,J,/(J),t(J)} (16), 

while the general integral is found by eliminating h between 
this equation and its differential with respect to J. 

The general integral therefore represents the envelope of 
the surface represented by (15), a, J, c being parameters sub- 
ject to any two connecting conditions. 

As m^ , m, are supposed equal, a necessary condition of the 
possibility of this species of integration is that 

>S»-4(5r-?7F)=0 (17), 

the value of tw is — , and the system (14) reduces to 



(18). 



dz — jpdx — qdy = - 



We conclude therefore that if (17) he satisfied and we can 
from (18) deduce a value of z in terms of a?, y, and three 
arbitrary constant^, the equation expressing that value will be 
a complete primitive, and the general integral will be found 
by making the constants vary in subjection to two arbitrary 
conditions. 

Ex. Let the given equation be 

xqr + ypt + a^ (5^ — rt) =^2. 

Here R = xq, 5=0, T^^yp, U=^xy, V^^pq. 
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The condition (17) is satisfied, and (18) becomes 

xydp —ypdx = 0, 

xydq — xqdy = 0, 

dz —pdx — qdy = 0. 

From tte two first of these we find 

p=^ax, q^hy, 

whence from the third, 

« = 1- - ■- + c. 

2 2 

This is the complete primitive, and the general primitive 
consists of all possible equations derived from this by makin 
a, J, c vary in subjection to two conditions. 

Ex. 2. Given 

Here the equation for m reduces to 

w? + 2pqm +jpV = ^> 
whence m = — ^g, and the system (18) gives 

(1 -\-p^ + 2 j' 



g 



dp 



(l+Z + S*)^ 



•\-pqdy + (1 4-^*) etc = 0. 
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Subtracting the upper equation multiplied by pg' from the 
lower one multiplied by 1 + g*, and dividing by 1 +]f + q^, 
we have 



whence 



'" + V(l+/+2')""' 



In like manner, 



2 =1. 



Hence determining p and q, 

, _ {x-a)dx-)r[y-h)dy • 

Therefore {x - aY + (y - J)' +[z- cf = 1. 

From this form of the complete primitive it is evident that 
the general integral will represent all possible tubular surfaces 
formed by the motion through space of a sphere of constant 
radius unity. 

Secondly, let m^ and m^ be unequal. 

Then since, in neither of the systems (I) and (II) is the 

condition -J- = -j satisfied, from neither system separately 

can values of _p and q be obtained which make dz=^pdx + qdy 
integrable. 

But, as will be shewn, any two integrals obtained, the 
one from the one system and the other firom the other, will 
give values of p and q which will render dz=:pdx + qdy in- 

B.D.E. II. 11 
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tegrable, and the general solution will consist of all possible 
integrals of the latter equation thus obtained. 

Or if the complete first integral of either system be com- 
bined with any particular integral involving an arbitrary 
constant obtained from the other, the two will fiimish values 
of jp and J which render dz=pdx + qdy integrable, and its 
integral will be a complete primitive involving one arbitrary 
function in its expressed form, another in the connexion of 
its two constants ; the general primitive being found in the 
usual way by making the constants vary as parameters in 
subjection to a single arbitrary connecting condition. 

In fact it may be shewn that if we attempt by the process 
of Charpit or Lagrange to integrate the partial differential 
equation of the first order w =/ (v), deduced we will here sup- 
pose from the system (i), we virtually construct the system 
(ii) in the auxiliary equations upon which the process of 
solution turns.. I have obtained a direct proof of this proposi- 
tion, but I think it preferable and at the same time suflScient, 
to direct attention to the prior ground upon which it rests in 
the relations of the systems of partial differential equations 
(11), (12) from which the systems of ordinary differential 
equations (i), (ii) are derived. 

Let P=0 represent any integral of the system (11), and 
Q = any integral of the system (12). Then we have 



dp " \dx) * \dy) 
dq^ ^\dxj \dy) 
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Hence we deduce 

_ y[dPfdQ\ _ dQ fdP\ dPfdQ\ _dQ fdP^\ 
Xdjp \dx) dp \dx) dg \dy) dq^ \dy)) 

={-©-. (f)}(i)-^©--(f)}(S 
M-(S-Kf)Kf)-{-©^^(f)}(f)- 



The second meml)er of this equation is identically 0. Hence 
dividing by F we have 

^f^^^^f^W^/^^^^^f^Vo fl9). 

dp\dx) dp\dx) dq\dt/) dgi\dy) ^ '' 

But this is the known condition under which the values 
of p and q deduced from the equations P = 0, ^ = make 
J2j=^cfe + g'c?y integrable; see Cfhap. xiv. Art. 13, Equation 
(36). 



We conclude then that if from the systems (i), (ll) we can 
dediLce two corresponding systems of integrals 

then will the first integrals of (1) le 

«'l=/lK)> «*2=/2W> 

while the second integral will consist of all possible relations oh- 
tained either 1st hy specifying the forms off , f and obtaining p 
and q a^ functions of x and y and integrating dz =pdx + qdy, 
or 2ndly, by specifying one of the functions f^ j^, leaving the 
other arbitrary, determining p, q, integrating dz = pdx + qdy^ 
and regarding the final constant of integration as an arbitrary 
parameter. 



^ — ady = 0> 

+ m^dx — ady = 0) 
+ m^dy — cdx = o) * 
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Ex. Given ar^-ls + ct-^-e {s^ - rt) = A, the coefficients 
being constant. 

HerejB=a, 8==h, T=c, U=e, V=h. 

Hence w,, m^ are the roots of 

w' — Jw + ac — cA = 0, 

and the systems (i), (ii) give 

edq + m^dx — ady = 0] 
edp + m^ 

edq + m^dx — ady = 0] 
edp 

Whence the first integrals are 

eq + WjOj — ay =/ (cp + m^y — ca:) , 

from which all possible second integrals are to be derived in 
the modes above explained. 

Let us take the second of those modes and give to the 
second of the above first integrals the particular form 

ep4- wi^y— ca?= (7, 

G being an arbitrary constant. From this, and from the 
other integral, left in its complete form, we have 

cx--m,y+C _ ay-- m^x+f, {(m^ - m,) y + 0} 
P" e '^'^ e ' 

whence, substituting in the formula dz—pdx + qdy^ inte- 
grating, replacing the arbitrary form 



jfi{t)dthy(m,^m,)4>(t), 
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and introducing an arbitrary function of C for the arbitrary 
constant, we have 



I Vex'' 



xy + ^+ Cx+(f> [{m^ '-'rn^)y+ C] 



+ t(C) 



for a complete primitive. The general primitive consists of 
all possible relations obtained by eliminating G between the 
above equation and 

when the forms of <f> and -^ are specified. 



Second Investigation. 

6. If from the equation 

Iir + 8s + Tt+U{s^^rt) = V. (20), 

we eliminate r and t by means of the equations 

dp = rdx + sdy, 
dg[ = sdx + tdy^ 
the result will be 

[Rdf - Sdxdy + Tda? - U {dpdx + dqdy)] a 

^Rdpdy+Tdqdx-Udpdq-^Vdxdy (21). 

There are different considerations (all of them however in- 
volving, as I have been led to think, a more or less explicit 
reference to some theory of the genesis of the given partial 
differential equation) which indicate that its solution depends 
upon that of the equations obtained by equating to the part 
affected and the part not affected by 5, viz. upon the solution 
of the equations 

Rdj^ - Sdxdy + Tdx^ - U{dpdx + dq^dy)^0 (22), 

EdpdyAr Tdxdx - Udpdq - Vdxdy ^ (23). 
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Without entering into these considerations let ns inquire 
what consequences may be deduced from these equations 
assumed to be true. 

It is seen that these equations are connected by a remark- 
able reciprocity with the partial differential equations (7). 
They will in fact be converted into these equations if we 
change 

dx, dy, dp, dq, U, V, 8 ' (24), 

into 

du du du du yr rr a /^rx 

-^'-^' d^' Tp' ^' ^'-^ (^^) 

respectively. From this formal connexion it follows that if 
we multiply (22) by \ and add to the result (23), we shall be 
able to determine \ so as to permit the resolution of the equa- 
tion thus formed into linear factors. Ultimately we shall, as 
appears from Art. 3, reduce the system (22), (23) to an equi- 
valent system of the form 

[-- Rdy - m^dx + TJdq) {-m^dy- Tdx+ Udp)=0, 

{'- Bdy - mjix + TJdq) {-m^dy-- Tdx+ R^) = 0, 

m, and m^ being determined by the equation 

or, changing the sign of m, 

(- Bdy + m^dx + Udq) {m^dy - Tdx + Udp) = 
(— Bdy + m^dx + Udq) {m^dy — Tdx + Udp) 

m^ and m^ being as in the former investigation roots of 



-^[...(26), 



Equating to the corresponding factors of the first mem- 
bers we have 



— Bdy + m^dx + Udq = 

— Bdy + m^dx + Udq 

m^dy--Tdx-\- Udp = 
m^dy — Tdx + Udp 



= 0) 
= 0)' 

= 0j 
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The first of these, m^ and m^ being different, is resolvable 
into 

Udq-Iidi/ = 0, dx = 0; 
the second into 

Udp-Tdx = 0, e?y = 0, 

and it is obvious that neither of these can lead to the given 
partial differential equation (1). The first of them combined 
with the equations 

dp=rdx + sdy, dq=^$dx + tdi/ (27), 

leads in fact to the partial differential equation 

B- Ut=^0 (28), 

the second in like manner leads to 

T-Ur^O (29). 

But equating to the non-wrresponding factors of the first 
members of (26) we have 



— Rdy + m^dx + ^dq^ = 
m^dy— Tdx+ Udp 

m^dy - Tdx^- Udp = 

— Rdy + m^dx + Udq 



::!• 

::i- 



Now these systems when completed by the equation 
dz = pdx + qdy agree with the systems (i), (ii) deduced in the 
previous investigation. 

It remains to shew that these systems actually lead to the 
given partial differential eq[uation (1) directly. Eliminating 
from either of them, combined with the system (27) the dif- 
ferentials dx, dy, dp, dq, we shall have as the result 

U{Er + Ss+Tt+U{s^-'rt)-r}=^0 (30), 

which, rejecting the factor U, as from (13) we rejected F, is 
the differential equation proposed. 
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Gr<mnd of the Reciprocity above noticed, 

7. The reciprocity above noticed is not of a primary cha- 
racter, but is founded upon two prior laws which I shall pro- 
ceed to demonstrate. 

If from the partial diflferential equations of the system (7) 
we eliminate Fand substitute the resulting equation in the 
place of the first equation of the system we shall obtain the 
equivalent system 

\a^J dq dp \dpl 



\\dx) dp \ 



du\ du\ ^ 



dy) dq^ 



...(31). 



J 



These equations are both symmetrical and it will T)e ob- 
served that they are convertible the one into the other by 
changing 



into 



du du du du jj Q j^ 
Jq' Tp' ^' dy' ^' ^' ^ 



du du du du y. q tt 

d^' Ty' ^' d^' ^' "^' ^ 



(32), 



(33), 



respectively. This is a law of reciprocity which connects 
solely the differential coefficients of u and the coefficients 
TJy 8, F of the original equation. 

Again w = is by hypothesis a solution of the given partial 
differential equation. Regarding it however simply as an 
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equation which is true and the truth of which is consistent 
with that of the equations 

dp=^rdx-\:8dy^ ' ^^^^^ 

dq = sdx + tdy) 

and differentiating it first with respect to a?, secondly with 
respect to y, we have 

du du dz du dp du dq __ 
dx dz dx dp dx dq dx" ^ 

du du dz du dp du da 
dy dz dy dp dy dq dy ' 

equations to which we may give the form 



du\ __ du du 
dx) dp dq 



\dy}~ 



du\ du . du 
dp aq 

Now this system is of the same form as the system (27) 
and will agree with it if we change 

du du du du , . 

■"^'■"^' ^' d^ ^"^^^^ 

into 

dp, dq, dx, dy (36), 

respectively — a change which does not affect the coeflScients 
of the given equation^ and which is therefore the expression 
of a law of reciprocity distinct from that last noted. The 
combination of these two laws does however lead to the 
law exemplified in the researches of the previous Article; 
see (24) and (25). 

The question here arises whether it would not have been 
better to employ from the first the symmetrical forms (31) of 
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the partial differential equations of the first order and second 
degree upon which u depends, than the unsymmetrical forms 
(7). It was indeed from the symmetrical forms that the 
chief results of this paper were originally obtained, but the 
unsymmetrical forms lead to the same end in a simpler way, 
and therefore they have been made use of in the present 
memoir. 

It may be proper to notice, in concluding this section, that 
the symmetrical forms in ordinary differentials would have 
emerged in place of the unsymmetrical ones of (22) and (23), 
if the quantity s^ — rt had been retained instead of «. The 
equations 

dp = rdx + sdt/j dq^ = sdx + tdy, 

enable us in fact to reduce the given equation (20) to the 
form 

Bdp^ + Sdpdq + Td4 - Vid'pdx + dqdy) 

= (s'-rO [Bdf^ Sdxdy-^^ Tda? - U{dpdx + dq^dy)]. 

Hence arises the symmetrical system 

Bdp^'^-Sdpdq + Tdg" - V {dpdx ■\- dq^dy) = 0, 

Rdr^-Sdxdy + Tda? - U {dpdx + dqdy) = 0, 

which is connected with the system (31) by the single law of 
reciprocity expressed in (35) and (36). 



Postscript. 

8. At the time when the above Investigations engaged my 
attention I was totally unaware that the subject of them had 
been discussed by Ampfere (Journal de VEcoU Polytechniqtiey 
Tom. XI.) and recently by Professor De Morgan {Cambridge 
Philosophical Transactions, Vol. ix. Pt. IV.). I feel it there- 
fore incumbent upon me to state why after acquainting my- 
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self with the results of their labours, I offer this paper for 
publication. 

The method of Professor de Morgan so far resembles the 
first method of this paper, and that of Ampfere the second, 
that while the former makes the solution of the problem 
depend directly upon that of simultaneous partial differential 
equations of the first order, the latter makes it to depend 
directly upon the solution of simultaneous ordinary differ- 
ential equations of the first order. The formal connexion of 
these methods by the law of reciprocity is, I believe, esta- 
blished for the first time in this paper. The system of partial 
differential equations of the second degree (7) has not, so far 
as I am aware, been given before. 

But a point which I think of deep importance is the follow- 
ing. By connecting, as in this paper, the differential equa- 
tions of the second degree, whether ordinary or partial, by an 
indefinite multiplier which is afterwards .determined so as to 
admit of the resolution of the system into its component linear 
elements, we assure ourselves that each step of the solution 
offers a complete sequence to that which has gone before, and. 
it only remains then to separate the different elements and 
determine whether they are relevant or irrelevant to the end 
in view. That any such distinction exists has not, so far as 
I am aware, been noticed before. And it seems to me the 
more important that it should be noticed because the solution 
of partial differential equations in cases far more general than 
those above considered seems to depend upon the satisfaction 
of simultaneous differential equations of a degree higher than 
the first. I have in fact by an application of the Calculus of 
Variations arrived at the conclusion that the theory of the 
solution of all partial differential equations of the second 
order, whatever the number of variables may be, is very inti- 
mately connected with the satisfaction of a system of differ- 
ential equations of the type 

jP= representing the given partial differential equation, x 
and y any two of the independent variables, and r, «, t the 
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second diflferential coefficients of the dependent variable with 
respect to x and y. 

I may perhaps at some future day resume the subject, to- 
gether with an inquiry into the theory of the solution of the 
partial diflferential equation of this paper, when the conditions 
under which the auxiliary equations (i), (ii) are supposed to 
be integrable are not satisfied. 

9. Note. It may be desirable to establish directly the 
converse form of one of the results of Proposition iv. For 
this object we shall shew that the equation of the envelope of 

= ^ (a?, y, a, &, c) (1), 

where a, 5, c are connected by any two conditions of the 
forms 

-^ (a, J, c) = 0, X («» *> ^) = ^» 
will satisfy a partial diflferential equation of the form 

Er+Ss+Tt^-U{f-rt)^V (2), 

in which also 

Diflferentiating (1) we have 

__ d(l> d<j> da d(f> db d4> dc 

d^ d^' da d(f} db d<l> dc 
^ dy da dy db dy dc dy 

and by the nature of an envelope these reduce to 

p-% »-i w- 

Again diflferentiating these equations with respect to x and 
y, and writing for simplicity 

dadx ' dbdx ' dcdx ' 

dady ' dbdy ' dcdy ' 
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we have 

da? dx dx dx^ 

__ d^^ sda jydh ^dc 
dxdy dy dy dy^ 

dxdy dx dx dx 

t — <^*^ . A' ^ ,Ti>^ , ri'^' 
"^ «y dy dy 

Hence we find 
V dxdyj \ da?)\ df) 

\ dy dy dyl\ dx dx dxj 

_(A^ + Bf+Cp)(A'^ + B-f + C'^] 
\ dx dx dx) \ dy dy dy I 

= {AB-A'B)(pf-^f) 

^ ^ \dy dx dx dy) 

+ {Ba-B'C)(f^-f^\ 

\dy dx dx dyj 

' \dy dx dx dyJ 

Now since a, J, c are connected by two conditions, so that 
h and c are functions of x and y onlj as being functions of a, 
we have 

da db da dh ^ dh dc dh dc ^ 

dy dx dx dy^ ^ dy dx di dy ' 

dc da ^ d^_/x 
dy dx dx dy" * 
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Thus the above equation reduces to 

V dxdy) V dx'jy dfj~^' 



or 
d^ 



dxdy dsif dof dy* \dxdy/ "* ^ '* 



This equation is of the general form (2). Its coefficients 

-j-Y- , &c. are determinable as fimctions of a?, y, z, p, q when 

the form of the complete primitive (1) is given. For this 
purpose the complete primitive with the two derived equations 
(3) suffice. 

Again, comparing (4) with (2) we have as the conditions 
of their equivalence 

R 8 T ^TT F 

d^ dxdy da? da? dy^ \dxdy/ 

conditions which suppose jB, 8^ T^ U, V connected by the 
relation 

8'^4:{BT^UV) = 0. 



^^--^. 
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CHAPTER XXX. 



ADDITIONS TO CHAPTER XVII. 



[The present Chapter consists of additions to Chapter xvil. 
Art. 1 was intended to follow Chap. xvii. Art. 1.] 

1. The theory of the solution of linear differential equa- 
tions in a series flows very beautifully from their symbolical 
expression. It is usual in treating this subject to assume the 
form of the series, and deduce from the differential equation 
the law of its coefficients ; but the symbolical form of the dif- 
ferential equation determines in reality the form of the solu- 
tion as well as the law of derivation of its successive terms. 

Let us begin with the binomial equation 

f,{D)u^f,{D)e^u = 0. 
Operating on both sides with {/^(Z))}"^, we have 

in which *(^^=='f(i- 

Hence {1 - <^ (Z>) e^]u = {/, {D)Y^O. 

Now {.;^(-D)}"^0 will be determined by the solution of a 
linear differential equation with constant coefficients, and will 
be necessarily of the form 

in which -4,5, (7, ... are arbitraiy constants, and P, Q, 5, ... 
are functions of the independent variable. 
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We have then 

therefore u = {l-<f>{D) e'^p (AP +BQ+CB+...). 

Now let us represent (f> {D) e*"^ by p ; then 
^, = (1 -p)-^ {AP+BQ + CB + ...) 
= {i + p + p^ + p'+ ...) {AP+BQ+CB+ ...) 
= A{l + p + p^ + p'+...)P 
+ B{l+p + p' + p'+...)Q 
+ C{l+p + p' + p'+..,)B 
+ ... 

Represent the first line of the above expression by w,, 
then since 

p*^ = ^ (Z>) e*"^^ (Z>) e'^...m times 

= €'"^^<^(i? + wr)<^(Z> + »wr-r) 4>{D + r), 

we have 

t^, = J[ {P+e'-^^CZ^ + r) P+ 6^^<^(2> + 2r) ^(Z> + r) P 
+ €«'-^<^(Z> + 3r) <^(i> + 2r) «^(Z> + r) P+ ...}, 

in which it only remains to perform the operations indicated 
by <f>{D + r), by ^(i? + 2r)^(Z> + r), ... on the ftmction P. 

Let us in the first place suppose the symbolic fanction 
/q (D) to be of the form (2> — a) (2> — 6) . . . ; then 

{fo{B)}''0=-A^^ + B€'^^ + .... 

Here P= e*^. Hence substituting in the above expression 
for u, and observing that/(Z>) e^=f{n) e'*^, we find 

or, since e^ = a?, 
u^ = Ax'' {I + ^(a + r) a?'"+ ^(a + 2r) ^(a + r) aj*'+ ...} ; 
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and 

u^^Ax'' {1 + (f>{a + r) x"" + <f>{a + 2r) <f>{a + r) 0^ + ...} 
+ Bx'' {1 + (l>{b + r)of+4>{b + 2r) «^(S + r) 0?*"+ ...} 

the solution sought 
Consider now the general equation 

Here we have, representing fjji hy ^m(-^)> 

{l + ^,(2))6»+.... + ^.(i?)e«»}M = {/.(i?)rO; 

therefore 

«={l + 0.(2?)e'+ <}>^{D)e"Y{fom-'0. 

Here we have first to determine {yQ(Z>)}"^0, then to deter- 
mine the effect of the operation represented by 

{i+^,(2>)6»+ <A»(i?)^r 

upon this. 

Now {/o(-D)}*^0 is given by the solution of a linear diffe- 
rential equation with constant coefficients, and will therefore 
be of the form 

AP+BQ + CB+..., 

A, B^ (7, ... being arbitrary constants, and P, Q, B, ,., func- 
tions of 0, 

Again, since 

{i+^.(i))6»+ +MD)^r 

^ 1 

i + <l>^{Dy+ + MD)^^' 

it may be shewn by a process of actual symbolical division, 

B.D.E. II. YL 
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attending to the laws of combination of symbols, that the 
expression may be expanded in the form 

To determine the functions F^{D), F^{D), we may 

proceed as follows. From the equation 

{1 + ^.(Z))€»+ + ^.(Z))€-»r 

= F,{D)+F,{D)^ + F,{D)e»+... 
we have 

l^{l+<f>,{Dy+...+<l>^{D)e«»]{F^{D)+F,{D),»+F,(D)e»+...] 
= FSD)+{F,iD)+4,,iJ))F,iD-l)}^+ (1). 

Hence F^{D) = 1, 

F,{I)) + <f>,{D)F,{D- 1)^0; 
therefore F^{I>) = - ^. {D)F,{D-1), 

and so on. Hence FJB), F^{D), are determined in suc- 
cession. The general law is as follows: the coefficient of 
e"^ in the second member of (1), when m is greater than 1, is 

F„{D) + 4,,{D)F^,{D-l) + <f>,(D)F^{I) - 2) + ... (2), 

whence 

F^{D) = -<f>,iD)F^,{D-l)-4,,iD)F^{B-2)-.... 

"By this formula the successive values of F^{D) can be 
deduced from those of F^^ (D) , F^^^ (i>) , . * . . 

Combining the above results we obtain thus for u the ex- 
pression 

u = {l + F,{D)^ + F,{D)e'^+....}{AP+BQ + ...} 

^A{P-\rF,[D)e'P+F^{D)^'P+....} 

+ B{Q + F,{D)e'Q + F,{D)e^^Q+...,} 

~t" ...» 

Let us in applying this expression first suppose that the 
factors of f^ (D) are real and unequal, so that f^ {D) is of the 
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form (Z) — a) {D — h) (Z) — c).... Further, let us suppose that 
no two of the quantities a, J, c, .... differ hj an integer. 

Then {/,(i>)}-^0=^6«^ + jBe»^+...., 

« 

whence we may assume 

Thus the expression for u "becomes 

+ jB{€W + ir^(2)) €(^«)«+ir (D) e^^+ ...:} 

"T" • • •• J 

or, since F{I>)e^ = F{m)^, 

u =A{^ + F^{a + 1) 6'"*"» + i?;(a + 2) 6"^*'«+ ....} 
+ B{^ +F,{b + 1) e'"^" + i?; (6 + 2) e<*"'» +....} 

"T" • • • • 

Hence, replacing «• by a;, 

« = ^{x» + i^Ja + 1) a;'*' + i?;(a+ 2) af** + ....] 

"l ...» 

In (2) replace in like manner D by a + t and we have, 
putting t for m, 

+ ^2(a+ -^i-aC^ + *'- 2) + = 0, 

or, if jFJ(a + i) be represented by u^+i, 

w«+»+^i(a+f)Wa+i-i + ^2(^ + w«+i-2+ = ^- 

Put m for a + i, thus 

««m+ ^1 W ^m-i + <^8 W w^a+ = 0. 

This agrees with the law established in [there is no refer- 
ence in the manuscript, but the law intended appears to be 
that given in Chap. xvii. Art. 9.] 

YL— 'I. 
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Secondly, suppose that r of the factors of f^ {D) are equal 
and of the form D — a. 

Then [fo{I>)]'^0 contains a term of the form 

^(c,+c,e+c,e'+ +cr.,er-'). 

Hence the corresponding portion of u is of the form 

k + F,{D)^ + F,{D)^+....^^{c, + c,0+ + Cr_,e^-') 

=^L'^+e^-^'^%{D +,a + l) + €<^J^i^,(i) + a + 2)4-...l v.. .(3) 

where v stands for c^ + cj9 + c^6P +.... + c^_l^"^ 
Now i^<(Z) + a + i)t; 

= |j?;(a + l) + i^/(a + l)i> + i^;'(a + l)^^+....|t;, 

which on performing the differentiations becomes a polyno- 
mial of the form 

A, + A^e+ +A,^^e^\ 

We see thus that (3) will assume the form of a series of terms 
€**^, e^""^*^^, .... each multiplied by a polynomial of the (r — 1)**^ 
degree in 0. Or arrangmg the terms otherwise it will con- 
sist of a series of terms of the form 

JBo + B^0+ +5,_,^-^, 

in which B^, B^, „..; B^_^ are series involving e"^ e^"^^^^, 

e^""^^^, Or lastly, changing e^ to x, the portion of t« in 

question is of the form 

j?,+5,(ioga;) + +B^^ aog^^^ 

J?Q, 5j, .... 5_j being polynomials in each of which the lowest 
power of X is aj**, and the successive powers increase by unity. 

This establishes the assumption in [there is no reference 
in the manuscript; probably Chap. xvil. Art. 10 is to be 
supplied.] 
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Thirdly, let /, (D) contain r factors D--a^, D — a^^.^.D — a^ 
in which a^, a^, .,.ar differ from each other by integers, toge- 
ther with other factors. 

The portion of u corresponding to the factor B — a^ will be 

in which 

K {.D) *"•» = - 1^, (2)) ^F^, {D) €<-)* 

Thus F^(P) 6*^ consists of terms of the form 

>i>i{D)i<>F^{P)^'^, 

i being one of the numbers 1, 2, .... w. Hence i']^ (2)) e*""^* 
will consist of terms of the form 

4>,{D)^F^{D)4'-^, 

j being one of the numbers 1, 2, ... w. Continuing this until 
* +y + ifc+ ... = m, we see that F^(P) e"*^ will ultimately con- 
sist of terms of the form 

^,(i))e%(Z))€^<^*(i>) €*«..., 

e,y, fc, ... receiving arbitrarily any of the values 1, 2, ... li, 
and i+y+i+ ... being equal to m. 

Thus the portion of w derived from A^^ will consist of 
all possible terms of the form 

Afi>m e^<l>j{D) ^<I>,{1)) 6*^ .... (6«0 
= ^<^,(Z)) <^,(Z) - i) ^(Z) - ^ -^l .... 6(-^*^^ 

- f,{D)f,{D^i)f,{D^i^j).J • 
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Let * = a, t+j = fiy i+j+Ic+ excluding the last 

term = fi ; and let the symbolical numerator which involves 
only direct functions be represented by/(i)), and we have 

Jf{D) 



/.(i>)/o(^-«)/.(^-i8) ...../;(i>-/*) 



Am+ot)9 



in which a, fi, ....fi,m are integers ascending by differences 
not exceeding w. 

[A few lines of the manuscript here are obscure, and I 
venture to express in other words the idea which seems to be 
involved. 

Let D — a, denote one factor oif^{D), then the correspond- 
ing factors in the denominator of 

/,(^)/.(^-«)/.(i?-i8).../.(i>-/*) ^*^' 



are 



{D — a^ {D— a^'-a) {D — a^ — fi) (5). 

Now if a, is not greater than a,, then a^ + fi is less than 
a^ + m; hence no factor in the expression (5) can be identical 
with D — m — a^, But if a, is greater than a,, then one fac- 
tor in the expression (5) may be identical with D — m — a^. 

Hence it follows that the denominator of the expression (4) 
mai/ contain D — m — a^ to the power r — 1, but not to a 
higher power.] 

And, since 

we see that u will contain r sets of terms together of the form 

^ + 5(loga:)+(7(loga?)^+ +£:(loga?rS 

-4, 5, (7.... being polynomials in x. 

This establishes the rule in [there is no reference in the 
manuscript; probably Chap. xvii. Art. 10 is to be supplied.] 
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[There is no hint in the manuscript as to the position 
which Article 2 was intended to occupy ; and the reasoning 
does not seem fully developed.] 

2. Prop. The solution of the equation 

f,{D)u+f^{D)^u+ +f,{D)e^u=^0 

l)eing expressed in the form 

{1 + F, {D) €» + F,{D) ^+ .... +} (/,(i>)}-'0, 

it is not necessary to introduce new constants in interpreting 
jF'j(i>), ...,; it suffices to interpret particularly if only uni- 
formly and consistently. 

For let 

{f,{D)ro^AP^BQ+ ; 

and in interpreting , 

F^{P)i-^(AP-[-BQ+ ) ' 

let a new constant be introduced which was not in the inter- 
pretation of 

F^^{D) 6(«-^)« (^P+ BQ-\- ). 

Now FJiip) ^J + <i>,{D) ^0 F^^{D) ^'V> 

+ ^,{D) ^F^iD) e("^»+ = 0, 

therefore 

iT, (i>) e«.» ^^ (i>) ^F^^ {D) e(-')» - . . .. 

hence the new constant comes from {^(-D)}""*0, and the term 
containing it must be -4'P, or J5'Q, ...., where A\ JS* .... are 
constants. Suppose it A'P ; 

then as derived from this, 

F^ (2)> 6<"^'" U (2>)1- = - ^1 (i>) e'^'P, • 

-<l,,(P)^»Fn.{D)^AP. 
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Thus F^^i (2>)6<-^» = -<k (2?) ^F^, e(««-)» 



The law of derivation is exactly the same as in the deri- 
vation of Fi [D) e^ from F^^ (D) €<*-')•, 

[Art, 3 seems intended for a reconstruction on an extended 
scale of part of Chapter xvii. Art. 3.] 

3. We proceed to consider more fully the theory of the 
binomial equation 

Now the possibility of solving the equation depends upon 
the nature of the symbolic function <f>{D). It is perhaps the 
most general account of the present state of the theory to say 
that there exist certain primary forms of this function which 
render the equation solvable, and that to each of these pri- 
mary forms an infinite number of the forms are reducible by 
general theorems of transformation. As these theorems 
admit of a statement which is independent of the form of the 
function ^(-0), we shall establish them first, 

■ 

Pkop. II. The function ^ {D) in the equation 

can without otherwise changing the first memher of that equa- 
tion he \st affected with any constant factor ^ or 2ndly con- 
verted into ^{I) + a), or Srdly converted into {<f) (— D)}"^. 

First. Let U=f(e^)^ and in the equation 

tt + ^(2))e«'w=/(6»), • 

- d d 

let €*=a''6*'. Then tz = 3^ > and the equation becomes 

M + a^(i>)€'*'M=/(aV), 
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ft ' 

in which D — -^,. Thas ^(i)) has been affected by a con- 
stant factor a« 

Secondly. In the same system let u = ^v. Then 

or ^v + ^(f> {D + a) e^v =/(€*), 

therefore t? + ^ (Z) + a) e*^ v = e'^^fi/) . 

Here <t>{P) has been changed into^(i> + a). 

The result of this transformation may be conveniently ex- 
pressed by the following theorem. 

The equation 

u + ^{D) e^'^u = Cr 

will he converted into 
hy the relations 

Thirdly. In the same equation let 5 = — 5' ; then 

dd" dff' 
and we have 

in which ^ = ttz? • Hence 

au 

■ tt + e-'^^{-(2)-r)}M=/(e-^); 
therefore e^u ■+ <f>{r — D)u= 6'^/(e~*0> 

whence u + {^{r-D)r€^u = {^(r-D)}-'^f(e-'). 
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In this equation let u = €*^t?. Then by the last theorem, 

v+{<f>(r i))r^v = {<}> (- i>)r/(«"*)- 

Thus 

is converted into 

f + {^ (- D)r^v = {^ (- 2))r/(€-^, 

in which D = -^ , Jy assuming 

The above transformations leave the index r in the first 

ff 
member unchanged. If however we assume ^ = — , whence 

Zfh^^'jff i '^^ should have 

By combining this with the previous results we see 
that it is possible to convert i(i?) into ^(ai) + J), and into 

But the most important transformation of the function ^ {D) 
is that which is established in the following proposition. 

[The proposition referred to is Prop. III. of Chap. xvii. 
Art. 3.] 

[Article 4 was intended to follow the words "or subse- 
quently in the derivation of w" in Chap. XVil. Art. 4,] 

4. It becomes therefore important to establish rules for 
the treatment of the constants which in these different ways 
arise. 
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Now the entire process of solution consists of three stages, 
namely : 

1st, the determination of V by the equation 

2ndl7, the solution of the transformed equation 
3rdly, the determination of w by the relation 






U — ^r 



Let us consider these separately, supposing ^ {D) to con- 
tain a single factor -jz — y which is made to disappear in the 

generation of '^{D)^ so that a and h differ by a multiple of r\ 
Thus the given equation is of the form 

u-^^{D)er*u=U (6). 

The transformed equation is of the form 

m which n = P,^^^v, F=P,-p^K 

First, suppose a — h = nr, where n is positive. 
Thus 

u=(J) + a) {D + a—r) ... [p + a — nr+r)v, 

F=|(i> + a)(i) + a-r) ... (i) + a-wr + r)| V. 



^re 



188 SYMBOLICAL METHODS. [CH. XXX. 

Hence 
w= (D + a) ... {D + a-nr + r) {l-i|r (2>)€ V^^ 

= (i) + a)...(i) + a-wr + r){l-^i|r(Z))€'^rI?7,+ Ci^^ 

+ C;€-<--')«+ + CLe"<""'"^l-.-.(7), 

where t^ is a particular value of 

j(Z)+a).... (Z) + a-wr + r)| JTl 

The part containing the constants will consist of terms of 
the form 

{D + a) ... {D + a-nr+r) {l-y}r{D)€'^}~' Ce'^'''^^ 
= (D + a) (Z) + a - r) ... (Z) + a -Tir + r) \l + y}r{D) 

= C{D + a){D + a'-r)...{D+a-nr + r) L'^^-*'^^ 

Now all these terms vanish up to the one containing 
g-(a-nr)d . thcrcforc we have to perform the operation 

C{D + a) {D + a — r) ..-. {D + a — nr + r) on 

|i|r(J))i|r(i>-r) ...i|r(Z)-yr) €-<"-"*•)« 

+ ylr{D)yfr{D-r) ... i|r(Z)-.yr -r) €-t^'*'^)^ + ....|, 

where ^ = w — i — 1 ; that is, we have to perform the operation 

C{D + a) {D + a-r)... {D + a-nr+r) on 



.-(«-*5» 



6« 
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Now ir{D) -^(Z) - r) ...ir{D --jr) e'^^"^^ 

=zy^(nr — a) y^{nr^r-' a) ... y^{nr'-jr — a) e"^*"*^ 

therefore we obtain 

= J5(7(2> + a) ... (Z) + a - wr + r) jl +ir{D) 

Thus this expression is the same in form for all values of 
t. Therefore all the terms containing an arbitrary constant 
in (7) are equivalent to only one term. 

Secondly, suppose a-^b^^-nr. 

Then w = |(i> + J) (i> + J-r)...(2> + J-wr + r)l v, 

F=(Z) + J)(Z) + i-r)... (Z)+J-nr + r) U. 

Here there are no constants in V. But u contains n arbitrary 
constants not in v, and as there is no subsequent process in 
the method for destroying these or reducing them to mutual 
dependence, it is necessary that the relations connecting them 
should be sought by comparing the solution with that given 
by the method of development in series. 

Note. It would be better to reduce (6) to the form 

before the demonstration. 

[Article 5 was intended to follow Chap. XVli. Art. 7. 

There is a memoir by Professor Boole on the subject of this 
Article, entitled On the Differential Equationa which deter- 
mine the form of the Boots of Algebraic Equations. The 
memoir occupies pages 733 — 755 of the Philosophical Trans- 
auctions for 1864,] 
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5. If we agree to regard ba primary those forms of bino- 
mial equations which are integrable but not through any 
reduction effected by the Propositions of Art. 3, and to which 
equations through the application of those propositions other 
equations are reducible and so made integrable, it becoipes 
very important to enquire what these primary integrable 
forms are. It does not appear at present possible to give a 
general answer to this question, but so far as is known, such 
forms if belonging to differential equations of a degree higher 
than the first stand in a remarkable connexion with the 
theory of algebraical equations. By the study of this theory 
Mr Harley was led to the conclusion that y defined as an 
implicit function of x by the algebraical equation 

/-ny + (n-l)a? = (8), 

n being greater than 2, satisfies the binomial differential 
equation 

in which e^ = x. In this expression the factors of the nume- 
rator are equidifferent, as of the denominator, their common 

w — 1 
difference being , but the equation is not resolvable by 

Propositions IT. and in. into forms, the integrability which 
had before been recognised. 

The above result first reached by induction was confirmed 
by Mr Cayley by the aid of Lagrange's theorem. 

To the form (8) all algebraic equations of the third, fourth, 
and fifth degrees are known to be reducible. 

Mr Harley has subsequently frftmd that y considered as a 
function of x defined implicitly by the equation 

y«-^yn-i^ (n-I)a? = a 
satisfies the symbolical differential equation 

tT' [{n - 1) Df'y - (n - 1) (nZ) - n - 1) [nD - 2Y^^y 
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the factorial notation according to which 

[m]* = in(m — 1) (m — 2) ...(w — » + l) 

being here adopted. 

These results are implicitly involved in a more general 
theorem which I shall now demonstrate. 

Theorem. If y^y^ Vn, «^^ *^^ ^ 'foots of the algebraic 

equation 

and if the rrC^ power of any one of these roots he represented 
hy u, and log a hy 0, then u as a function of 6 satisfies the dif- 
ferential equation 



w — 



\ n n \n n J 

IP? 



-e*^w = 0. 



And the complete integral of the above differential equation 
will be 

«*=C'^r + <?^r + 0„yn''. 

Let y^ = Zj then the given equation may be expressed in 

the form 

ft-i 

z = b + az'' y 
in which J = — 1. Hence, by Lagrange's theorem, 

m m n—l J m 

do 



+osf«*")*i'V*"- 



the general term being 



TTt.drF)!'^. 



192 



SYMBOLICAL METHODS. 



[CH. 



which on effecting the differentiations and adopting the £eic- 
torial notation becomes 



m 



V?* + (n — 1) r 



n 



-1 



TW 



f^i 



or 



and this expression will be found to represent the first term 
as well as tne others of Lagrange's expansion provided thiat 
we interpret the form 

[pr by 1, and [p]-' ^7 fT^ • 

Further, the above general development includes the n 

particular developments of u or y** arising from the giving 

to J* its n particular algebraic values. In this way it repre- 
sents the mr power of each of the n roots y^, y,, y» in 

succession. 

Now representing the above general term by u^^^ we shall 
have 



m 



«*r = 



m + {n — l)r 



n 



m 



'^r-n — 



"w + (w — 1) r 



r— 1 tii^-r 

6» 



n 



n 



n 



[r-nj-'' 



r—n-l tn-r 



J" 



+1 



Therefore, after reduction and replacing 6 by — 1,. 



u. 



'm+{n — l)r 



n 



- ^VM- ') 



u 



r.^n 



\rr 



(9). 



It follows therefore that the complete series of which the 
general term is w^ will if represented by u satisfy the diffe- 
rential equation 
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w — 



i)+ 1 



n n 



[i>j'» " V « = 0...(I). 



If we integrate the equation in a series (Chap. xvii. Art. 9), 
^he initial terms of the value of u will be 

- Co+ 4«+ Gji^ + G^^dr\ 

"the succeeding terms being formed from these by the law 

(9). Hence, if the arbitrary constants G^G^ C^«-i be so 

determined as to make the above initial terms agree with 
the first n terms of the Lagrangean expansion in any of its 
particular forms, the succeeding terms will also agree, and 
the Lagrangean expansion will thus become a particular inte- 
gral of the equation (I). The aggregate of such particular 
integrals, each afiected by an arbitrary constant, will therefore 
also be an integral of the difierential equation, and will, in 
fact, constitute its general integral, subject to exception only 
in the case in which for a particular value of m the integrals 
y^jy^j Vn^ cease to be independent. 

For instance, if m = — 1, and we reduce the equation to 
the form 

(y-i)»»a^-x + i = o, 

it is seen that except when n = 2, we have 

Here then the solution 

n^C\y,^+Gjf;r + Q^n" (10) 

ceases to be general for it becomes 

«=(^x-a.)yr^-h(c.-a)y,-'-+(CL.-c.)r-,, 

and virtually involves but n — 1 arbitrary constants. . 
If, however, we give to the integral the form 

B.D.E. II. 13 
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the last term of which becomes a vanishing fraction when 
w = — 1, we find for the general value of u in this case 

+ CL (y,-" log y^ + y,-" logy, +y»"'logy,), 

and in this way we may proceed in failing cases generally. 

Lastly, it may be observed that in certain cases the differen- 
tial equation (I) admits of reduction to an order lower by 
unity than its own. And in particular this happens in the 
failing cases above noticed. Thus, if in (I) we make w = — 1 
the equation will be expressible in the form 

i)(i)-l)...(i)-n + l)t^ 

^L n ^ J 

whence, operating on both members with (i) — w + 1)"*, we 
have 

The general integral of this eq^uation will be expressed by 
(10) provided that a proper relation be established between G 

and the constants G^, C^^, (7„. If we choose to determine 

G so as to give to the integral the particular form y~\ we shall 
find on substituting for u its Lagrangean development making 
7?j = — 1, J = — 1, and calculating the coefficient of a^^ or 
^«j)tf in the first member of the differential, 



^^[n_-3] 



n 



Hence, if ti be greater than 2, we have (7=0. It follows 
therefore that if n be greater than 2, the equation 
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in which ^^a has for its general integral 

y^^y^i y«_i being any n — 1 roots of the equation 

It may be useful to notice the forms which the above 
results assume when is changed into — 5, and therefore D 
into — 2> ; see Art. 3, 

It will be found that (I) becomes 



u — 



L 71 n\ \n n) 



of which the integral is therefore 

y^, y^i yn being the roots of the equation 

1 



y"--3^'"^ + l = (11); 



a 



and log a being denoted by ; 

while as the equivalent of (II) we have 



u — n 






rie^^ = (IV), 



n n 

of which, supposing n greater than 2, the integral is 

y^j^j, y,».j being any n — 1 roots of the same algebraic 

equation. 

13—2 
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Mr Harlej's results may readily be deduced from the 
above. Thus it will be found that the equation (11) re- 
duces to 

r - n< + (n - 1) a? = 
if we make 



n 



-1 / ,n4 -L (n-l)» as 



Hence, making a?= e^ and representing -j^, by D\ we have 
for the transformation of (IV) 

n 



2) = -" 

n 



h^'' 



1. -^ 
u={n^l) " e * t. 

Substituting and multiplying the result by €», we find 



n-l 

t 



-(^7 






m 



et^-^X'^rrO, 



which is Mr Harley's first equation. 



m 



If in (I) and (III) we make 1 = a, whence w = w — no, 



n 
1 



and at the same time change a into ab **, and y into y&~*, we 
-shall obtain the following somewhat more general jstatement 
of their united import. 
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The differential equations 

n-l 



1 



D-a 



»*-i 



(?--) 



^t4 = 0, 



"-3^ m 

are both satisfied hy tJi6 general integral 

tr^n ^1 J y,, ...• y» are <^e roofe o/^^Ae algebraic equation 

^-^^^"^.+ 5 = 0, 
]9rovided that for the first equation a = e*, and for the second 

' If 91 = 2, the above equations assume the forms 

f -(.-•)■ 

6. [The two principal papers hy Mr Harley on the dif- 
ferential equations exhibited on page 190 are the following : 

(1) On the Theory of the Transcendental Solution of 
Algebraic Equations, Quarterly Journal of Mathematics, 
Vol. V. pages 337... 360. 

(2) On a certain class of Linear' Differential Equations. 
Jfano^e^^r Memoirs. Third Series. Vol* il. pages 232... 245. 
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In a letter bearing date January 13, 1864, Professor Boole 
pointed out to Mr Uarley that his second equation might also 
be deduced from the general theorem discussed in Art. 5. 
Employing the above notation the deduction may be pre- 
sented in the following form. 

The equation (11) will reduce to 

r-nt"-^ + (n-l)a? = 0, 
if we make 

y = (n - l)"*aj""/, a = i (n - ifaF^ 

and for the transfonuation of (III) we have 

1 I t 

£»=-(«_ l)»c» 2> = nZ>', 



TW fH ^ 



t« = (» — 1) * € * 1*'. 

These substitutions being effected we arrire, after some 
slight reductions, at the following equation, 

n" [(n- 1) iX- ml^-^iXii'- (n- 1) [n2>' -.w-l]Vtt' = 0, 

which, making m = l and u = t, gives 

w*[(n- 1) i>'-l]"-*i)'<- (n-1) [«/y-2]V< = 0, 

an equation which admits of reduction. In fact, operating on 
both members with (2/ — 1)~\ and determining the constant, 
as in the former case, by the aid of the Lagrangean expansion, 
we find 

n"-^ [(w - 1) DY^t ~ (tlD' - w - 1) \nU - i\^eU = [« - IJ'-^e^', 

which is Mr Harley's second equation. 

The references and deduction here given were to have 
been added to the memoir which is cited in page 189, ac- 
cording to Professor Boole's desire; but by some accident 
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they were not printed, and the omission was not discovered 
until after his death. 

Mr Harley has lately succeeded in obtaining the fol- 
lowing extension of Professor Boole's theorem. 

The differential equation 



''^ [^i] 



L n ax n j {ji ax n J 
is satisfied by the wt*^ power of any root of the equation 

y- - icy""*' + a = 0, 
u being considered as a fenction of cc. 

From this he deduces the following ; the differential equa- 
tion 



w" 



I r ^ dx rj idx] ^ 



I ^\r Vn d m ^1* - 



is satisfied by the m^ power of any root of the equation 

3(* - wy""** + (w - 1) x = 0. 

For the materials of this Article I am indebted to Mr 
Harley.] 
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CHAPTER XXXI. 



THE JACOBIAN THEOKT OF THE LAST HULTIPLIER. 



1. A SYSTEM of n differential equations of iSxe first order 
and degree containing n + 1 Tariables admits of n integrals 
of the form 

Mj,Wj,, ... w„ "being independent functions of the original vari- 
ables. When w -- 1 of these integrals have been found they 
enable us to eliminate n — 1 variables, with their differentials, 
from the given system of equations, and so to obtain a single 
final differential equation of the first order between the two 
remaining variables. The final equation admits of being 
made integrable by a factor, and its solution so foimd would 
constitute the n^^ and last integral of the system. We pro- 
pose in this Chapter to develope the theory of the above 
integrating factor as established by Jacobi. The term ' prin- 
ciple of the last multiplier,' which is more usually employed, 
seems objectionable; for the essence of Jacobi's discovery 
consisted not in demonstrating the existence or the nature of 
the last integrating factor, but in the peculiar form of the 
method which he gave for its determination, and in the rela- 
tions which are implied in that form. The discovery may be 
briefly said to consist in this; viz. that instead of forming by 
means of the n — 1 known integrals the final differential 
equation between two variables and applying methods analo- 
gous to those of Chap, v., to determine its integrating factor, 
we construct antecedently to all integration a linear partial dif- 
ferential equation of the first order, any one integral of which 
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will enable us to assign an integrating factor of the final difier- 
ential equation, whatever the order of the previous tntegrationa 
may have been. Again, this partial differential equation de- 
pending for its construction only upon the form of the system 
given, we can often by examining it affirm beforehand that if 
all the integrals but one of the system be in any way found, 
the final integral will be deducible by quadratures. This 
happens in the case of the most important of all systems of 
differential equations — that of Dynamics. 

Further, an ordinary differential equation of the n^ order 
being reducible to a system of n differential equations of the 
first order, Jacobi's tneory may here also enable us to pre- 
dicate the possibility of the last integration when the previous 
integrations have been effected. 

Beginning with a single differential equation of the first 
order reduced to the form 

dx ^dy 

in which X and F are fimctions of the two variables x and y, 
we know by Chap. v. that the integrating factor fi will be 
given by the solution of the partial differential equation 

^-^-» (■>. 

the form of which should be carefully noticed. 

Consider next a system of two differential equations of the 
first order expressed in the general form 

dx^dy^dz ^ ^ 

X, F, and Z being functions of the three variables a?, y, a, and 
suppose one integral, represented by 

^(aJ, y, «)=c (3), 

to be known. The function ^ (^r, y, z), or, as we shall express 
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it for brevity, ^, will obviously satisfy the partial differential 
equation 

Xp^+Y^ + Z^ = (4), 

ax ay az 

of which indeed the given equations form the Lagrangean 
auxiliary system ; see Chap. xiv. 

If from the given integral we determine « as a function of 
X, y and c, and substitute its value in the first of the given 
differential equations, viz. 

dx ^dy 

the latter will be converted into a differential equation be- 
tween X and y. But we may leave to the equation its prior 
form, provided that we regard X and T as functions of the 
variables x and y, both explicitly as they appear therein, and 
implicitly as they are involved in z. And this being so, the 
equation (1) will become 

d{fiX) , d{iiX) dz , difiY) , difiY) dz ^ 
dx dz dx dy dz dy 

The values of -7- and -y- in this equation must be found 

ax ay 

from the known integral (3) ; they are 

dz d<f> d(l> dz ^ d<f> ^ d<f> 
rfaj"" ^x ' dz^ dy^ dy ' dz^ 

substituting which we have 

d(jiX) d<l> d(jjLX)d<l> d(fiY) d<t> d(jiY)d<l>_^ .^. 
dx dz dz dx dy dz dz dy 

This then is the partial differential equation for determin- 
ing fi. But the construction of this equation supposes <f> to 
be known. We propose to shew that fi can be determined by 
a process in which the only partial differential equation to be 
solved can be constructed without the knowledge of ^. 
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Since by actual differentiation 

^( A ^\ ^ — ( A ^^\ ^^ ^^ ^^ ^ 
dx \ dz) dz \ dx) dx dz dz dx ' 

it follows, writing fiX for -4, that 

d{fjbX) d(f> d{fiX) d<f> ^ d A y^^\ ^ ( Y^^\ 
dx dz dz dx^dxy dz) dz\ dxj' 

Similarly 

d{iiY) d4> d{ii.Y) d^^d r d4\d f^d^\ 
dy dz dz dy dy V ^/ dz \ dy) ' 

Lastly, we have 

"dzx dz) dz\ dz)' 

Now adding the. last three equations together we see that 
the first member of the result vanishes by (5) : we have thus 

-sKS)-l«)-l('^f)-- 

The second line of the first member is equal to 
and therefore vanishes by (4). There remains then 



d_ 
dx 



i^Hi'-^yu^th- 
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Hence if we put 

we have 

d(MX) ^ djMY) ^d{MZ)_^ ^,^. 

dx dy dz 

If then hy the solution of this equation a value of M dia- 
tinctfrom he founds the function -rr will he an iniegratinff 

dz 
foLCtor of that final differential equation which remains when z 
has been eliminated from the system (2) hy means of any 
known integral ^ = c. 

It will be observed that the equation for M is analogous in 
form to the equation for /i in the previous system. And this 
suggests the form of the general theorem. 

Thus proceeding to the case of a system of three equations 

dx _^dy ^dz ^dt 

we see that if 

be a known integral, yft therefore satisfying the equation 

then the system 

dx ^dy __,dz 

will virtually involve only the variables x, y, «, since t 
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-through the known integral becomes a function of a?, y, z. The 
equation (6) now becomes 

d(MX) d(MX) dt , d{MY) , d{MT) dt 
efo? ^i( dx dy dt dy 



d{MZ) d{MZ) dt^^ 
dz dt dz ' 



or putting 



dt __ <?^ ^ f?^ 
(2s; dx ' dt * 



djMX) dyfr d{MX) dyjr djMY) dyft d{MY) dyjr 
dx dt dt dx dy dt dt dy 

^d{MY) d'^ d{MY) dylt ^ 
dz dt dt dz ' 

and this is equivalent to 

i^/Kf ) /(«f ) / (^^§) 

dx dy dz dt 

and therefore becomes on rejecting the term in the second 
line by (7), and putting 

djNX) d(lTr) d(NZ) iVT) _ ., 
dx dy dz dt "* ^ ^" 
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If from this equation a Tslue of N distinct from be ob- 

N 
tained, then if=^, and therefore 

N 



dyjt d(l> 
1U da 



This is the final multiplier, i. e. the integrating factor of 
the final difierential equation between x and y which remains 
when z and t have been eliminated from the given system by 
means of the two known integrals. In calculating fi from 
the above formula we must proceed as follows. The value of 

-^ must be found from any given integral -^ = c ; but that of 

•=2 must be found from another integral from which by means 

of the former one t has been eliminated. Thus the general 
forms of the integrals will be 

^ (a?, y, z, c) = c'. 

Lastly, the values ^^-jt j -^ found as above, and that of 

N given by any solution (distinct from 0) of the partial dif- 
ferential equation (8) having been substituted in the expres- 
sion for fly we must eliminate z and t from that expression by 
means of the two known integrals. The resulting function 
of cc, y, c and d will be the integrating factor sought. 

The reasoning above employed is in its nature quite inde- 
pendent of the number of tlie equations of the original sys- 
tem. The general theorem to which it leads may be thus 
stated. 

Theorem. The system of n differential equations 

^ Y^ ya "•••!; 
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being given, if a system of « — I integrals 

be so reduced by elimination that the variable y^ shall not 
appear in <^g, the variables y^, y^ shall not appear in ^3, and 
so on, then the integrating factor fi of that final differential 
equation between x and y„ will be given by the formula 

M 






in which Jf represents any integral distinct from of the par- 
tial differential equation 

d(MX) djMY,) .d{MY^)_. 
dx dy^ dy^ 

In applying this theorem the expression for /jl must be 
freed from all the variables except x and y„ by means of the 
given integrals. 

This is Jacobi's theorem. On account of its great importance 
I propose to give another demonstration of it founded upon 
the Cfalculus of Variations. 

2. Second demonstration founded upon the Calculus of 
Variations, 

It will be most convenient to present the proposed system 
of differential equations under the symmetrical form 

dx^ dx^ __ dx^ 

the independent variables being x^ x^, a?, of which 

Xj, Xj, .... X^ are any functions. We have thus n — 1 diffe- 
rential equations, and we are to seek the integrating factor of 
the differential equation which remains when by means of 
» — 2 known integrals w — 2 of the variables with their diffe- 
rentials have been eliminated. 
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Suppose P= c to be any integral of the system, then JP 
satisfies, and it suffices that it satisfies, the partial differential 
equation 






(a). 



Now if in place of a?j, a?,, ....a?» we introduce a new sys- 
tem of independent variables t^j, t^, •••• t««, which are functions 
of the former, then we shall have 

^dP^^dP ^yr dP 

^jdP ^ ^dP rr dP 

ZJp Z^, .... U^ being functions of u^^ Wj, ....w^. And by the 
theory of the transformation of multiple integrals. 



where 



H= 



du^ 



du^ 
dx^' 



du^ 

dXn 

dx^ 



The foregoing equation we may express in the form 
SI Xi'-T' dXjdx^..*dx^=^ I -A-r-'du^du^..:du^. 



Hence, representing by 8 an operation of differentiation 
which affects only the form of P as a function of 0;^, arj, ... a?^ 
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or of u^, Ui, .... w„, and not the independent variables them- 
selves, we have 



S I Xi -1 — dx^ dx^ .... dx^ = 2 I -^F ~7 — du^ du^ . . . du^^ 

and therefore integrating by parts and e(]^uating the portions 
on each side which remsdn under the sign of n-fotd inte- 
gration, 

S I V-^ SP dx^ dx2 .... dx^ 

= 2 I ~r-(-M8Pdu^du^...,dUn. 

Whence again transforming the integral in the first member 
^ f" dXi jjp du^ du^ . . .. du^ 

and this being true quite irrespectively of the form of P, 
we have 



I2— *=2— f-*V 

H ^ dxi dui \HJ * 



In this equation Jacobi's theorem is virtually contained. 
For let the given equation be multiplied by any factor. Then 
changing in the above Xi into -MZi, and Ui into MUi, we 
have 

1 ^ d{MX,) ^^ a /MU,\ 
S dxi dui \ H J * 

Hence, if Jf be determined to satisfy the equation 

axi 

B.D.E. II. 14 
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we shall have 

si(^=» (»)•■ 

This is wholly independent of the relations connecting 
t*j, u^y...u^ with iTj, a?,,....a?„. Now choose the n — 2 variables 
Wp w,, .... w^ so that Wi=^Cj, W2 = ^2> ••••^»»-s~^«-2 shall be 
integrals of tne given partial difi'erential equation (a). Then 
that equation transformed becomes 

of which the auxiliary ordinary equation is 

^ndu^^- U^du^=0. 

At the same time the equation (J) becomes 
d fM ^ \ d (M 



du^^ 



(l^")-^.(5^.) = ''- 



if . . . 

Hence -jf is the integrating factor of the preceding diffe- 
rential equation between u^^ and u^^ 

Jacobi's theorem in its most general form is thus seen to 
be the following 

Theorem. If the system of differential equations 

dx^ _ dx^ _ dxn 



be transformed by the introduction of a new system of vari- 
ables M^, «*3, ....w„, so chosen that 

shall be integrals of the given system, then the final differen- 
tial equation between u^^ and u^ shall have for its integrating 
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factor -g:, in which M Is any function satisfying the partial 
differential equation 

djMX,) djMX;^ d{MX:) _ 

dx, ^ dx^ •••• + — tZ ", 

and JS" stands for the determinant 

du. 



dx. 



dx^'' 



du^ 
dx^ 



da^ 



du^ 
dx. 



The form of Jacobi's theorem obtained by the previous 
demonstration may be deduced from the above by choosing 
for u,^,u^ two of the original variables, for example a?^^, a?^, 
and transforming the integrals w^, w^, ..„u,^ so that u^ shall 
contain only x^..,x^^u^ shall contain only x^.,.x^j and so on. 



Examples, 

3. Jacobi has established by means of the above theorem 
the very remarkable theorem that in^ any ordinary dynamical 
problem the forces depending not iipon the time but upon the 
material constitution of the system, if all the integrals but 
two of the dynamical equations are found, the two remaining 
integrals can be found by quadratures, 

1st. In a dynamical system of free points the forces act- 
ing upon which depend only upon the position of the points, 
we have if we represent the entire system of rectangular co- 
ordinates taken in any order by a;, y^z,... and the correspond- 
ing resolved forces divided each by the corresponding mass 
by X, Yj Zf... the system of equations 



d^ 



df 



14-2 
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or putting 

dx _ , 3,y _ , 



• • • 



, da? _ e7y _ dx* _ dy 
X y X 1 

Now as X, F... do not contain i we may consider first the 
system 

dx __dy _ dx' __ dy' 

and it is evident that if we can find all the integrals of this 
system, t will be given by the equation 



{dx 



x' having been first converted by means of the supposed in- 
tegrals into a function of x. 

To determine the last multiplier of the system last writt^ 
we have first the equation 

d jMx') d (My') djMX) d{MY) ^^ 

dx dy "" daf dy' '" ' 

which since X, Y... do not contain x\y\.. is satisfied by 
Jf = a constant Giving to the constant the particular value 1, 
we see that if 

are n — 2 integrals of the system, and if by means of these 
we eliminate n—2oi the variables and construct the differen- 
tial equation between the two remaining variables, the inte- 
grating factor of that equation will be -^, in which -ff is the 
functional determinant of i^^, w^, .... w^. 
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2ndly. Suppose the system subject to a material connex- 
ion which establishes an equation of condition among some 
or all of the co-ordinates. If we represent the co-ordinates 
taken in any order and multiplied each by the square root of 
the corresponding mass by a:, y,... the corresponding resolved 
forces by X, F, ... and the equation of condition expressed by 
means of the above modified co-ordinates by = 0, the diffe- 
rential equations will be 

the transformation above employed reducing all the equations 
to the same type. [See the next Chapter.] 

Making 

dx __ , dy _ , 

di^^' di^^'"' 
the system becomes 

dx dy 

and the Jacobian equation for M becomes 

— TZ T 7 •••• + 




dx dy '"* d^ 






dy 

Now ^ does not contain a?', y'.... Let us inquire whether 
it is possible to determine il also as a function of x, y .... 
without x'y y.... so as to satisfy the above differential equa- 
tion. 
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The equation would become 

dx ^ dy "* \dx dx dy' dy "'*J^ > 

or if we write 

' ^ A. * A. — ;^ 
^^^y dy ^' 

and from this we must eliminate X. 

Now since = 0, we have by differentiating and putting 

dx , 

and again differentiating 

d^d^ d^dy' _ 
^dx dt "^ dy di'^ ' 



or since 



dx' y-.'.d'k 

'^ dix^^y df^ — ^^^^dxdy^"" 
dx dy 



+ X 



{(l)'-©"--}=«- 
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and differentiating with respect to x\ 



or if we make 



Similarly 



ax dx 



^'f*l«=«- 



Therefore 



d(f> ^ d<f) d(l) ^d(f> 
ax ax ay ay 

fd\ d(f> dKd(l> \n-n 

^Kdx'di^dy'dy'^'''')'^^^' 

or 

. , T . ^. dX dd> , dX dd> , 

and now elimmatmg j~f 'f^ + 'j~i "T"^ '*' 

we obtain QhM— MZ Q = 0, 

which is satisfied "by Jf = Q» , 

4 [Among Professor Boole's manuscripts I found five 
pages in German, forming part of a memoir, which was pro- 
bably intended foF Crelle s Mathematical Journal. The 
memoir was to have discussed two applioations of the Calcu^ 
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lus of Variations; one to the Jacobian Theory of the Last 
Multiplier, and the other to the Solution of Pfaff's equation 

X^dxi + X^dx^+ +X^dx^=0. 

But there is only a single paragraph relating to the second 
application. 

The manuscript contains the same demonstration of the 
Jacobian Theory of the Last Multiplier as in Art. 2 of the 
present Chapter; after this demonstration some remarks occur 
of which the substance wiU now be given.] 

It is worthy of notice, that Jacobi in the 36th volume of 
Crelle's Journal^ deduced by the aid of the Calculus of Varia- 
tions the result on which the preceding demonstration of the 
Theory of the Last Multiplier depends. Li fact, he shewed 
that if V denotes any function of 



*^l' *^8' ••" •^•»' ^ 



dz 



»' "' dx,' 



••• 



dz 
dx^^ 



and Fbe transformed by the introduction of a new system of 
independent variables w^, «*,, ... u^, then the following rela- 
tion holds. 



'dV d dV 



dz dx^ , dz 
dx^ 



d dV\ 

dx,^d^\ 

dx^f 



^ c?(AF) d d{LV) d d (A F) 

dz du^ 7 dz 



du. 



du^ jdz 
d -r- 
du^ 



where 



A = 



dx^ 
du^ ' 



dx^ 



dx, 
du^ 



dx^ 
du. 
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Jacobi applies this result to the transformation of the ex- 
pression 

d^V d'V d'V 

dx' ^ d'^'^ de' 

4 

But neither Jacobi himself, nor any x)ther person, so far as I 
know, has drawn attention to the application of the result 
which I have given here. 

[The substance of the single para^aph relating to the 
second application of the Calculus of Variations will now be 
given.] 

Clebsch has earned the thanks of all who are interested in 
the higher parts of the Theory of Differential Equations, since 
he has performed the same service for Pfaff's problem as 
Jacobi did for the Theory of Partial Differential Equations of 
the first order, and thereby for the equations of Dynamics. 
But while I recognise the great importance of the results, I 
consider it desirable to give a simpler deduction of the system 
of partial differential equations therein involved, and on which 
the other results depend. 
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THE DIFFERENTIAL EQUATIONS OF DYNAMICS. 



[It will be seen that this is only a fragment of the Chapter 
which was to have appeared under this title.] 

I do not propose in this Chapter to discuss the origin and 
interpretation of the differential equations of motion or to enter 
into those details of their application which are found in all or- 
dinary treatises on Dynamics. But they constitute a system 
analytically so remarkable from the forms in which it is 
capable of being expressed, and from the general methods of 
integration which emerge out of those forms, that they are 
well deserving of a special attention. 

Referred to rectangular co-ordinates the differential equa- 
tions for the motion of a system of points free or connected 
are 

d^x' _ y d(f> d^t 



THE DIFFERENTIAL EQUATIONS OF DYNAMICS. 219 

Here m is the mass at the point (a?, y, 2?), tri that at (a?', y\ z'), 
Xy Yy Z the resolved forces at (a:, y, z) tending severally to 
increase those co-ordinates, and so on. Lastly 

^=0, '^ — 0,... 

are the equations of condition each of which may involve all 
the co-ordinates, and X, /^... are indeterminate multipliers. 

The above is usually termed the first Lagrangean form of the 
differential equations. In applying it we must either elimi- 
nate \, /^... from the given equations, and then by the equa- 
tions of condition just so many of the co-ordinates with their 
differentials, or we must retain X, /^, ... as variables so conditioned 

that the values of -^ , -;i^,-«« in the system shall satisfy iden- 

tically the differential equations involving -^, -^v* de- 
rived from ^ = 0, -^ = 0,... viz, the equations 

^ = ^ = 

de^' df ^'••' 

The first Lagrangean system may by a slight transforma- 
tion be reduced to a form in which all the equations are of 
one type, viz. of the type which they would have if all the 
masses were equal to unity. 

For taking the first equation of the system and dividing 
by m* we may express the result in the form 

d^{rrfix) X ^ d^ dyfr 



= — 1+^ ,. 1 . +/^ 



• •• 



^^ m* d{m^x) d{mSx) 

from which we see that if a;, y ... had been taken to represent 
the entire system of co-ordinates taken in any order and mul- 
tiplied each by the square root of the corresponding mass, and 
Xy F... the corresponding resolved forces taken in the same 
order and divided each by the square root of the correspond- 
ing mass, th6 system of equations would have been 
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all being of one type. In general investigations this form 
is to be preferred. 

From the first Lagrangean form another known as the 
second Lagrangean, and from this again a third known as 
the Hamiltonian are derived. The second Lagrangean form 
is properly speaking an expression for the effect of a trans- 
formation of co-ordmates in the most general sense upon the 
original system, i.e. of a transformation which in place of 
0?, y, ... me entire system of given co-ordinates substitutes 
a new system of variables f, ^, ... the expressions of which as 
functions of a?, y, ... are known. It is not necessary that this 
new system of variables should be co-ordinates in the proper 
sense of that term, determining three by three the positions 
of the several masses; it suffices that they should in their en- 
tirety determine and be determined by the co-ordinates given. 

The second Lagrangean form may be established as 
follows : 

Differentiating the equations <^ = 0, '^ = 0,. . . with respect to 
any one of the new variables f we have 

dj>dx d<f>dy __ 
dx d^ dy d^ *" "" ' 

d^jrdx d^dy _a. 
dx d^ dy d^'" ' 

whence if we multiply the equations of the given system by 

dx dy in 1 

-7^ , -^% ,... and add, we have 
d^' d^' 

dx d^x dy d^y ^ y.dx y^ 

dlW^dlW""^^dl'^^dl'" 
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CHAPTER XXXIIL 



ON tHE PEOJECTION OP A SUEFACB ON A PLANE. 



[The following memoir was found among Professor Boole's 
manuscripts; a Title and Introductory Remarks were to have 
been prefixed, but with this exception the memoir appears to 
be finished for publication. It is sufficiently connected with 
the subject of Differential Equations to find a place in the 
present volume. 

The memoir by Sir John Herschel to which allusion is 
made is entitled, On a new Projection of the Sphere ; this was 
read before the Royal Geographical Society of London on 
the 11th of April, 1859, and was printed as part of the Journal 
of the Society, Vol. xxx. 1860, pages 100... 106. A chart of 
the World on Sir John flerschers projection has been pub- 
lished by A. and C. Black of Edinburgh. 

The history of the subject will be found in Chapter xxiil. 
of the Coup d'ceil historique sur la Projection des Cartes de 
OSoffraphie... Par M. D'Avezac, Paris, 1863. 

For the materials of this introductory notice I am indebted 
to Sir John Herschel.] 

1. Let aj, y, z be the rectangular co-ordinates of any point 
on the given surface; x', y' the co-ordinates of the correspond- 
ing point on the plane of projection. Let the equation of the 
given surface be 

or, for simplicity, 
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The condition of projection upon which Sir John Herschel's 
investigations are founded, and which we shall adopt here, is 
that of the similarity of corresponding infinitesimal areas on 
the surface and on the plane. The object of the problem then 
in general is the discovery of the mode in which a;', %/ depend 
upon a?, y, and z in accordance with the above condition; its 
object in any particular case is the determination of a?', y as 
functions of a;, y, z. 

Regarding then a?', y' as ultimately functions of x, y^ z we 
have 

ax =^-7- dx + '—-du + ^r dz, 
ax dy ^ dz 

in which dx, dy^ dz are not independent, but are connected 
by the condition 

dFj dF ^ dF ^ ^ 
— d^+dy+-r-dz = 0. 
tjujLf ^y ^^ 

Now for brevity write 

ddd __ doi ^ - dal __ 
dx^ ^ dy"^ dz ^ ' 

^'=a' ^' = J' ^' = c' 
dx ^ dy ^ dz ' 

dx ^ dy ^ dz ' 
then 

dx' = adx + hdy + cdz (1), 

dy =adx + l'dy + c'dz (2), 

0=^Adx + Bdy+ Cdz (3). 
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Now the condition of. the similarity of infinitesimal cor- 
responding areas may be resolved into the two following 
conditions, viz.: 

1st. The equality of their corresponding angles. 

2ndly. The proportionality of their corresponding sides. 

And these conditions we shall introduce separately. 

1st. Assuming any point x\ 'if on the plane of projection, 
let a?' alone vary, and the infinitesimal line generated is dx^ 
while (since dy' = 0) (2) and (3) become 

a*dx + h'dy + c'tfo = 0, 

Adx + Bdy+Cds^O, 

whence, if we write 

L^Bc-CV, M^Ca'-Ac', N=AV-Ba\ 

. dx dy dz .. 

we have T^~M^N ^^' 

so that the direction cosines of the infinitesimal line on the 
surface F corresponding to the line dx' on the plane {x% y[) 
will be 

(i'+if'+JV^)*' {p^M'^N'f' (i^+jf^+iv^)* ^ ^" 

In like manner, if y alone vary, we shall find for the 
direction cosines of the infinitesimal line on the surface F 
which corresponds to dy' on the plane 

{L-'+M^+N'")^' (,L"'+M"+2r')i' {L'^+M'+N")^ ^ '' 

where L'=Bc -Cb.M'== Ga-Ac, N' = Ah -Ba. 

By the first of the conditions of similarity the fingle be- 
tween these lines on the surface must be a right angle since 
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daa and dy' aie at right angles. Hence we have, &om 
(5) and (6), 

LL' + MM' + NN' = (i (7). 

2ndl7. The ratio of the length of the element did to the 
corresponding element on the simace is 

dx' 



Jda^+^ + dz'' 
or, by (1), 

adx+idy + cdz 
'j^ + df+di' 

and therefore by (4) 

aL + hM+cN 



'JL' + M^ + N" 

equating which to the corresponding expression for the ratio 
of the length of dj/ to that of its projection on the sui&ce, 
we hare 

aL + hM+ cJSr ^ a'L' + b'M" + c'lST . . 

'^IF+W+N'*~ 'fU^TM^TW 

Now if we substitute for L, M, N, L', M\ jV' their values, 
we shall find 

aL + lM+cN=A{h'c-lc') + B{c'a-ca') + C{a'b-ab'), 

a'L'+ h'M'+ c'N' = A {be'- b'c) + B {ca!- c'a) + C {ab'- a'b), 

and the second members of these equations differ only in sign. 

Thus (8) may be expressed in the form 

[a {b'c - be') + B {da -ca')+C {ab'- o'5)l 

X I \ T ^ t}- = 0...(9). 
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But tte first factor of the first member of this equation 
being the determinant of the system 

adx + hdy + cdz = 0, 

a e&j + t't?y + c tfo = 0, 

Adx+Bdy + Gdz = 0, 

expresses when equated to zero the condition that if in the 
system (1), (2), (3) dy* vanishes dx* shall also vanish; and 
da! and dy being independent, this condition cannot be satis- 
fied, so that (9) reduces to 



= 0, 



whence 

r^j^M^j^ir^U^jSP^lP^^ (10), 

and this, with (7), will fully express' the conditions of simi- 
larity. 

2. If we multiply (7) by 2 V— 1, and add and subtract 
the result £rom (10), we obtain the equivalent system 






(i' + ii V- !)• + (if' + if V^^i)' + {ir + ^V-i)' = 

(i' - L V^)' + {¥"- M4^y + {W- N-Z^l) 
Now i'±iVri = ^^'_^^' 

ay dz dz dy 



(dF d£ _dF df\ J— 
\Sy dz dz dy) 



dF d{x[±yW--l) dF d{x' ±y'^r^) 
dy dz dz dy * 

a D. E. II. 16 
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Writing thdn 



we have 



V^-L^/^l^ 



dFdu ^dFdM 
dy dz dz dy^ 



^ j—T dFdv dFdv 



In the Bame wajr 



M'^M^P^^ 



dF du _^dF du 
dz dx €^ 3S^' 



Jlf'-.JfV^ = 



dFdv ^dFdv 



nr + Nj^^ 



dF du ^F du 
dx dy dy dx^ 



XTf XT /^-T <iF dv dF di) 



Substituting which in the system (11) there result 



fdFdu^ 
\dy dz 



dFduV /dFdu 
dz dy) \dz dx 

,(dF da 
\dx dy 



dFdu 
dx ^ 

dF du 



1 



"N 



/OF do _d£ 
\dy dz dz 



dy dz 



dFdvV 
dz dy) 



(dF do 
\dz dx 



.(dFdv 



\dx ^ dy 



^d^'duV^ 
dy dx) 

_dFdoy 
dx dz) 

dFdoV 
dx) 







^....(12), 



= 
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to which we may give the somewhat more convenient form 

- {(S)'-(f)'^(S]0-(|)'-^©} 

\dx dx dt/ d[y dz ^J ***"^ '* 

\dx dx dy dy dz dzj~^ *' 

These are partial differential equations of the first order» 
serving to determine u and v as functions of a?, y, z. 

But it is not necessary to solve the equations in their 
general form. For, x, y, and z being connected by the equa- 
tion of the surface, the above e(]^uation3 may always be so 
reduced as to involve only two independent variables. As 
latitude and longitude determine the position of a point on 
the earth, so two co-ordinates of any given species will deter- 
mine the position of a point on the given surface, and these 
co-ordinates, when fixed upon, become the independent varia- 
bles of the problem. 

Let 8 and t represent such co-ordinates, and let their ex- 
pressions in terms of x, y, z give 

« = ^1 (»» y> «)» < = ^« (»> y> «)j 

which equations combined with that of the given surface will 
reciprocally determine x, y, z as functions of s and t Then 
1st the differential coefficients of F which in the equations 
(I), (II), are functions of a?, y, z may be transformed into func- 
tions of 8 and t ; 2ndly, we have 

du du ds du dt 
dx d8 dx dt dx^ 

du ^du da du dt 
dy ds dy dt dy^ 

du __ du ds du dt 
dz ds dz dt dz ' 
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and as ;7- > -r-.** are known functions of a?, y, z^ thej also 

are expressible in terms of s and t. The result of these sub- 
stitutions will then be to convert (I) into a partial differential 
equation in which u is the dependent and s and t the inde- 
mdent variables, and this equation being, like (I), of the 
rst order and second degree in the differential coefficients of 
w, will be of the form 

For V we shall have an exactly similar equation with the 
same coefficients. 

The above equation is, by the solution of a quadratic, 
resolvable into two equations of the form 

du •). du _^ du •. du 

dl^^'dt'^' d^^^'dt"^' 

To these correspond the respective auxiliary equations 

dt + \ds=0, dt-\'\ds = (13). 

If the integrals of these are 
respectively, then we have 

Now V being determinable by an equation of the same 
form as w, it follows that of the above two values of u one 
must be assigned to v, so that the solution of the problem will 
be contained in the system 

or in the system 

u = <l>{T), v = y{r{S). 
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The particular forms of the arbitrary functions <j} and y^ 
will depend solely upon the nature of the problem under con- 
sideration. 

• 

One other point remains to be noticed. The first mem- 
bers of (12) are essentially positive, being composed of 
squares ; so are then the first members of (I) , (II) ; and so, 
if the intermediate transformations are real, is the first 
member of the equation whose coefficients are P, Q, R, 
Hence the quadratic determining X^, \ will have imaginary 

roots of the form a + )8 J— 1. Ultimately therefore it will suf- 
fice to integrate one equation of the system (13) and then to 

deduce the solution of the other by changing ^A-l into 

3. Application of the above formulce when the given sur- 
-face is an oblate spheroid, such as the earth. 

Let the plane of the equator be that of projection, the 
centre being the origin. Let the co-ordinates x, y pass 
through the meridians of and of 90° respectively, and z 
through the poles. The equation of the surface will be 

t±t.t^X (14) 

where a is the earth's equatorial, b its polar radius. ' Let also 
the latitude of the point a;, y, z be represented by «, the 
longitude by t. We have 

dFjix dF_2y dF_2z 

dx~c^' 'dy~ €?' dz~'W' 

and substituting in (I), 



'aj^ + y 



v'^mK) 



'<t 



X du y du z du^ __ 

^jj^ dx a? dy b^ dz) ' 
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or, if we represent ^ by A", 



(-•-^■-»'^{(l)"-(|)'-(s)] 

\ die ^ dy dz) 



(15). 



Now BS 3P, y are rectangular co-ordinates in the plane of 
the equator^ and x passes through the first meridian, we have 



^ = 



X 



iamt 



Again, representing in the annexed figure the meridian of 
the point P, or {x, y, z) 
touched by the straight line 
QB in the same plane, we 

have g.M"=V ^4y, MP=^z. 

Therefore if Vaj' + y^=r, the 
equation of the meridian is 

r^ z^ 
that of the tangent 

r\ z being current rectangular co-ordinates of the tangent. 
Hence 




c^z 



lez 



tan 0^^ = TT- = -7= 
But GQR = latitude. Therefore finally 



8 = tan"^ 



fi'z 



Vi'"+ 



1/' 



, « = tan"^^ 

X 



(16), 



and we must now transform (15) so as to make s and t the 
independent variables. 
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From the above equations combined with (14) we find 

__ ah cost _ ahaint atan^ 

^""VA« + tan«? ^"VFT^I' ^^h^W+t^s'''^ ^' 

and substituting in (15), 



Beck's 



( 'du , . du ^ duV ^ ,,^. 

— fcos^^ + sm^^ + tan^^j =0 (18). 



Again, 



du ^du da du dt 
dx da dx dt dx^ 

du du da du dt 
dy da dy dt dy^ 

du ^du da du dt 
dz da dz dt dz' 

jT da _ ' — h^zx — sin a cos a cos t JM 

^">/5*T/(aj'+/ + AV) oJT 

where H^ A* + tan"^. In like manner 

da _ — sin g cos a sin tjH 



dy 


ah 


da 

dz 


h cos'* JM 
a 


dt 
dx'^ 


— sintjM 
ah 


dt 
dy'^ 


cos tsfS 

ah ' 


dt__ 
dz"^ 


0. 
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Hence 

du JH I . ^ A'^ * .^^\ 

-71- = ^S- — sins C0S5 c^o^t-z — sin t-j- , 
dx ah \ ds dtj 

du Jli f . . ^du , du\ 

^i- = ^^S- — 8in« cos 5 sin^-T- +cos<-y- , 
dy ah \ ds dtJ 

dz ah \ day 

Substituting these values in (18), and dividing by the com- 

JT 

mon factor -^^^ we have on reduction 
an , 

(|)V co8«« {1 + (A«- 1) coaM'(gy=0, 
which is resolvable into 

^ - ;/:=:i coss {1 + (A*- 1) cos'«} ^ =o, 

^ + 7=nrcos5{l + (A«-l)co3^}^=^0, 

partial differential equations of which the integrals are in- 
cluded in the common formula 

Now I 



C03 «{1 + {h' — 1) COS**} 

_ f ds , /•. _ it\ [ C08* <^ 
~ J COS* ■•" ^ '^' j 1 + (A»-l) cos's 

Jcos* ^ 'J A'— (A* — 1)81 



sm's 
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rs / e^lz a » i t Since e=s — »— I 



— 6 sin d 
sins 



» ^ /TT* «\ . e, 1-6S1 



Hence 



»-*h{(ri5i-:/-(f^i)l*'^l- 



or, changing ^ (^) into ^ (e*), 



'-*{(TiS;)''»(!+l)-^}- 

Let r and 9 be the polar co-ordinates of that point in the 

Elane of projection which corresponds to the point whose 
ititude and longitude on the surface are s and t ; and let 

U + esin^y U 5^/ 
then the complete solution as£iumes the very simple form 
r^^^^ ^ (56*«>^, r6-«>/=^ = t (/S€^«^ (III). 

Of particular deductions the most interesting is that which 
atises from the supposition that the parallels of latitude are 
projected into circles round the pole. This requires that r 

B.D.E. lu 16 
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should be independent of t, a condition whicli is Batisfied in 
the most general manner by assuming 

we then find 

whence, on multiplication and division, 

u 

whence, A and S being new arbitrary constants derived from 
C and C 

If we observe that 9 and t should vanish together, we have 
-B = 0, and the equation 0= ±nt shews that the surface of 
the sphere will be projected into a sector of a circle, the arc 
of which is to the circumference of the circle as n : !• Thus, 

if w = 2 , the sphere is projected upon a quadrant, and so 

on. 

The other equation gives 

If 5 = we find r = -4, whence A is the distance of the 
equator from the pole in the plane of projection, and if that 
distance, which is arbitrary, be assumed as the unit, we have 



= tan |f^ + 1)17 J^l^^)^ 
\\4. 2j} \l-\-eam8j 



for the distance from the pole of that parallel whose latitude 
is 6. We may give to this expression a better form by 
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assuming i? = ^ + «, and introducing an auxiliary quantity q 
determined by the equation 

ecosp = cosj. 
We have then 

The following table gives the values of r for the sphere and 
for*the spheroid whose eccentricity is '08 (which is about that 
of the earth), for each ten degrees of polar distance, for the 

values w = 1, and n = - . 

4. 



Polar 


n- 


= 1 


n 


1 

"4 


Distance. 










Sphere. 


Spheroid. 


Sphereii 


Spheroid. 


10<> 


•0875 


•0880 


•5439 


•5447 


20® 


•1763 


•1774 


•6480 


•6490 


300 


•2679 


•2694 


•7195 


•7205 


40^ 


•3640 


•3658 


•7767 


•7777 


60^ 


•4663 


•4682 


• -8264 


•8272 


60^ 


•5774 


•5792 


•8717 


•8724 


w 


•7002 


•7017 


•9148 


•9153 


S(fl 


•8391 ' 


•8400 


•9571 


•9574 


90<> 


1-0000 


1-0000 


1-0000 


1-0000 


100» 


11918 


1-1904 


1-0448 


1-0445 


IIQO 


1-4281 


1-4250 


1-0932 


1-0926 


1200 


1-7321 


1-7265 


1-1472 


1-1463 


1300 


21445 


21357 


1-2101 


1-2089 


1400 


2-7475 


2-7340 


1-2875 


1-2859 


1500 


3-7321 


3-7114 


1-3899 


1-3880 


1600 


5-6713 


5-6372 


1-5432 


1-5409 


1700 


11-4301 


11-3581 


1-8387 


1-8358 



THE END. 
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